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ABSOLUTE YIELDS OF THE XENON AND KRYPTON ISOTOPES 
IN U* SPONTANEOUS FISSION! 


B. G. YounGc? anp H. G. THODE 


ABSTRACT 


The absolute abundances of the isotopes of fission-product xenon and krypton 
in six uranium minerals have been determined mass spectrometrically using the 
isotope dilution technique. The fission products were resolved into a U** spon- 
taneous fission component, a U2* neutron-induced fission component, and a U8 
neutron-induced fission component. Internal consistency in the analysis was 
achieved only when the Xe!® yield used for the U** Cherm neutron fission 
component was 20% lower than that reported by Purkayastha and Martin at 
]!29. This discrepancy in the mass 129 chain yield measured at |! and at Xe!29 


has not been resolved. 

Only one of the six minerals, Cinch Lake pitchblende, retained essentially 
all of its fission product inert gases throughout geological time. Inert gas losses 
from the remaining five minerals ranged from 20% to 75%. The absolute yields 
of the stable xenon and krypton fission products in U** spontaneous fission were 
determined from an analysis of the inert gases from the Cinch Lake mineral. 
This analysis showed that 95.9% of the fission gas in this case resulted from the 
spontaneous fission and only 4.1% from neutron-induced fission. 


INTRODUCTION 


Soon after the discovery of fission, Petrzhak and Flerov (1940), Maurer 
and Pose (1943), Scharff-Goldhaber and Klaiber (1946) demonstrated with 
ionization chambers and boron trifluoride counter tubes that uranium under- 
goes fission spontaneously and emits neutrons in the process. Later, Segré 
(1952) used enriched isotopes of uranium and showed that U** has a spon- 
taneous fission rate of about 23 times that of U**, the half lives for spontaneous 
fission of U*** and U** being 8.04 X 10" years and 1.87 X10" years respectively. 

In 1950, Macnamara and Thode, using mass spectrometric methods, identi- 
fied xenon and krypton fission products in the rare gases extracted from old 
uranium minerals. The pattern of xenon and krypton isotopes found for 
natural fission indicated an asymmetric mass yield curve similar to that obtained 
in the neutron fission of U** and U8. There were, however, certain marked 
differences. First, the mass yield curves in the xenon and krypton mass ranges 
are steeper for natural fission. Since there is much less energy available for 
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spontaneous fission than for neutron fission, the process is expected to be more 
selective (steeper curves). This was considered strong evidence that spontaneous 
fission of U** is involved in natural fission. Secondly, the abnormally high 
yield at masses 133 and 134 (fine structure in the mass yield curve) found for 
neutron fission of U?* is shifted toward the lower mass for natural fission, 
which is the same direction of shift found between neutron fission of U2 
and neutron fission of U*, Finally, the xenon-to-krypton ratio was found to 
be much higher in the natural fission gases. The low-mass hump of the mass 
yield curve is shifted toward higher mass with increasing mass of the fissioning 
nucleus. This results in lower krypton yields in the fission of U?** as compared 
with U*5 fission. From this and other evidence, Macnamara and Thode con- 
cluded that a large part of the fission gases in uranium minerals resulted from 
the spontaneous fission of U%, 

In 1953, Fleming and Thode showed that both neutron fission of U2 and 
spontaneous fission of U** occur in nature. They found that the fission isotope 
patterns of xenon and krypton varied from one uranium ore to another, 
depending on the proportions of the two fission processes which have taken 
place. Wetherill (1953) confirmed the existence of neutron and spontaneous 
fission in uranium ores. Since the mass 129 chain yield is about 1% for 
neutron fission of U*** and was found to be very nearly zero for spontaneous 
fission of U**, the yield of this isotope was taken as one index of the propor- 
tion of neutron fission of U*** and spontaneous fission of U8 which has occurred 
in a given mineral. The fission yield of Xe'?* in several minerals was less than 
0.02% which indicated that at least 98% of the fission gas was due to the 
spontaneous fission of U**. Other samples showed up to 23% neutron fission 
of U*, There are two main sources of neutrons for this process, spontaneous 
fission and (a,m) reactions on light elements. 

Recently the absolute yields of Mo*® and several iodine isotopes from U?% 
spontaneous fission were reported by Parker and Kuroda (1958) and Ashizawa 
and Kuroda (1957). 

In the present work, absolute abundances of the cumulative yields of xenon 
and krypton fission products from six uranium minerals were determined mass 
spectrometrically using isotope dilution techniques. The purpose of this work 
was to obtain the absolute yields of the stable xenon and krypton fission pro- 
ducts ior U*** spontaneous fission, to determine the leakage of fission gas from 
the minerals throughout geological time, to assess the spontaneous fission 
method of mineral age determination, and, finally, to investigate further the 
various fission processes which occur in nature. 


EXPERIMENTAL 


The uranium minerals were crushed and sieved to provide a homogeneous 
aggregate for the extraction and analysis of the fission product gases. The 
mineral aggregates were concentrated to a density greater than three by 
gravitational separation with diiodomethane to facilitate the preparation of 
fission gas samples suitable for mass spectrometric analysis. Chemical analyses 
were performed on samples of the mineral concentrates to determine uranium 
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concentrations. Lead isotope ages of the minerals were obtained by Farquhar 
(1959) and by Collins et al. (1954) of the geophysics group at the University of 
Toronto. 

The stable tracer isotopes Xe!*, Kr®°, and Kr*? were prepared by neutron 
irradiation of iodine and bromine salts in the N.R.X. reactor and standardized 
by diluting against a standard mixture of normal xenon and krypton. 

The apparatus used to extract and purify about 10~‘ ml? at S.T.P. of xenon 
and krypton from a uranium mineral sample is shown in Fig. 1. The vacuum 
induction furnace, powered by a 6-kw Ajax Northrup mercury arc converter, 


Induction 
Furnace 


To Pumps 





Sample 
Tube 


Ti C1 Te Ce Calcium 
Furnace 


Fic. 1. Inert gas separation system. 


could be charged with about 100 g of the crushed mineral concentrate. A 25- 
liter reservoir was connected at the points shown to receive gases evolved 
from the heated mineral and to mix them with a sample of tracer isotopes. 
The open traps T, and T, together with the activated charcoal traps C; and 
C. were used to partially remove the extraneous gases by condensation and 
preferential adsorption on cooled charcoal. The final purification was by 
calcium vapor in the calcium induction furnace, which could be isolated by 
means of mercury ventils. The extraction was considered complete if reheating 
the mineral to 1400° C for 8 hours evolved less than 0.5% of the fission products 
extracted in a first heat to 1300° for 8 hours. 

The inert gas mixtures were analyzed in a 10-in. directional focusing mass 
spectrometer using an electron multiplier for a detector. The multiplier out- 
put was integrated and amplified by a vibrating reed electrometer and dis- 
played on a recorder. A resolving power of about 625 allowed doublet separation 
of the hydrocarbon background from the xenon and krypton isotopes. The 
relative abundances of normal xenon and krypton reported by Nier (1950) 
were used to correct for instrumental mass discrimination. 


RESULTS 


The primary mass spectrometric data provide relative and absolute xenon 
and krypton isotopic abundances. In our analyses of these data, relative 
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xenon isotopic abundances only are required for the normal xenon correction 
and for the analysis of the fission products into their components. Absolute 
abundances, on the other hand, are required to determine absolute yields and 
to measure any loss of fission product gas from the mineral. The normal xenon 
and krypton corrections were made on the basis of Xe", Kr®°, and Kr®, 
since these isotopes are not formed in natural fission in appreciable amounts 
(Kennett 1956). 

The fission product abundances were normalized to 100 at Xe" for com- 
parison and are shown in Table 1(1). Table I(2) gives the total number of 
fission product Xe** atoms in each sample as determined by the isotope dilu- 
tion method. Since the Xe!®® abundance relative to Xe!* is considered to be a 
good indication of the percentage of neutron fission in the mineral (see above) 
the minerals were arbitrarily arranged in order of increasing Xe!*® abundance. 
The Xe’ is preceded in its mass chain by the long-lived I'*® and corrections 
were applied to the Xe!?® abundances for the retention of a small part of the 
129 mass chain as I'*® at the time of extraction. Table I(6, 7) gives the ratios 


TABLE I 
Absolute abundances of xenon and krypton fission products 








Mineral 
Beaver Great Bear 
Fission Cinch Lake Jahala Lake Lodge Lake Lake Eagle Mine Belgian Congo 
product pitchblende uraninite pitchblende pitchblende pitchblende pitchblende 
(1) 
Kr 0.854 + .030 0.728 + .022 0.865 + .044 1.71+.02 1.75+.05 2.56+.03 
kre 2.52 +.15 2.45 +.09 1.98 +.21 3.94+.09 4.23+.18 §.40+.13 
Kr 15.79 +.21 12.54 +.15 14.37 +.18 16.49+.17 15.43+.20 17.95+.20 
Xel29 0.460+.021 0.488 + .009 0.800 + .073 2.06+.01 2.29+.02 3.38+.02 
Xel3t 9.88 +.06 9.77 +.07 11.41 +.11 15.43+.12 15.85+.11 20.51+4.15 
Xelsz 57.97 +.27 58.18 +.48 58.44 +.33 59.76+.43 59.05+.59 61.39+.37 
Xel% 83.36 +.36 83.72 +.63 85.27 +.38 89.17+.75 89.834 .84 94.46+.50 
Xe =§=©100.00 +.42 100.00 +.94 100.00 +.41 100.00 +.95 100.00 + .83 100.00 + .65 
(2) 
Xe!* atoms 3.312 7.710 1.352 6.033 7.639 4.838 
(X10) 
(3) 
Lead age 1120+50 1740 +50 785 +60 1400 +50 1700 +30 635 +20 
(X 108 y) 
(4) 
Um, % 29.0+.3 30.4+.4 41.2+41.4 43.34.4 45.34.4 56.4+.5 
(5) 
Sample wt., 65.1 128.0 103.8 69.7 99.5 115.3 
g 
Normal Kr (8) 
ae 130 .640 -440 -096 .180 -082 
Fission Kr 
Normal Xe (7) 
oe -0050 -0240 .0061 .0055 -0066 


Fission Xe 
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of normal to fission product gas for krypton and xenon respectively. These 
ratios show the extent of normal gas contamination in each sample. The errors 
are the standard deviations based on seven mass spectrometer double scans. 
Duplicate analyses showed that the absolute abundances were reproducible 
within this experimental error. Lead ages and uranium contents of the six 
mineral concentrates are included in Table I(3) and Table I(4, 5) respectively. 


Absolute Yields; Spontaneous Fission 

The absolute yield calculation depends on the independent determination 
of the number of fissions which have taken place and the amount of a given 
fission product which has been formed. 

The total number of U** atoms which have undergone spontaneous fission 
in a mineral sample since deposition was calculated from the U** content, 
the age, and the rates at which the U?* nucleus decays by alpha emission and 
spontaneous fission according to the relation 


N = W(At/da) (e**" — 1) 


where W = the present number of U* atoms in the mineral sample, 


At = the U8 spontaneous fission decay constant, 
= 8.62+.30X10-"/y (Segré 1952), 
Aa = the U**® alpha-decay constant = 1.54X10-'°/y (Kienberger 
1949), 
T = the period during which fission products accumulated (assumed 


to be given by the lead age of the mineral). 


The fraction of the fission gas (xenon and krypton isotopes) in each mineral 
sample which resulted from the spontaneous fission of U** was obtained from 
a breakdown of the results for each natural fission gas sample into its various 
fission components. 

Fleming and Thode (1953) and Wetherill (1953) obtained fairly consistent 
results by assuming only two fission processes in nature, the spontaneous 
fission of U* and the neutron fission of U*, Since the mass 129 chain yield 
is about 1% for neutron fission of U** and very nearly zero for spontaneous 
fission of U** (the curve for spontaneous fission is much steeper at this point) 
they used the yield of this isotope as an index of the proportions of neutron 
fission of U** and the spontaneous fission of U?* which have occurred in a 
given uranium mineral. A preliminary analysis of the data in Table I showed 
that it was not possible to separate the six mineral patterns into a U*® neutron- 
induced component and a U**® spontaneous fission component satisfactorily. 
However, a satisfactory analysis was obtained when a third component, 
U8 neutron-induced fission, was considered along with the two processes 
proposed by Fleming and Thode. This analysis made use of the known relative 
abundances of the xenon and krypton isotopes for the two neutron fission 
processes separately and a comparison of these with the relative abundances 
obtained for the fission gases extracted from the various minerals (Table I). 

In this connection, the U** spontaneous fission pattern was determined by 
a successive approximation procedure, in which the Cinch Lake pattern, 
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corrected for its small neutron fission component on the basis of the Xe!” 
yield as described above, was used as the first approximation. A three-com- 
ponent breakdown, using the Xe'*!, Xe!*4, and Xe! relative abundances in 
the three components and the five remaining mixtures, produced a more exact 
Xe!” relative yield from the neutron fission components. 

The criterion for a satisfactory three-component analysis is that each 
mineral sample gives the same relative abundances of the xenon isotopes 
after subtraction of the two neutron-induced fission components. In other 
words, the spontaneous fission pattern (U?%*) obtained should be the same 
for each sample. It is seen from Table II that the xenon pattern agreement 


TABLE II 


Fission product patterns corrected for neutron fission 





Mineral 
Beaver Great Bear 
Fission Cinch Lake Jahala Lake Lodge Lake Lake Eagle Mine Belgian Congo 
product pitchblende uraninite pitchblende pitchblende pitchblende pitchblende 
Krés 0.519+.031 0.412+ .023 0.175+.048 0.25+.02 0.15+.06 0.00+.06 
Kre 1.94 +.16 1.92 +.09 0.74 +.23 1.36+.11 1.47+.22 0.85+.20 
Kr 15.10 +.22 11.68 +.16 12.82 +.20 13.414.21 11.73 4.25 12.46+.29 


Xel29 0.000+.019 +0.060+.015 -—0.113+.086 + 0.06+.06 +0.15+.06 +0.01+.13 
Xelst 8.31 +.10 8.31 +.05 8.31 +.11 8.314 .07 8.314.12 8.31+.16 
Xels2 57.54 +.17 57.79 +.21 57.58 +.20 57.48+.15 56.544.24 57.394.31 
Xelu 81.54 +.27 82.04 +.24 81.66 +.14 81.54+.15 81.54+.27 81.544 .23 
Xess =. 100.00 +.24 100.00 +.24 100.00 +.20 100.00 + .17 100.00 + .27 100 .00 + .32 


among samples is excellent. The disagreement among the krypton patterns 
is discussed below. 

The percentages of the three components, spontaneous fission of U”%, 
neutron fission of U?*, and neutron fission of U8, obtained from the analysis 
are given in Table III(1). Included also in Table III(2) are the ratios of the 


TABLE III 
Absolute yield of Xe!** and loss of spontaneous fission products 














Mineral 
Beaver Great Bear 
Cinch Lake Jahala Lake Lodge Lake Lake Eagle Mine Belgian Congo 
Component pitchblende uraninite pitchblende  pitchblende  pitchblende _ pitchblende 
(1) 
U2 spontaneous fission, 95.86 96.15 91.82 81.65 80.42 68.65 
% 
U5 neutron fission, % 4.14 3.85 8.18 15.19 17.22 25.23 
U8 neutron fission, % 0.00 0.00 0.00 3.16 2.36 6.12 
(2) 
Xess 
es 6.30+.38 4.38+.18 1.59+.26 4.78+.23 3.214.13 3.514.17 
Fissions 


(U8 spontaneous fission) 


(3) 
Xe! loss, % 0.0+6.0 30.5+6.6 74.8247.3 24.1+47.1 49.0+6.3 44.346.5 
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number of atoms of Xe" produced to the number of fissions (spontaneous 
fission of U***) for each sample. These ratios, expressed as percentages, can 
be interpreted as absolute yields only if fission products have not been lost 
from the mineral. This ratio for the Cinch Lake mineral agrees with other 
absolute yield estimates (Fleming and Thode 1953; Wetherill 1953) and there- 
fore suggests that inert gas losses from this mineral have been very small. 
If the lead ages are assumed to be the periods during which spontaneous 
fission products accumulated, then the low ratios of Xe! to the number of 
fissions for each of the other five minerals must be due to a loss of the rare 
gases during geological time, and these losses are given in Table II1(3). If 
we assume little or no loss of gas from the Cinch Lake sample during geological 
time, the data shown in Tables II and III should give the absolute yields of 
the xenon and krypton isotopes from the spontaneous fission of U**. These 
yields are given in Table IV(1). 


TABLE IV 


Absolute yields in U8 spontaneous fission 


Absolute yields 








Absolute yields (*Parker and Kuroda, 
(This work) tAshizawa and Kuroda) 
Mass Element Yield, % Element Yield, % 
(1) (2) 

83 Kr 0.0327 + .0028 

84 Kr 0.122 +.012 

86 Kr 0.951 + .057 

99 Mo* 6.44.5 
131 Xe 0.524 +.031 It 0.44.1 
132 Xe 3.63 +.22 It 3.64.4 
133 It 1.54.3 
134 Xe 5.14 +.31 It 5.24.5 
135 It 5.14.5 
136 Xe 6.30 +.38 

DISCUSSION 


Tables 1(4) and III(1) indicate that there is a correlation between the per- 
centage of uranium in a mineral concentrate and the percentage of neutron 
fission. This relationship is shown in Fig. 2. The percentage of uranium in 
the mineral concentrate should be approximately equal to the percentage of 
uranium in the mineral itself, since most of the rock matrix was removed by 
the gravitational separation with diiodomethane. 

It had previously been assumed that the presence of the rare earth elements, 
which have large thermal neutron absorption cross sections, was largely 
responsible for the low proportion of neutron fission in the uraninites. The 
occurrence of U** neutron fission in uranium minerals suggests much larger 
effective neutron energies than had been suspected, and indicates that the 
rare earths may have had little effect on the amount of neutron fission in a 
uranium mineral. The fact that the Jahala Lake uraninite and Cinch Lake 
pitchblende, with about equal uranium concentrations in the concentrates, 
have the same proportion of neutron fission supports this conclusion. 
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Fic. 2. Relation between the amount of neutron fission and uranium concentration in a 
uranium mineral or ore concentrate. 


Tables I(4) and III(1) also suggest that the proportion of U?* neutron fission 
in the neutron fission component decreases with decreasing uranium con- 
centration in the mineral concentrate. The presence of two neutron fission 
processes, each with a very different energy dependence, allows some con- 
clusions regarding the effective energy of the neutrons in a uranium mineral, 
U8 has a neutron fission cross section which is less than 1 barn in the 1-Mev 
energy region and has a threshold at about 0.5 Mev. U?*, on the other hand, 
has a neutron fission cross section which varies from about 1000 barns in the 
thermal energy range to about 1 barn in the 1-Mev region. The ratio of U8 
to U* neutron fission in normal uranium is therefore very sensitive to the 
effective neutron energy and decreases by a factor of about 10° when this 
energy is reduced from 2 Mev to 0.5 Mev. Because of this sensitivity, and the 
effects of mineral age on the U?**/U?8 ratio, the correlation between uranium 
concentration in the mineral and the proportion of U** neutron fission in the 
neutron fission component is not expected to be quantitative. 

It may be concluded from these correlations that the materials which 
dilute uranium in minerals strongly absorb neutrons, but reduce only slightly 
their mean energy. 

Internal consistency in the spontaneous fission pattern obtained from the 
analysis of the natural fission patterns into their three components required 
the use of a Xe!” yield in U** neutron fission of 0.71%, which is about 20% 
lower than the I'*° yield reported by Purkayastha and Martin (1956). These 
authors first drew attention to this discrepancy when they extracted the 
Xe!2® and I'®° from two uranium minerals and found that the I'” abundance 
was about 20% lower than would be expected from the Xe!” content. No 
satisfactory explanation for this discrepancy has been found. 
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The ratio of krypton to xenon fission products from spontaneous fission 
(Table II) is about 20% higher in the Cinch Lake mineral than in the other 
five uranium minerals. This discrepancy is attributed to the fractionation of 
krypton relative to xenon as the rare gases diffuse from a mineral. The ratio 
of krypton to xenon obtained from the Cinch Lake sample should therefore 
be a minimum value for U*** spontaneous fission. Since the loss of fission gas 
from this sample must be small we have taken it to be the true ratio. The 
large losses of these gases by most uranium minerals prevent the general use 
of accumulated rare gas fission products for age determination. Nevertheless, 
the fractionation which accompanies loss should be a sensitive test of the 
validity of a ‘‘fission product age’”’, i.e. the period during which fission products 
have accumulated without loss. 

Recently, Parker and Kuroda (1958) and Ashizawa and Kuroda (1957) 
have reported the absolute yields of Mo®* and several iodine isotopes from 
the spontaneous fission of U**. These yields, along with the absolute yields of 
xenon and krypton from the present analysis, are shown in Table IV. The 
good agreemer.t between the xenon and iodine yields at masses 131, 132, and 
134 is in accord with the assumption that no fission product xenon and krypton 
had been lost from the Cinch Lake mineral. Segré’s (1952) value of the U* 
spontaneous fission half life has been used for all of the absolute yields obtained 
thus far, and therefore, the reliability of this half life does not enter into the 


comparison. 
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A *II“II ELECTRONIC BAND SYSTEM OF THE 
FREE NCO RADICAL! 


R. N. Dixon? 


ABSTRACT 


A series of red-degraded absorption bands has been observed between 2650 A 
and 3200 A and is attributed to a “MII transition of the NCO'radical. The bands 
probably represent a progression of the upper-state stretching vibration vj. The 
rotational structure of one band has been analyzed. Diffuseness in some of the 
bands indicates predissociation of the upper state, and is discussed in terms of 
the dissociation energy of NCO. 


INTRODUCTION 


A system of violet-degraded emission bands near 4400 A was observed at 
low dispersion by Holland and Style in the fluorescence of a number of aliphatic 
cyanates and isocyanates. These same bands were also observed in absorption 
during the flash photolysis of similar compounds (Holland, Style, Dixon, and 
Ramsay 1958), and showed rotational structure under high resolution. The 
vibrational and rotational analyses of these bands have been carried out 
(Dixon 1959, hereafter referred to as Paper I), and show that the bands are 
due to the A(?Z+)-X(*II,) transition of the NCO radical. During this flash 
photolysis experiment a series of red-degraded bands was observed in the 
region 2650 A — 3200 A. These bands are more complex and less resolved than 
the 4400-A bands, but the structure of strong bands in the series suggests that 
they are ?II-*II bands. They are attributed to the B(?II,)—X(?Il,) transition 
of NCO. 

EXPERIMENTAL PROCEDURE 


The ?II-*II bands of NCO were observed during the flash photolysis of 
HNCO vapor, using the same apparatus and experimental conditions as for 
the *=+—*II bands (Paper I). The bands at wavelengths greater than 2900 A 
were photographed on Eastman Kodak I-O plates in the third order of a 
21-ft concave grating spectrograph, with 12 traversals through the 50-cm 
reaction tube. The complete region from 2600 A to 3200 A was also photo- 
graphed in the fourth order of the grating, but lack of light-gathering power 
necessitated a reduction in the absorption path to two traversals through the 
reaction tube. 

Two plates of the strong sharp features in the spectrum were each measured 
twice on a comparator equipped with a photoelectric scanning device similar 
to that described by Tomkins and Fred (1951). The diffuse features were 
measured from a print of the spectrum. A hollow-cathode iron-arc lamp was 
used for wavelength calibration. 


‘Manuscript received September 24, 1959. 

Contribution from the Division of Pure Physics, National Research Council, Ottawa, Canada. 

Issued as N.R.C. No. 5461. 

2National Research Council Postdoctorate Fellow, 1957-59. Present address: Department of 
Physical Chemistry, The University, Sheffield 10, England. 
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VIBRATIONAL STRUCTURE 


The wavelengths and vacuum wave numbers of the band heads and absorp- 
tion maxima are given in Table I, and portions of the spectrum are illustrated 


TABLE I 


Band heads and absorption maxima for the B(*II1) — X(#II) system of 
‘CO, with possible assignments 











Feature Mair), A v(vac.), cm! Assignment 
VW 3167.14 31,565. 1 

VW 3156.48 31,671.7 

S head 3149 .873 31,738.12 Rz 000-000 
S head 3148 .388 31,753.09 R, 000-000 
VW 3142.34 31,814.3 

VW 3141.99 31,817 .7 

W head 3059 . 642 32,674.06 

S 3051.30 32,763 .4 Ro v; 

S head 3045. 877 32,821.72 Ri v3 
VW head 3038. 351 32,903 .02 

W b 2973.59 33,619.6 

S head b 2964.31 33,724.8 Ry 21 
Md 2959 . 50 33,779 .7 

S head b 2953 .92. 33,843 .4 R, 2} 
Md 2949 .2 33,897 

Md 2936.6 34,043 Re v3 
Md 2935.7 34,054 Ri v3 

WwW 2873.06 34,795 .9 

M b 2866 . 93 34,870.3 R 3p} 
VW b 2865 .50 34,887 .7 

VW b 2864 .02 34,905 .7 

Wd 2849.7 35,081 Rez vitv; 
Wd 2848.6 35,095 Ry vitv3 
VW \ 2843.51 35,157 .5 

VW 2843 .30 35,160.1 

W 2786.78 35,873 .2 

MW head 2783 .27 35,918.4 R 4p; 
W 2780.47 35,954 .6 

VW 2712.72 36,852 .5 

VW 2702.19 36,996 .0 

W 2690 .38 37,158 .5 

W 2677 .24 37,340.8 


Note: VW very weak, W weak, M medium, S strong, b broad, d diffuse. 


in Fig. 1. It is seen that the bands are grouped in regions about 1000 cm~! 
apart. These groups most probably arise from the excitation of successive 
quanta of the upper-state stretching frequency vj. The bands at 3148 A, 
3046 A, and 2954 A are approximately equal in intensity, and there is no trace 
of any bands to the red of the 3150-A group. The 3148-A band must therefore 
be the 000-000 band of the electronic transition. 

Upon closer inspection of the features in the groups of bands it is seen that 
in each of the 2950-A, 3050-A, and 3150-A groups there are two prominent 
heads which have the appearance of the R; and R: heads of *II-*II bands. 
However, the doublet splitting in these three bands is far from constant, and 
does not even vary in a linear manner. Furthermore, it is not possible to select 
a series of band heads in a vibrational progression with a constant second 
difference. This suggests that there are vibrational perturbations in the upper 
state. It is not difficult to understand how perturbations leading to such 
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irregularities could arise. In COz and similar molecules, and in the A (22+) 
state of NCO (Paper I), v; ~~ 2v2, and there is a Fermi resonance between those 
levels v;,¥2,03 and v;—1,v2+2,v3 which are of the same species. This may also be 
true in the B(?II) state of NCO. In addition, as this state is a II state, the 
electronic angular momentum Ah/2m may couple with the angular momentum 
lh/2m of the bending vibration such that A and / are no longer good quantum 
numbers. The resultant angular momentum is then characterized by 
K = !A+/|. The energy levels of such a II state have been discussed by 
Renner (1934). In Paper I it was shown that the X(?II) state of NCO is 
subject to such an effect, with additional effects due to the coupling of the 
angular momentum Kh/2z with the spin angular momentum Zh/2z. The levels 
of this latter system have been investigated theoretically by Pople (1959). Jt 
may be shown that the spin doubling in the ?II (vibronic) levels of even v2 is 
much smaller than in the zero level. Thus, if there is a Fermi resonance between 
the levels v;,v2,v3 and v;—1,v2+2,v3 in the B(?II) state of NCO, this would 
involve the interaction of levels with different effective spin-coupling constants, 
and would give a very complicated vibrational pattern. In consequence of 
these possibilities the vibrational assignments of Table I can be considered as 
only tentative. 

The spectrum shows a number of diffuse features, presumably as a result of 
predissociation of the upper state. The two pairs of diffuse features near 
2936 A and 2849 A show doublet splittings similar to that in the 000-000 band. 
It is unlikely that these are hot bands as there are no similar features near the 
3046-A and 3149-A bands. They are therefore assigned to v} and vi+v}. 

The following vibration frequencies are derived from the mean values of 
the frequencies of the R; and Re heads of the 000-000, 100-000, and 001-000 
bands; vj = 1047 cm~', v5, = 2303 cm~!. A complete determination of the 
stretching force constants in the B(?II) state of NCO cannot be achieved 
without data from an isotopic molecule. However, if the above frequencies are 
substituted into the equations for the bond stretching of a linear XYZ molecule 
with no interaction force constant ki: (Herzberg 1945, p. 173), then either 
Ren = 17.3X10° dynes/cm and kcg = 7.4X105 dynes/cm, or kcoy = 7.0X105 
dynes/cm and keg = 18.4X10* dynes/cm. If ky2 is included in the equations it 
is found that with these vibration frequencies real solutions for the force 
constants are possible for any ki: >—1.9X10° dynes/cm, at which value 
Ron = 11.2X10° dynes/cm and kceo = 11.9X105 dynes/cm. Thus, either one 
bond of NCO in the B(?Il1) state is much stronger than the other, or there is a 
negative interaction stretching force constant kiz. (which implies that the 
stretching of one bond of NCO in this state weakens the other bond). 


ROTATIONAL STRUCTURE 


Several of the bands in the spectrum show rotational fine structure. Although 
the 000-000 band has the two heads expected for a *II(a)-*II(a) band, a 
rotational analysis of this band was not possible, as the structure appears to 
be considerably perturbed. This could be due to perturbation by higher levels 
of the A(22+) state, which is only 9000 cm! below the B(?II) state. Examina- 
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tion of the other bands in the spectrum showed that only in the 3046-A band 
was the fine structure sufficiently extended for a rotational analysis. 

The 3046-A band has two subbands, each with one P and one R branch. The 
lines in the region of the subband origins are very weak. Consequently it was 
not possible to determine the J numberings from the positions of the zero 
gaps, and the J values of Table II were assigned using lower-level combination 
differences. In columns 6 to 9 of Table II the lower-level combination dif- 


TABLE II 


Vacuum wave numbers and assignments for the lines of the 3046-A band of NCO, and lower-level 
combination differences 


R(J-1) RxAJ—1) 
J P, Ri P, R, —P\J+1) —PAJ+1) A2Fi(J)t  A2F:'(J)t 
i 32,815.30 — 3.13 
2} 32,811.55 816.11 4.42 4.65 4.69 
3} «810.88 816.59 32,752.11 6.31 6.21 6.26 
44 809.80 817.30 752.63 7.68 7.76 7.82 
5} ©=—S 808.91 817.71 753.39 (9.25) 9.31 9.31 9.38 


63 808 .05* 818.36 32,743.32 754.09 (10.68) (10.61) 10.86 10.95 
73 807 .03* 818.74 742.78* -754.77 12.49 (12.45) 12.41 12.51 
83 805 . 87 819.24 741.64* 755.24 13.90 14.12 13.96 14.07 
94 804.84 819.77 740.65 755.93* (15.63) 15.59 15.51 15.64 


103 803 .61* 739.65 756.51 (17.30) (17.19) 17.06 17.20 
113 802.47* 738.74 757 .06 18.74 18.61 18.76 
123 801.23 737.77 757.68 20.29 20.17 20.33 
133 800.01 736.77 758.24 21.94 21.71 21.89 
14} 798.75 735.74 758.76 23.54 23.26 23.45 
153 797.61 734.70 759.31 25.03 24.82 25.02 
164 796.22 733.73 759. 68* 26.62 26.37 26.58 
173 794.88 732.69 760.15 (28.11) 27 .92 28.14 
183 793 .42* 731.57 29.69 29.49 29.73 
193 792.01 730.46 31.02 31.27 
203 790.81 (Head 729.35 32.54 32.81 
214 789.32, 9.13 821.72) 728.06 34.12 34.39 
224 787.85, 8.65 727 .02 35.67 35.95 
233 786.30, 6.11 725 . 87 37.22 37.52 
24} 784.89, 4.54 724.84 38.77 39.08 
253 723.68 40.32 40.64 

41.86 42.20 


264 781.46 
274 780.31, 9.96 43.41 43.76 


*Blended with lines of the 1-0 band of the A(*II) — X(#Z~) system of NH. 

tCalculated from rotational constants of Paper I. 
ferences obtained from this band are compared with values calculated from 
the rotational constants of Paper I. The combination differences from the 
subband of longer wavelength are available to sufficiently high J values to 
prove that this is the subband *II, <— *IIy (the ground state is known to be an 
inverted state). 

The rotational term values for doublet states of Hund’s coupling case (a) 
have been calculated by Hill and Van Vleck (1928) for any magnitude of the 
spin-coupling constant A: 

(1) F(J) = Bi(J+})?—A*]4[B*(J+4)?+44 (A —4B)A*}}. 
With the addition of a centrifugal stretching term this expression was found 
adequate for the representation of the rotational term values of a *II level 
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(00'0) and a 2A level (0170) in the ground state of NCO (X, ?II,). Equation (1) 
shows that the mean energy of the F,(/J) and F2(J) levels is independent of 
the magnitude of A: 


(2) AFi(J)+Fx(J)] = Bi(J+3)?—A’]. 


There are insufficient lines in the R branches of the 3046-A band for the 
determination of an accurate B’ value from upper-level combination dif- 
ferences. Consequently the upper-level term values were computed by adding 
the rotational term values of the appropriate lower levels (calculated from the 
data of Paper I) to the measured frequencies of the lines in the band. The mean 
value for each J of the resultant 7T\(J) and 73(J) was plotted against 
[(J +3)?—1]. The scatter of the plot was too great to allow a determination of 
the upper-state centrifugal stretching constant D’. The term values 7’(J/) 
were therefore corrected to allow for a value of D’ put equal to D” (= 1510-8 
cm~'). Then from the intercept and slope of the line: 


(3) vo = 32,781.13; cm™!, B’ = 0.3765+0.0002 cm-!. 
If .1 > B, the square root of equation (1) can be expanded to show that 
(4) F,(J)—Fi(J) = A|A —2B|+(2B*)/(A|A —2B})[(J +4)?—A?] 


and thus the rotational levels of the two sublevels of a *II level (A = 1) can 
be represented using effective B values of B+ B?/|.1 —2B{. Values of the left- 
hand side of equation (+) were obtained from the separations of lines in the 
two subbands and the lower-level combination differences, according to the 
equations 


(5) FJ) — FJ) = [Fy (J+) — Fy (J +1) -[PaVJ) — P2(J)] 
(6) = [Fy (J-1) -— F'\(J-1)]-[ Ri) — RoJ)], 


and were plotted against [((J+3)?—1]. This gave a straight line of intercept 
31.6 cm~! and slope +0.012) cm~!. Thus A’ = —30.8 cm—!. The value of the 
slope of the plot calculated according to equation (4) with these values of 
A’ and B’ is 2B’’/\|4’—2B’| = +0.0090 cm-'. Therefore, although (Bo en— 
Bi.) is of the correct sign and order of magnitude, the rotational term 
values of this upper level cannot be accurately represented by equation (1). 
Possibly this lack of agreement between observation and calculation is due to 
a Fermi resonance, as discussed above in the section on the vibrational 
structure of the bands, but as was pointed out above, even the 000-000 band 
appears to be perturbed. 

The value of B’ = 0.3765 cm~! determined above cannot be a ‘true’ B value 
for the level, as the rotational structure does not obey equation (1). In addition, 
this B’ value is an effective value for an excited vibrational level of the B(?11) 
state of NCO. Even so it can be used to give an approximate value of a mole- 
cular dimension for NCO in this state. The individual bond lengths ryc and 
reo cannot be determined without data from an isotopic molecule, but it is 
possible to give an upper limit to the distance between the nitrogen and 
oxygen atoms consistent with this B value. To calculate this limit the carbon 
atom is placed at the center of gravity of the molecule, and the distance 
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corresponding to B’ = 0.3765+0.0030 cm? is (tyro +rco) <2.45+0.01 A. The 
molecular dimensions determined for the A(?2+) and X(?II) states of NCO 
from the values of Bo given in Paper I are: 


(7) A(?3+), (trc+roo) < 2.369 A, 
(8) X(?I1), (mwc+reo) < 2.408 A. 


DISSOCIATION AND PREDISSOCIATION 


The only published kinetic study involving the NCO radical is on the 
thermal decomposition of potassium cyanate (Shushunov and Serdyuk 1953). 
The initial states in the pyrolysis are thought to be 

KNCO — K + NCO, (9) 
NCO —N + CO, (10) 


of which the first appears to be the rate-determining step. As potassium cyanate 
decomposes at about 700° C the activation energy of its decomposition is prob- 
ably of the order of 70 kcal/mole. Thus the dissociation energy of NCO into 
N+CO cannot be greater than 70 kcal/mole (3.0 ev). An independent estimate 
of the dissociation energy of N—CO may be made from the heat of formation of 
HNCO (36.5 kcal/mole in aqueous solution, Bichowsky and Rossini 1936). If 
the heat of solution of gaseous HNCO is assumed to be 10 kcal/mole, and the 
NH bond strength is 100 kcal/mole, then the dissociation energy of N—CO is 
64 kcal/mole (2.8 ev), with an error probably no greater than 20 kcal/mole 
(0.8 ev). 

The possible electronic states of linear NCO that can be formed from the 
low-lying electronic states of N+CO, or O+CN, were discussed in Paper I. 
These are indicated in Fig. 2, which is a schematic correlation of the three 


k ca!/mole 


si N+co NCO NG+0 J !20 





Fic. 2. Correlation diagram for NCO and its dissociation products N + CO or NC + O. 
The possible electronic states arising from each combination of dissociation products are 
indicated in braces. 
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known states of NCO with the dissociation products, taking account of the 
non-crossing rule. The dissociation energy of N—CO has been arbitrarily made 
69 kcal/mole (3.0 ev) in this figure. 

The observed diffuseness of some of the bands in the spectrum above 33,709 
cm~! (4.2 ev) is presumed to be due to predissociation of the B(?II) state of 
NCO, although the possibility that the diffuse features are unresolved groups 
of lines cannot be ignored, especially as the spectrum appears to be perturbed. 
If this predissociation is produced by the ‘2- state, leading to N(*S,)+ 
CO('=+), it follows that the dissociation energy of N—CO < 96 kcal/mole 
(4.2 ev). However, it would be surprising if a “‘forbidden’’ 42-,?II predissocia- 
tion results in a noticeable diffuseness in the spectrum. On the other hand, if the 
predissociation is produced by the ?2~ state, leading to N(?D,) +CO(‘2+), then 
the dissociation energy of N—CO is less than or equal to 41 kcal/mole (1.8 ev). 
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THE SPECIFIC HEAT OF COPPER FROM 20° TO 300° K' 


Douctas L. MARTIN 


ABSTRACT 


The specific heats of commercially pure cold-rolled copper and of annealed 
and heavily cold-worked 99.999% pure copper have been measured in the tem- 
perature range 20° to 300° K. When results are averaged over the whole 
temperature range of measurement the specific heat of the pure cold-worked 
copper is about 0.15% above that of the pure annealed sample while results for 
the commercially pure cold-rolled material lie in an intermediate position. 
Results on a given sample are reproducible within 0.05%. The entropy of pure 
annealed copper at 298.15° K is 7.92+0.04 cal/°K g-atom. 


INTRODUCTION 


Measurements of the specific heat of copper in the temperature region 20° 
to 300° K are too numerous to detail. Early measurements (see, for example, 
Regnault 1840) suggested that there was a difference in the specific heats of 
cold-worked and annealed copper but careful later comparisons (Eucken and 
Werth 1930; Maier and Anderson 1934) showed that the difference, if any, 
was not greater than the experimental error (~ 0.2%). Later work by 
Giauque and Meads (1941) produced results systematically higher than other 
recent measurements (loc. cit. and Dockerty 1933, 1937) which they attributed 
to the perfection of their sample, arguing that since annealed copper is softer 
than cold-worked material the Debye temperature was lower and the specific 
heat higher. However, they failed to distinguish between elastic and plastic 
properties and in fact the elastic constants of cold-worked material are lower 
than those for annealed copper. Hence, on this argument one would expect the 
specific heat of annealed copper to be lower than that of cold-worked material. 

The primary purpose of the present work was to assess the performance of 
a calorimeter but at the same time it has been possible to resolve the above 


discrepancy. 


EXPERIMENTAL 


The measurements were made with a continuous heating semiautomatic 
adiabatic calorimeter basically similar to that described by Dauphinee et al. 
(1954). This earlier apparatus was constructed so that it could be placed 
between the pole pieces of an A. D. Little electromagnet, and the cramped 
design resulting from this requirement greatly hampered the attainment of 
adiabatic conditions so that above about 150° K it was always necessary to 
surround the cryostat with a liquid bath maintained within a few degrees of 
the calorimeter temperature. Figure 1 is a diagram of the redesigned cryostat 
and it will be seen that the three adiabatic shields are now separated from one 
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Contribution from the Division of Pure Physics, National Research Council, Ottawa, 


Canada. 
Issued as N.R.C. No. 5463. 


Can. J. Phys. Vol. 38 (1960) 


17 











18 CANADIAN JOURNAL OF PHYSICS. VOL. 38, 1960 


another and the rest of the cryostat. The calorimeter and shields are assembled 
in a ‘cage’ which is then screwed onto the bottom of the can and electrical 
connections established by means of copper ‘jumpers’ soldered between the 
two rings of teflon-insulated terminals. 

The electronic control of the shields and the measuring circuit are exactly 
as described previously except that time measurement is now by a stop clock 
driven from a standard frequency supply and triggered by a cam and micro- 
switch on the balancing motor shaft of the Leeds and Northrup potentiometer 
recorder measuring the thermometer resistance. 

A slightly different calorimeter vessel from that used previously (Dauphinee 
et al. 1954, 1955) has been used in this work and can be seen in Fig. 1. The 
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Fic. 1. The cryostat, shield system, and calorimeter. 
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sealing of the calorimeter is carried out in a helium atmosphere and the 
Wood's metal is manipulated with a non-tinning molybdenum spatula, so that 
the weight of metal remains constant. The only conclusive test for the vacuum 
tightness of the seal is measurement of the calorimeter equilibrium time after 
assembly. (The tubular part of the calorimeter was made from copper sheet 
0.004 in. thick with a hard-soldered lap joint. For such a thin wall this method 
is much superior to the customary turning.) 

The single thermocouple control of the adiabatic shield presupposes a 
uniform surface temperature of the calorimeter vessel. The large area of the 
surface covered by the heater contributes towards satisfying this requirement 
but it is also necessary to have a specimen which is a good conductor fitting 
closely the inside dimensions of the calorimeter vessel. 

If a suitable specimen is used, the results are independent of heating rate 
over a range of at least 2 °K/hr to 30 °K/hr. The surrounding bath temperature 
now has no influence on the results. 

The three probable causes of error in these measurements are in the deter- 
mination of temperature, in the power measurement, and in the non-attainment 
of adiabatic conditions. The copper temperature scale of Dauphinee and 
Preston-Thomas (1954) was used. Experience with a larger number of copper 
thermometers indicates that their original estimate of errors was too optimistic 
and should be doubled. Resistors and standard cells have been checked against 
standards. Determination of the calorimeter equilibrium times at various 
temperatures indicates that errors due to thermal lags in the calorimeter are 
negligible. Drift checks with the adiabatic shield temperature offset from that 
of the calorimeter show that errors due to non-attainment of adiabatic condi- 
tions during normal operation are unlikely. The limiting factor in shield 
control is stray thermal voltages in the thermocouple circuits. Hence results 
above 80° K are expected to be accurate to +0.2%, the accuracy decreasing 
to +2.0% at 20° K. The reproducibility of results on a given specimen is of 
the order of 0.05%, as will be detailed in the next section. 

The copper samples were prepared from two different sources of material. 
Two of the samples were prepared from 99.999% pure copper obtained from 
the American Smelting and Refining Co., New Jersey. This was melted into 
graphite crucibles by induction heating under high vacuum. The first sample 
(the ‘annealed sample’) was cooled slowly over a period of about four hours 
and measured in the ‘as cast’ condition, no machining or other operations 
having been carried out. After measurements were completed the sample was 
etched and found to consist of a large crystal with much smaller crystals at its 
edges. The second sample was heavily cold worked below room temperature, 
using a hydraulic press, until its length was increased by 50%. It was then 
turned to shape, care being taken to prevent heating. X-ray measurement 
showed that the high-angle lines were very weak or absent and low-angle 
lines showed considerable broadening and some displacement. It is inferred 
from these results that the lattice is severely distorted as a result of the cold 
working. By approximate analysis of the line profiles a plot of fraction of 
sample with a given strain against strain was obtained. From this the average 
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strain in the sample was deduced to be about 0.2%. The weights of these 
samples were 65.092 g (annealed sample) and 71.424 g (cold-worked sample). 

The other samples were turned from commercially pure cold-rolled copper. 
Spectrographic analysis showed the chief impurities to be silicon and magnesium, 
both estimated to be less than 0.01%. The first of these samples was solid and 
the second was hollow (made in two parts) so as to change the mass without 
altering the external dimensions. The weights of these samples were 71.138 
and 38.625 g respectively. After completion of measurements the solid sample 
was etched and the average grain size found to be of the order 0.05 mm. 
X-ray measurement on the same stock showed line broadening which, by the 
analysis indicated above, was found to correspond to an average strain of 
about 0.1%, indicating considerable cold work. 


RESULTS 


In this continuous heating apparatus the time elapsed during a compara- 
tively small temperature increment is measured (normally less than 0.5° K in 
the present work). In this way a great number of ‘points’ are recorded and, 
except in regions of special interest, only a sample of these are evaluated and 
plotted. The smoothed values in Table I are read from large-scale graphs. 

Results on the various samples examined during this work are compared by 
summing the smoothed specific heat results for equal temperature increments 
for one sample and comparing with a similar sum for another. In this way 
mean systematic deviations between samples may be found. The method is 
equivalent to a comparison of enthalpies. 

As a check on the reproducibility of the apparatus the heavier commercially 
pure sample was measured on two occasions and the results independently 
evaluated and tabulated (columns I(a) and 1(6) of Table I). The deviation 
averaged over the whole temperature range is 0.05%. The lighter commercially 
pure sample (column II of Table I) is systematically high with a mean devia- 
tion of about 0.25%. It is not known whether this is due to the mass difference 
or to the extra machining operations or both. 

The results obtained for the very pure copper are almost the same as those 
for commercially pure material. Again averaging over the whole temperature 
range, the cold-worked sample has a systematic average deviation of about 
0.15% above the annealed sample and about 0.07% above the commercially 
pure sample of approximately the same mass. 


DISCUSSION 


The results for the annealed sample are compared with other work in 
Table II. Giauque and Meads (1941) remark that Dockerty’s work on com- 
mercial copper shows an average deviation of 0.018 cal/°K g-atom below their 
results over the temperature range 50° to 300° K. Satterthwaite et al. (1954) 
have since found for commercial copper an average deviation of 0.017 cal/°K 
g-atom below Giauque and Meads’ results, thus closely confirming Dockerty’s 
work. The present results for commercial copper deviate (taking the average 
of I(a) and I(d)) about 0.020 cal/°K g-atom below Giauque and Meads and 
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TABLE I 
Specific heat of copper (units cal/°K g-atom) 


99 .999% pure copper 
Commercially pure ee 





cold-rolled copper Annealed 

_ Cold- —_—-——- 
7. °K I(a) 1(d) II worked Cp Cy Op 
20 _ —_— _— _— 0.113 0.113 324 
21 0.132 0.132 0.131 0.132 _— _ 
22 _ _— 0.153 — — _- _- 
25 0.230 0.231 0.231 0.229 0.230 0.230 317 
30 0.406 0.407 0.408 0.405 0.404 0.404 313 
35 0.637 0.636 0.637 0.636 0.634 0.634 311 
40 0.896 0.897 0.901 0.896 0.893 0.892 311 
45 1.178 1.182 1.185 1.180 1.179 1.178 311 
50 1.470 1.471 1.477 1.472 1.467 1.465 313 
55 1.769 1.769 1.774 1.770 1.764 1.761 313 
60 2.056 2.061 2.064 2.061 2.055 2.051 314 
65 2.330 2.342 2.341 2.336 2.334 2.328 314 
70 2.592 2.598 2.601 2.597 2.589 2.581 315 
75 2.841 2.834 2.848 2.849 2.841 2.831 315 
80 3.074 3.071 3.089 3.081 3.074 3.061 316 
85 3.289 3.287 3.301 3.291 3.287 3.272 316 
90 3.481 3.488 3.496 3.498 3.488 3.470 316 
95 3.657 3.662 3.678 3.668 3.665 3.645 315 
100 3.811 3.821 3.838 3.824 3.826 3.802 317 
110 4.108 4.111 4.118 4.109 4.116 4.086 317 
120 4.350 4.358 4.367 4.358 4.354 4.317 318 
130 4.563 4.572 4.579 4.572 4.565 4.520 318 
140 4.748 4.749 4.763 4.744 4.749 .4.697 317 
150 4.899 4.902 4.909 4.896 4.903 4.844 316 
160 5.030 5.035 5.042 5.031 5.029 4.963 317 
170 5.143 5.147 5.156 5.145 5.138 5.064 318 
180 5.243 5.245 5.256 5.242 5.237 5.155 318 
190 5.325 5.328 5.347 5.331 5.327 5.238 317 
200 5.401 5.404 5.420 5.409 5.399 5.302 318 
210 5.472 5.471 5.487 5.475 5. 466 5.362 318 
220 5.535 5.531 5.542 5.534 5.527 5.415 319 
230 5.592 5.586 5.593 5.587 5.579 5.460 319 
240 5.645 5.632 5.643 5.636 5.625 5.498 321 
250 5.687 5.681 5.690 5.685 5.672 5.537 321 
260 5.724 5.728 5.737 5.731 5.718 5.571 319 
270 5.764 5.772 5.778 5.773 5.761 5.611 318 
273.15 5.776 5.784 5.790 5.784 5.773 5.620 _ 
280 5.798 5.811 5.815 5.807 5.794 5.636 319 
290 5.828 — 5.844 5.835 5.824 5.660 320 
298.15 5.850 —_ — 5.858 5.843 5.673 — 
300 5.855 — _ 5.863 5.846 5.675 324 





Note: Atomic weight = 63.57; J = 4.186. 


are therefore within about 0.05% of Dockerty and of Satterthwaite et al. 
Hence it is concluded that the present apparatus yields specific heat results in 
very close agreement with those obtained elsewhere. 

From our experimental results it has been shown that the specific heat of 
annealed copper is systematically very slightly lower than that of cold-worked 
material. This conclusion was reached by a comparison of the results obtained 
from the measurement of a number of different samples of copper with the 
same calorimeter and therefore does not depend upon the absolute accuracy 
of the measurements. Now Giauque and Meads compared their measurements 
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TABLE II 


Comparative values for specific heat of copper 
(units cal/°K g-atom) 





This work Dockerty Giauque and 
T,°K 99.999% annealed (1933, 1937) Meads (1941) 


20 0.113 = 0.114 
25 0.230 _— 0.233 
30 0.404 0.401 0.412 
35 0.634 0.628 0.637 
40 0.893 0.900 0.901 
45 1.179 1.188 1.191 
50 1.467 1.489 1.493 
55 1.764 1.791 _ 
60 2.055 2.076 2.074 
65 2.334 2.345 _— 
70 2.589 2.607 2.615 
75 2.841 2.859 — 
80 3.074 3.090 3.091 
85 3.287 3.304 — 
90 3.488 3.498 3.500 
95 3.665 3.669 — 
100 3.826 3.828 3.847 
110 4.116 4.115 4.138 
120 4.354 4.364 4.380 
130 4.565 4.571 4.580 
140 4.749 4.749 4.752 
150 4.903 4.904 4.904 
160 5.029 5.026 5.042 
170 5.138 5.145 5.163 
180 5.237 5.242 5.263 
190 5.327 5.328 5.350 
200 5.399 5.395 5.430 
210 5.466 5.459 5.499 
220 5.527 5.523 5.558 
230 5.579 5.576 5.613 
240 5.625 5.622 5.670 
250 5.672 5.669 5.719 
260 5.718 5.714 5.755 
270 5.761 5.752 5.788 
273.15 5.773 5.763 —_ 
280 5.794 5.789 5.817 
290 5.824 5.824 5.837 
298.15 5.843 5.850 5.848 
300 5.846 5.857 5.850 


on a single sample of annealed copper with results of other workers and con- 
cluded that the specific heat of annealed copper was systematically higher 
than that of cold-worked material. It appears therefore that there must have 
been some systematic error in Giauque and Meads’ results. 

Measurements of the effect of strain on elastic constants have been made 
by Kawai (1930), Smith (1953), and others. For copper the apparent Young’s 
modulus decreases for small elongations but increases again for greater strains. 
It seems clear that the latter effect is due to the development of preferred 
crystallographic orientation as a result of elongation in which case the elastic 
anisotropy results in an apparent increase of Young’s modulus. Hence the 
over-all effect of cold work is a reduction of elastic constants, which Smith 
(loc. cit.) attributes to the elastic oscillations of relatively mobile screw 
components of dislocation loops. 
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The effect of this change of elastic constants on the specific heat is only 
likely to be appreciable at very low temperatures where the dominant vibra- 
‘tional modes are of such long wavelength as to include a number of disloca- 
tions. (There is also a contribution to the specific heat due to the vibration 
of mobile defects (Granato 1958).) Thus in the temperature range covered in 
these measurements no very large difference in the specific heat of annealed 
and cold-worked copper is to be expected. Any residual of the low-temperature 
effects would increase the specific heat of cold-worked material above that of 
an annealed sample. 

The specific heat at constant volume for the annealed copper sample has 
been calculated using the thermal expansion data of Rubin ef al. (1954) and 
the elastic constant data of Overton and Gaffney (1955). From these C, values 
Debye temperatures have been evaluated after subtraction of the electronic 
specific heat (Corak et al. (1955). Figure 2 shows the results and also the limiting 





wy 50 100 150 200 250 300 


TEMPERATURE °K 


Fic. 2. The Debye temperature for copper. Smoothed experimental results are represented 
by the solid line, which has been extrapolated to the 6 value of Corak et al. and Rayne. 


value at the absolute zero taken from calorimetric results below 4.2° K of 
Corak et al. (1955) and of Rayne (1956). (Strictly there is a further correction 
to be applied to the Debye temperature values due to the thermal expansion 
of the copper. This may be calculated using the Griineisen constant for copper 
and is greatest at the high-temperature end where the Debye temperature 
values would be increased by about 6° K.) There are two main features of 
interest in this curve: firstly, the very shallow minimum and, secondly, the 
fact that there is no sign of a decrease of the Debye temperature at the high- 
temperature end, indicating a negligible anharmonic contribution to the 
specific heat even at this temperature. 
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The entropy at 298.15° K has been calculated from the results for the 
99.999% pure annealed sample and is 7.92+0.04 cal/°K g-atom. A value at 
20° K was obtained from the extrapolated @p plot in Fig. 2 and the elec- 
tronic specific heat as determined by Corak et al. (1955). 


CONCLUSION 


Averaged over the temperature range 20° to 300° K, the specific heat of 
heavily cold-worked 99.999% pure copper is about 0.15% above that of the 
annealed material. The reproducibility of measurements on a given sample 
is within 0.05%. 
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THE SPECIFIC HEAT OF A LITHIUM-MAGNESIUM ALLOY. 
THE MARTENSITIC TRANSFORMATION! 


Douctas L. MARTIN 


ABSTRACT 


The specific heat of an alloy of lithium with 0.95 at.% magnesium has been 
measured from 20° to 300° K. The specific heat of the body-centered cubic phase is 
less than that calculated from the Kopp—Neumann rule. A specific heat anomaly, 
due to the martensitic transformation, is observed and has almost the same heat 
content as the corresponding anomaly for pure lithium but the high-temperature 
end is of significantly different shape. 


INTRODUCTION 


Up to 70% magnesium may be dissolved in lithium without changing the 
body-centered cubic structure of the pure parent metal. This alloy system 
therefore provides a means of studying the effect of a solid soluble ‘‘impurity” 
on the physical properties of lithium metal. 

Zener (1947, 1948) showed that the rapidly increasing free energy with 
decreasing temperature for body-centered cubic metals and alloys might lead 
to a solid-phase transformation. Such a transformation in pure sodium and 
lithium and some lithium-—magnesium alloys was discovered and studied by 
Barrett and co-workers (1947, 1948, 1950, 1956) using X-ray and metallo- 
graphic techniques. More recently measurements have been made on the 
electrical, elastic, and thermal properties of sodium and of lithium of two 
different isotopic compositions (Dugdale and Gugan 1958, 1959; Verdini 1959; 
Martin 1958a, 6, 1959a, 6) in which effects due to the transformation were 
found. The interpretation of the electrical resistance data on lithium was 
made difficult by the effects of small impurity concentrations on the residual 
resistance and measurements were therefore extended to dilute alloys of 
magnesium in lithium in order to investigate the effect of controlled amounts 
of impurities (Dugdale and Gugan 1959). The present work examines the 
effect of a small addition of magnesium on the thermodynamic properties of 
the transformation. 


EXPERIMENTAL 


The sample was prepared from the “‘sodium-free”’ grade of lithium obtained 
from the Lithium Corporation of America (principal impurities Na, 0.005%; 
K, 0.01%; Ca, 0.01%; N, 0.02%; Fe, 0.001%; atomic weight = 6.945) and 
“spectroscopically pure’’ magnesium obtained from Johnson Matthey. A 
master alloy of approximately 10 at.% magnesium was prepared. by heating 
the lithium and magnesium to 800°C in a low carbon steel crucible. The 
master alloy was then diluted with lithium in a stainless steel container to 
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make an alloy containing approximately 1 at.% magnesium. The calorimeter 
sample was extruded from the diluted alloy. The exact composition of the 
sample was determined as follows. The residual electrical resistances at liquid 
helium temperature of pieces taken from the ends of the sample were com- 
pared with the residual resistances of alloys for which the composition had 
been determined gravimetrically. The composition was found to be 0.95+0.05 
at.% where the tolerance indicates both the uniformity of the specimen and 
the estimated absolute accuracy of the measurement. From this the ‘‘average’”’ 
atomic weight is 7.110. The weight of the sample, after buoyancy correction, 
was 4.451 g. 

The measurements were made in a continuous heating adiabatic calorimeter 
(Martin 1959c). The estimated accuracy of the measurements is +0.3% above 
80° K, decreasing to +2% at 20° K, after allowing for the uncertainty in the 
exact composition of the sample. 


RESULTS 


Provided that the alloy is not cooled below about 90° K the body-centered 
cubic structure is unaltered; smoothed results of measurements under these 
conditions are given in column 2 of Table I. On cooling the alloy below this 
temperature a partial transformation to a faulted hexagonal close-packed 
phase occurs, the amount transformed increasing as the temperature is 
lowered. On heating such a partially transformed sample a specific heat 
anomaly is observed in the region where reversion to the high-temperature 
phase takes place. Specific heat results obtained after cooling the sample to 
20° K are given in column 3 of Table I for a sample annealed at about 300° K 
and in column 4 for a sample not annealed following a previous cycle of 
transformation and reversion. 


DISCUSSION 


The results for the body-centered cubic (b.c.c.) phase will be considered 
first. The “calculated” results given in column 5 of Table I have been obtained 
from specific heat results on lithium (Martin 1959a) and magnesium (Craig 
et al. 1954) assuming the Kopp—Neumann rule (i.e. that the thermal capacity 
of the alloy equals the sum of the thermal capacities of its constituent parts). 
The experimental results are some 2% lower than the calculated results and 
are also lower than the specific heat of pure lithium. This situation is not at all 
uncommon in alloys and it is interesting to seek a physical explanation. 
Assuming that small quantities of magnesium dissolve in lithium in a divalent 
form and that there is no change in band shape then the electronic specific 
heat of the alloy would be higher than that of pure lithium. Hence the observed 
reduction in specific heat is most probably due to a change in the lattice 
contribution. On a simple Debye model the Debye temperature would be 
lowered by the reduction in elastic constants on alloying (Verdini 1959) but 
it would be increased by the accompanying decrease in lattice parameter 
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TABLE I 
Specific heat of a lithium — 0.95 at.% magnesium alloy 


After cooling to 20° K 





Not Calculated by Kopp- 








Temp., °K B.c.c. phase Annealed annealed Neumann rule 
21 _— ke ~ 
25 — .173 — 
30 — . 285 
35 — .424 — 
40 _ 587 — 
45 — . 767 — 
50 — .963 —_ 
55 — 1.170 — 
60 — 1.388 — 
65 ans 1.608 = 
70 — 1.823 — 
75 — 2.040 — 
80 — 2.254 — 
85 - 2.465 — 
90 — 2.667 2.667 _ 
95 — 2.859 2.864 - 

100 3.026 3.045 3.052 3.107 
110 3.329 3.552 3.665 3.417 
120 3.611 4.047 3.969 3.707 
130 3.870 4.005 3.930 3.963 
140 4.098 4.141 4.118 4.195 
150 4.298 4.316 4.305 4.397 
160 4.479 4.482 4.480 4.581 
170 4.647 — — 4.755 
180 4.788 - — 4.900 
190 4.918 — — 5.036 
200 5.040 — — 5.153 
210 5.149 ~ — 5.256 
220 5.243 — — 5.350 
230 5.342 = —_ 5.438 
240 5.428 _ 5.521 
250 5.499 — — 5.601 
260 5.579 -- — 5.680 
270 5.654 — — 5.752 
273.15 5.675 — — _- 

280 5.722 — — 5.822 
290 5.785 — — 5.889 
298.15 5.833 —_ —_ 5.941 
300 5.843 — —_ 5.952 





Note: Atomic weight = 7.110; J = 4.186. 


(Herbstein and Averbach 1956). Insufficient detail is available to make a 
quantitative estimate of the first of these effects. The second would be negli- 
gibly small in comparison with the observed change (at 100° K) from a Debye 
temperature of 393° K for pure lithium to 402° K for the approximately 1% 
alloy (where these values have been calculated after subtraction of the elec- 
tronic term but without correction for thermal expansion). The effect must 
therefore be due to the change in the vibration spectrum on alloying. 
Turning now to the experiments relating to the martensitic transformation, 
Fig. 1 shows the experimental results for a typical anomaly caused by reversion 
to the high-temperature phase. The shape of this anomaly and another, 
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Fic. 1. Thermal capacity of calorimeter containing 4.451 g lithium—0.95 at.% magnesium 
alloy. 


O body-centered cubic phase; 
A after cooling to 20° K (heating rate 20° K/hr). 


obtained under different conditions, are shown in Fig. 2. The heat content of 
the anomaly observed after cooling an annealed sample to 20° K is 9.3 cal/ 
g-atom. For pure lithium under similar conditions a figure of 9.4 cal/g-atom 
was observed (Martin 1959a). It seems probable, therefore, that the presence 
of a small percentage of magnesium in lithium does not significantly alter the 
course of the martensitic transformation. However, it is already known that 
the temperature of equilibrium between the two phases is increased with the 
addition of magnesium and this is reflected in the position of the reversion 
anomaly. Figure 3 shows the progress of the reversion with increasing tem- 
perature and was obtained by integration over the anomaly shown in Fig. 1, 
thus assuming that equal areas of anomaly contribute equally to the reversion. 
This may not be so, as discussed elsewhere (Martin 1959a), due to the effects 
of stored elastic energy and cold work. Figure 2 also shows an anomaly obtained 
after cooling to 20° K a sample which had been previously transformed and 
then heated to 170° K, i.e. the sample was in a cold-worked state due to the 
previous transformation. The heat content of the anomaly is slightly reduced 
and the peak shifted to a lower temperature, but the general shape remains the 
same—in contrast to the corresponding results for pure lithium (Martin 
1959a). On the whole it might be said that the shape of the anomalies observed 
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Fic. 2. Specific heat anomalies in lithium — 0.95 at.% magnesium alloy due to martensitic 
transformation. 


After cooling to: Heated to: 
20° K 170° K 
20° K 300° K 


100 


Fic. 3. 


Measurements made in this sequence. 


8.9 





20*K/HR HEATING RATE 


130 140 150 


110 120 
TEMPERATURE °K 


Progress of the reversion in lithium—magnesium alloy. 


160 


Jac,.d T, cal/g-atom 
9.3 











30 CANADIAN JOURNAL OF PHYSICS. VOL. 38, 1960 


for this alloy are more like those of sodium than those of lithium (Martin 1959a, 
b). In particular there is no sudden ending of the anomaly at the high-tem- 
perature end, but there are some time-dependent effects in the specific heat, 
as shown in Fig. 4, from which it will be seen that a change in the heating rate 


14.0 





= 
N 
43:07 
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TEMPERATURE °K 
Fic. 4. Time-dependent effects in the reversion region. 

Runs at 20° K/hr heating rate: Runs in which heating rate is changed at 145° K: 
V body-centered cubic phase, body-centered cubic phase @ 20° K/hr, 
A after cooling to 20° K. O 2° K/hr; 

after cooling to 20° K @ 20° K/hr, 
0 2° K/hr. 


alters the shape of the anomaly. The results of a check experiment are also 
given in Fig. 4, showing that a change in heating rate does not significantly 
affect the results obtained with the b.c.c. phase. 
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CONCLUSION 


The results would be consistent with the fact that the addition of a small 
amount of magnesium to lithium does not significantly alter the course of the 
martensitic transformation, apart from raising slightly the temperature of 
equilibrium between the phases. There is some difference between the alloy and 
pure lithium in the final stages of the reversion to the high-temperature phase. 
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COMPARISON OF FLOW LINES IN VARIOUS TYPES OF 
RHEOLOGICAL BODIES! 


Sistr CHANDRA Das 


ABSTRACT 


An attempt has been made to investigate the nature of flow in Bingham bodies 
compressed by means of parallel rigid approaching plates. The results have been 
compared with those of plastic and viscous flow under the same conditions. It is 
shown that the flow lines in all three cases are very similar. 


INTRODUCTION 


In recent years attentions have been focused on Bingham bodies as an im- 
portant case of rheological substances, but not much work has been done on 
its theoretical aspects. The case of Bingham flow through a tube has been 
discussed by Buckingham (1921) and Reiner (1926), while Reiner and Riwlin 
(1927) have done the integration in the case of a rotating cylinder. In a series 
of papers Oldroyd (1947a and 1948) has obtained the solutions for Bingham 
bodies for cylinders of various shapes moving through each other. He has also 
discussed the two-dimensional case (19475) in some specified instances. 

In the present paper an attempt has been made to investigate the relation 
of the Bingham flow in two dimensions in the case of approaching parallel 
plates with other types of flow (plastic, viscous). A yield condition similar 
to that valfd in plastic flow has been assumed, and the flow lines have been 
drawn. The result is that the flow lines for a Bingham substance are very much 
similar to those for plastic and viscous materials. 


GEOMETRY OF EXAMPLE. YIELD CONDITION 
We consider the case of a plate compressed between planes y = a and 
y = —a; there is no motion in z-direction. Therefore the principal stress in 
the z-direction is 1/2(¢:+¢2) where o; and a2 are the principal stresses in the 
xy-plane. The criterion of yield is taken as 


(1) 01—02 = 2p 


where » is the yield stress. It may be noted here that equation (1) satisfies 
both Tresca’s and von Mises’s yield criterion. 

If oz, oy, Try are the components of stress referred to x and y axes, the equa- 
tion of equilibrium can be written in the absence of a body force as follows: 


00 Or 
wae ae ie 
Ox + ay 0, 


(2) 
Otay 4 Oey _ 
Ox tay = 0. 


'Manuscript received September 22, 1959. 
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The yield criterion (1) becomes 
(3) 1/4(o,—0y)? +1iy = v. 
Taking the plates to be rough at the surfaces of contact and ruling out the 


possibility of their transmitting a shearing stress in excess of the yield stress 
in simple shear, we can formulate the boundary condition as follows: 


(4) Try = sky at y = ¥a. 


The stresses o, and a, can be easily eliminated from the equations (2) and 
(3), and we get an equation for rz, which is 


(5) Orr Sos Or, ‘nd +2 ae: (v?—7? y 
“Ox? ay axdy cit 


A solution of this equation which satisfies the boundary condition (4) is 


saa 


(6) Tey ™ a 
The other stress components in this case are 
% vx ( ”) 
or = —pt+—+2v\1—-4), 
(7) a a 
vx 
ore 


where p is a constant. 
The stress deviators (denoted by s) can be written as 


2\} 
(-2) 
a 


Ss: = 2 
2\4 
(8) Sy = =ihy%) 5 
Try = ae ‘ 


a 


SOLUTION FOR VARIOUS TYPES OF FLOW 
The solution up till now has been quite similar to that in the case of a narrow 
strip of plastic material if it is squeezed between two rough rigid plates which 
are kept parallel to each other. We now propose to find the velocity field in 
the case when the material behaves like a Bingham substance, as the flow 


starts after the yield stress is reached. 
Writing U and V for a and #, we reformulate the equations (8) in the case 


of a Bingham substance as follows: 
2\4 
aU 
=" (1-%) “30 Ox 


2\4 
~o(i-2)' = 265 


dy 


(9) 


Sy 
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where ¢ is a function of the co-ordinates; the bulk viscosity term has been 
neglected. Further, 

P 
(10) ty = 2 = Frty( 22420) 
according to whether both y and [(dU/dy)+(dV/dx)] 2 0. Here 7 is the 
plastic viscosity and y, 7 are both taken to be positive. 

Here we have retained the stress-strain relations for the normal stress 
deviators as in the case of plastic flow. Only for the tangential stress we make 
use of the Bingham law which may be taken as an extension of the law de- 
scribing a Newtonian fluid with the addition of a yield stress. As we are only 
interested in the study of flow in a Bingham fluid arising out of shear stress, 
it may be assumed that the bulk viscosity is negligible in comparison with 
the shear viscosity. We have also retained Tresca’s or von Mises’s yield con- 
ditions for plastic flow since we do not have sufficient experimental evidence 
to indicate the contrary. The yield condition arises from simple common 
sense when the elastic limit is exceeded and the matter begins to flow or tends 
to flow. In the case of a Bingham substance where some kind of flow also 
exists, it seems natural to assume that the same condition will hold. Inte- 
grating equations (9) and (10) we have 


ws yVo 
a 
,_ Vox 2 f) : 
(11) U= : ae ay—5 ify > 0, 
Vox v ( »*) 
o. ee ee tS . 
U= — ayt+5 ify <0, 


if the plates are approaching with speed Vo. 
The flow lines for the above distribution of velocities turn out to be 


3 

tS] 3 yok ] = const. for y > 0, 
(12) ; . 
5 2 y° 

xy te 3y te = const. for y < 0, 


where, following Oldroyd (1947), we take S as a dimensionless parameter 
depending on the material and velocity of the boundary as well as the distance 
between the plates which is given by 


(13) oes 


It is to be noted that to have the boundary condition as specified in the 
problem, some type of arbitrary normal pressure should be applied in the 
boundary. This has been roughly shown in Fig. 1. 








DAS: FLOW LINES IN RHEOLOGICAL BODIES 35 











We now draw the lines of flow given in (12) with S = 3, and obtain the 
result shown in Fig. 2. In the same condition, the flow lines for a plastic 
material would be as shown in Fig. 3. Finally, for the sake of comparison, the 
lines of flow have also been drawn in the case of Newtonian fluid and are 
shown in Fig. 4. 





Fic. 2. Bingham flow. 
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Fic. 3. Plastic flow. 
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Fic. 4. Viscous flow. 
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DISCUSSION 


From the diagrams of the flow lines in various cases it is seen that the 
latter for a Bingham body differ very little from those in a plastic material. 
The same conclusion was arrived at by Oldroyd from his calculations of the 
flows within different types of cylinders. It seems that this small difference 
of flow patterns in the various types of rheological behavior is a general 
result. This follows from a consideration of the boundary and yield conditions 
in the case of plastic and Bingham flow. In the above calculations it turned 
out that these are the same in both these cases, and therefore the separation 
of flowing and nonflowing parts of such rheological bodies is also the same 
in the two cases. Equally, the boundary stresses are the same and the only 
difference lies in a slight change of the shape of the flow lines. 
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A NOTE ON DIFFRACTION BY AN INFINITE SLIT! 


R. F. MILLAR 


ABSTRACT 


The two-dimensional problem of diffraction of a plane wave by a narrow slit 
is considered. The assumed boundary values on the screen are the vanishing of 
either the total wave function or its normal derivative. In the former case, a 
differential-integral equation is obtained for the unknown function in the slit; in 
the latter, a pure integral equation is found. Solutions to these equations are given 
in the form of series in powers of e (where e/z is the ratio of slit width to wave- 
length), the coefficients of which depend on loge. Expressions are found for 
the transmission coefficients as functions of e and the angle of incidence; these 
are compared with previous determinations of other authors. 

A brief outline is given for the treatment of diffraction of a cylindrical wave 
by the slit. 


INTRODUCTION 


The problem of diffraction of a plane wave by the infinite slit separating 
two half-planes has received considerable attention in recent years. The extreme 
cases in which the wavelength is much less than, and much greater than, the 
slit width have been treated in some detail. In the latter circumstances, with 
the exception of the work of Morse and Rubenstein (1938), Hoop (1955), and 
Nomura and Katsura (1957), little attention appears to have been given to 
other than normal incidence. 

It is the purpose of this present note to extend previously given results for 
the long wavelength region to oblique incidence in two dimensions. The 
method employed is that of Bouwkamp (1953) and Miiller and Westpfahl 
(1953). No new mathematical difficulties are encountered, although the alge- 
braic manipulations are somewhat more tedious. 

Expressions are obtained for the first few terms in the series for the trans- 
mission coefficients when the wave function or its normal derivative vanish 
on the screen. A comparison is made with the results of the Mathieu-function 
analysis of Morse and Rubenstein (1938), the calculations of Hoop (1955), 
and the short wavelength asymptotic solution (Millar 1958). The diffraction 
of a cylindrical wave by the slit is briefly considered. 


FORMULATION OF THE PROBLEM 


Let the diffracting screen occupy the plane z = 0 in a rectangular Cartesian 
co-ordinate system (x,y,z). The edges of the slit (of width 2a) are parallel to 
the x-axis, with the origin at the center of the slit. The propagation vector 
of the plane wave e* °8¢ + = 8") (incident from z < 0) makes an angle a with 
the positive y-axis, and lies in the y-z plane; thus all field quantities are 
independent of x. The total wave function will be denoted by ¢, the notation 
employed being essentially that of Bouwkamp (1953). The time-dependence 
factor e~‘* is suppressed throughout. 


1Manuscript received August 4, 1959. 
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(1) @ = 0 on the Screen 
If ¢ = 0 on the screen, then 
a C cos a+z sin a) __ ,iky cos a—z sin “+. ¢,(y,—3), 
$1(9,2), 
The function ¢$1(y,z), defined only for z > 0, is given by 


ay Ww 
VA 
oS 


() aye) = 552 fF a0. OHPRVIO-y)*+2') ay’ 
Since 


iky cos a 





= iksinae 


0 
ae $1(y,2) Be 


when —a < y < a, the following differential—integral equation is obtained for 
the determination of ¢:(y’,0): 


@) (+5) J ¢1(y',0)H5” (kly—y'|) dy’ = 2k sin ae” °°8* 


when —-a<y <a. 

Equation (2) may be integrated, and the two arbitrary constants deter- 
mined in terms of the unknown function by consideration of the integral, 
and its y-derivative, at y = 0. Since ¢1(—a,0) = ¢:(@,0) = 0, it is found that 


JF 610070) 228° Cely—y') cos bye (eIy') +B HME’) ay 
= 2(e%" *—cos ee asin ky)/(k sina), 


where —a < y < a, and the integral involving the Hankel-function derivative 
is to be interpreted as a Cauchy principal value. When a = 32, ¢1(y’,0) is 
an even function of y’, and eq. (3) reduces to the form obtained by Sommer- 


feld (1954). 


(it) d6/dz = 0 on the Screen 
If d¢/dz = 0 on the screen, then 


5g cos a+z sin O16 tk(y cos a—z sin a) hake ake z <0, 
~ \g2(y,2), z>0, 
where ¢2(y,z), defined only for z > 0, is determined by 


(4) daly) = 24 [33 2600) HP AVI-9')*+8°I) dy’ 


The equality of ¢2 and e“” °*¢ in the slit implies that 


(5) Se 5q! $2(9'0) Ho” (kly—y'|)dy’ = 21 et” 8 
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SOLUTIONS OF THE INTEGRAL EQUATIONS FOR ka < 1 
Because the determination of (0/02’)$2(y’,0) is rather more readily accom- 
plished, eq. (5) will be considered first. The substitutions y = acosu, y’ = 
acosv, ¢ = ka transform eq. (5) into 


(6) f F2(v) H$” (elcos u—cos v|) dv = 2i et 82 8 
v0 
where 
5 ; 0 
(7) F,(v) = asinv = ¢2(a cos 2,2) ; 
Oz z=0 


andOQ<ucr. 
For small e, H{”(e|cos u—cos v|) may be written as 


(8) H3” (elcos u—cos v|) = =! (Wot vie’ + woe’ + vse + ...) 
where 
Yo = pt+log|2(cos u—cos v) |, 
(9) (-1)" a 1 
| vn = 2 nl? (cos u—cos v)"[p—1—3... —7 +log|2(cos u—cos v) |] 
for x > 1, 
(10) pb = log tve—}71, 
and log y = 0.5772157 . . . is Euler’s constant. 


A solution is assumed of the form 
(11) F,(v) = fothietfoer+ Beis 


where the f, depend only on », a, and log e. (For normal incidence, fon41 = 0.) 
If the above expansions are inserted into eq. (6), the equating of like powers 
of ¢ leads to the following infinite set of integral equations: 


|| fovo ds = T, 
0 


oF 

J fivodv = im cosa cos u, 
| 0 
| 

! 7 x 
| 2 2 
| J fevode = —%3mCOS acos u— J fos de, 
| 


(12) | | favo do = alee cos" a cos! u— f fava a, 
| 0 0 
| 


| 7 7 7 
J favo dv = = cos‘ a cos‘ u- f fobs dv [ fovs dv, 
i Jo 0 0 
| ar . 7 
J fovo dy = za cos* a cos® u- f favs do— f five dv, 
| J90 0 0 


Pe ; 
Js dv = gan cose costu— f fas ao— J favs do— J for dv. 








one eee NE 


ee 
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If it is now assumed that 


Son = > bon, » COS 2vv, 
yen0 


(13) . 
[fants = ai bony1,» cos (2v-+1) », 


41 


ts QBs 5s 


with 6,,, independent of v, then, since ¥o possesses the expansion 


c@ 
COS MU COS MV 
p—2 S> Sosmucos me 


(14) 
m=1 m 


Yo = 


it is found, 


after the evaluation of certain integrals (see Appendix), that 


{ - 3 
1 ‘ 
Jeo. = p’ bi.0 = —1COSQ@, boo = - 
| 1 1 i i 
es = 5 cos” a- 14 i), b3.0 = 3 cos” a—7 (p+3) cosa, 
a = -—cos’ ae oO: 
| 21 = 8 a 16 © Sa, 


2), 


Bei ol 208 (yn) SS pa hhh ee 
[one 64 2p/ 32 *128\p? "2p 


|, _ costa € ian 2s 2), 
meee 24 T\3t a5 16 +16 \3 2p 











| 3 
b42 = _ 2.3 0 7 ee 24 (1 (1 +3) ’ 
| jcosta | i 1,13 
| bs,.0 = _* 1 Ae (p+4) cos’ a—Fa (7-1)438 COs a, 
(15) ee i 4 fom ay (p+ #) cos a, 
5. . < —if0s'a , icos' a , i cos a 
sy 384 768 3072’ 


6 2 
cos a l cos’ a 3 3 4 cos a@ 
° ~~ 2304p ( 1) 384 +(1 Ma 2p°/ 768 


ik E 
+763 1 


6 4 2 
cos @ 1 ,1)\cos a 49 ‘) cos a@ 
be. = 768 -(144) 512 +(r-4 4125) 256 


20.3 
2, 3.) , 


1 
~ 512 (»-! 8, 5), 


be ¢ = £08 @_COs* a _ (i 8s ok (i. *) 
; * «9601920 \p 2048 2048 \p 5/’ 
cos’a@ cos’ a cos" et ee 

3840 7680 30720 73728p 61440° 





| bes = 
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Equation (2) (or eq. (3)) for ¢; may be solved in a similar manner. If 


(16) eF,(v) = ¢1(a cosv,0), 

and a solution is assumed of the form 

(17) Fi(v) = Zotgietgee?+..., 
where 


i 
| 


Zn = es An,» sin (2v+1)v 
(18) j = 


Ba 3 
qsin a(p—3+4 cos’ a), 











OR 
os “ 

a2, = ig a eee +cos’ a), a3,0 = 2 —. 

sina cosa ,, 
131 = _sinaoma ($+cos" @), 
| i sin @ 5 ) 7 cosa, cos* #1 
d4.0 = esp ~ 2( 84 cos? a +i¢§- 3 + 12 |" 
zsin a cos’ a cos’ a 
eS ae oe —— a, ia 3 .), 
| isina 6 ee a ) 
4.2 = — F100 (1440s at3 cos @ }, 


| 
| sl °OS O23 2 1 
(19) } 4,0 ees (eye cos a—3 cos* «) . 


512 


ll 


| SiN @ COS a@ 2 4 

125.1 = 7680 (3—cos a-—2 cos a), 
Sin @ COS @ 2 4 

Gs.2 = aaa (+4 cos a+cos a), 


oe 
a Steal ag 1p Fe a 607 4 (o Se) it 
p 19 oe | 
+(2- it) cos'at 35 | 
_ tsina /* a ‘ee cos 9) 
| 26.1 = “519 [o*-2 g +96 20.4(4 go) 008 @+ Gg +39 | 
isin @ 


64 Pe 
a6.2= ne sp 243 = COS tate cos ‘ats cos® a} 

















| hei 
1 sin a 
| 26.3 = ee. cos” ats cos‘a+cos .) ; 





n 
| . 
/Zentt = Dy Gentry Sin (2v+2)y, o= 0,1,2..., 
\ v= 
it is found that 
{ 
|do,o = —tsina, d1,0 = }sinacosa, 
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FAR-FIELDS AND TRANSMISSION COEFFICIENTS 


In the following, attention will be confined to the far-field region. However, 
it may be noted in passing that, for /(y?+2°) > a, it is legitimate to expand 
the Hankel functions appearing in eqs. (1) and (4) in powers of ¢, and so 
representations of ¢1, ¢2 in this region may be found. 

Let a point in the far-field region be specified by polar co-ordinates (9,8), 
where @ is measured from the positive y-axis. Then, according to eq. (1), 


ett ett? 
20 ; 
( ) gily z)~ V/ (27) Tt (Rp) 
or, in terms of the quantities introduced in eqs. (16), (17), and (18), 


‘7 
. . —i 6 
€ sin of sin v ¢;(a cos v,0)e “* °° “°° "do, 
0 





Qn\} We ie 
(21) o1(y,2) ~ (2*)'. ae > Tinvila,O;p)e"*”, 
where 
| T}? = 7 40,0 sin 6, 
| 
6 
| ) = out (i 0,0 cos” 9— =" cos i ; 
| 4(1) sin@} 7 
| QI5 =- 4 92 do,9 COS* 6— 2 cos £ (anotae 1) cos” 0 
(22) |; +"3" cos 0+ ia, | ; 





| sng} t 41.8 coc? 
| 77 4 9216 20.9 C08" 4 758 © 0 


1 d2,0 , 42,1 4 1 3,1 2 
-é(? 3 *4833) con: Cog (1.423!) 00s : 


1 \ 2 (5.4 ° 
+3 (a4,0+a4,1) cos’ A— =o 0~ ine | ; 


The plane-wave transmission coefficient t" (a), defined to be the ratio ol 
the energy transmitted through the slit per unit area to the energy flux in 
the incident field per unit area normal to the direction of propagation, is 
given by 


(23) t (a) —Fim| > TBslaaipye™ | 
n=(0 
so 
) 8. 4 2 \,e/ 5 
(24) ¢'(a) = 35 sin’ aylt+2 e( 1~70~F con a)+£| 5 cos 6 


17 1 ; 2 
-——6§ ) cos or 109 —3365-+2886°—240°) |+... ¢ 
32 192 

where 6 = log }ye. 

Note: The transmission coefficient may also be obtained by a direct calculation of the 
energy flow across a semicircular cylindrical surface in the far-field region. Here ¢; is of order 
é as « tends to zero, and this method of calculation, together with equations (22), will yield 
one more term in the series for ¢ ‘)(@) than appears in eq. (24). 











44 CANADIAN JOURNAL OF PHYSICS. VOL. 38, 1960 


In a similar manner it may be shown that 


Be. Sa \* sites ‘ 
(25) 2(y,2) ~ (25) _ DL, Ti? (a,oip)e™, 
where 


i 1 
0 a bo.0, 


7?) = 70.9 co5 sets ail: 3 b2,0, 





8 
| (2) bo,0 coe 6 
T4 = ~ 728 cos ae tial 0+ (be, ot} be, 1) 
(26) ; 
+7 bs,0 cos 6— 3 54,0, 
77(2) ee boo 1b1.9 5 — - 
Ts = 469g COS O+79g COS O—\ 5 b2.0+b2.1) Top 
an — 
—1(3b3,0+03,1) —— ~ + (2b4,04+04,1) cos” wet bs,0 cos 0— } be. o. 
The transmission coefficient, ae is given by 
(27) t(a) = — Im 7. 72 (aaip)e” 
n=0 


so 


c 2 
(28) ta) = F }1-£ cos 2a+=>> : él +5 (q—3) cos’ a+8 cos tid 


q 
dae wn 2 5 5 
= {| ee S40 (y-2) et (o-§) £0845 cog haley 


with g = m’°+46’. 


The first two terms of eq. (24), and the first three terms of eq. (28), have 
been found, for a = 3x, by Bouwkamp (1953), and Miiller and Westpfahl 
(1953); see also Bouwkamp (1954). Nomura and Katsura (1957) have found, 
for normal incidence, the first two terms of eq. (24), and the first term of 
eq. (28). Hoop (1955) has calculated ¢‘ (a) to the same order in € as eq. (24), 
and ¢®)(a) up to and including terms of relative order e‘. 

The variation of the transmission coefficients, ¢ (a) and ¢ (a), with a and 
ka is exhibited in Figs. 1 and 2. For ka < 1, the curves are determined by 
eqs. (24) and (28); for larger ka, the asymptotic values (Millar 1958) are 
plotted. As a standard of comparison, the curves of Morse and Rubenstein 
(1938) for a = 0° and 90° have been included. 

It is apparent that eq. (24) for (a) is a good approximation for ka < 1; 
agreement is less satisfactory in the case of ¢® (a). For fixed ka, the accuracy 
of approximation in both cases appears to decrease as a decreases from 90° 
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a) 
¥ (a) 











Fic. 1. Transmission coefficient ¢™ (a) as a function of the angle of incidence, a, and ka. 
————— Present theory for small ka; asymptotic theory for larger ka. 
—--— Theory of Morse and Rubenstein for a = 90°. 


and, in the so-called ‘“‘resonance region’’, near ka = 1, neither the present nor 
the asymptotic method appears to be particularly well suited to the analysis 
of the problem. 

The question of the convergence of series such as those which appear in 
eqs. (24) and (28) is still a matter for conjecture. Because of the presence of 
log ka in the coefficients of powers of ka, it is possible that the expressions 
given above are only asymptotic as ka — 0. 


DIFFRACTION OF A CYLINDRICAL WAVE 


Diffraction of a cylindrical wave by the slit may be treated in the same 
manner as the case of an incident plane wave. Suppose a line source whose 
location is specified by polar co-ordinates (r,p) radiates, from z < 0, a field 
H) (kR), where R is the distance from the source. Then the equations obtained 
are once again (2) and (5),with the right-hand sides replaced by 


2i < ‘ (RR) and 24 HS(kR) 


z=0 z=0 
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0 2 4 6 8 1.0 12 14 1.6 1.8 2.0 
ka 


Fic. 2. Transmission coefficient ¢'?’(a) as a function of the angle of incidence, a, and ka. 
— Present theory for small ka; asymptotic theory for larger ka. 
—-- Theory of Morse and Rubenstein for a = 0° and a = 90°. 


respectively. These functions may be rewritten as 27kr sin @ H‘P)(RS)/S and 
2iH'» (kS), with S? = r°+y?—2ry cos @ and -a < y <a. If r > 4, it is per- 
missible to expand each function as a convergent power series in € (= ka), 
the coefficients of which depend on kr, cos ¢, cos u(= y/a), and the analysis 
proceeds as before. For example, 


H(kS) = H?(kr)+cos ucos ¢ Hi (kr)et ... 
and, if kr > 1, 


2 4 cecil 
H\?(kS) ~ (2 ) e*'-**"(1 ji cosucos d-e+...). 
kr 


a } 
27 lim [= e *tr (RS) 


kro ~ 


is formed, and ¢ is equated to +a, then the power-series expansion of the 
right-hand side of eq. (5) is obtained. 
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APPENDIX 


EVALUATION OF CERTAIN INTEGRALS 
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THE RADAR CROSS SECTION OF A SEMI-INFINITE BODY! 


H. Brysk? 


ABSTRACT 


The concept of cross section as applied to a semi-infinite scattering body seems 
to require some clarification. The need for careful formulation of the problem 
arises because of the simultaneous occurrence of two characteristic lengths tending 
to infinity: the range from the radar to the target, and the size of the target. The 
infinite range assumption in the definition of the cross section allows the incident 
wave to be approximated asa plane wave in the case of a finite scatterer. For a 
semi-infinite body, it is customary to retain the plane-wave incidence, and 
introduce ad hoc arguments to dispose of the awkwardness due to the infinite 
extent of the scatterer. A return to the basic definition of a cross section, and 
examination of its motivation, lead here to an unequivocal formulation for the 
cross section of a semi-infinite body. Its consequences are pursued in the physical 
optics approximation. In particular, the nose-on backscattering from a body of 
revolution is exhibited, and results are computed for the paraboloid and the cone 
(which turn out to agree with the traditional ones). The broadside backscattering 
from a cylinder is also calculated, and the difference in this case between mono- 
static backscattering and the return in the backward direction from an incident 
plane wave is discussed. 


INCONCEPT 


The return backscattered from a radar target is expressed formally by the 
range equation 


(1) a P,/P, = GG, o/(4r)*R! 
where P,; = power scattered into the receiver, 
P, = power transmitted, 
G,;. = gain of receiver, transmitter, 
r = wavelength, 
R= range (minimum distance from radar to target), 
o = radar cross section. 


(The discussion will be restricted to monostatic operation in order to simplify 
the mathematical manipulations; the generalization to bistatic operation does 
not involve any additional conceptual problems.) Equation 1 can be applied 
with great generality provided that the possibility is admitted that the gains 
and the cross section vary with the range. It is most commonly used to describe 
scattering by a finite target at very long range, in which case the gains and the 
cross section become range independent and the factors in equation 1 separate 
naturally as distinct physical entities. The concept of a cross section is usually 
referred to the long range limit. From the above definition, it then follows that 
the cross section is 
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(2) o = lim 4rR’|S*/S'| 
Raa 
where 
S*:' = scattered, incident Poynting vector 


or, using the relationship between the electric and magnetic field vectors in 
free space, 


(3) o = lim 42R*|H*/H'|’ 
Roa 
where 
H*:' = scattered, incident magnetic field. 


It is customary to assume that the incident field is a plane wave; this 
amounts to removing the source to infinity. In fact, electromagnetic boundary- 
value problems are usually formulated for an incident plane wave. Instead of 
the limit process indicated by equation 3, what is effectively done is to go 
separately to the limit R— © for one factor in the expression (namely 
|H'|?), and later carry out the limit for the remaining factors, as if there were 
two distinct and independent limits involved. This procedure simplifies the 
mathematical manipulations, and it can be justified provided that the range 
is sufficiently long compared with all the characteristic lengths of the problem 
(wavelength, maximum linear dimension of target, ratio of projection area of 
target to wavelength). 

The above approximation breaks down when one dimension of the body has 
to be considered infinite. This need not impair the usefulness of the cross- 
section concept, although it does necessitate greater care in setting up the 
computation. The usual practice in this situation (for example, Kerr 1951) is 
to retain the plane-wave incidence, and to formally cut up the body into the 
analogue of Fresnel zones, each section then being considered as a finite body. 
A qualitative argument is then made for ignoring all but the first section (on 
the basis of approximate destructive interference from neighboring zones) and 
for discarding most of the contribution of the first zone (ascribing to end 
effects those terms that depend on a dimension of the zone and claiming that 
they are cancelled by corresponding end terms from the next zone). This 
approach has some intuitive appeal, but the reasoning is a bit inelegant. 

In reality there is no such thing as a semi-infinite body. The problem can be 
clarified by going back to the physical situation that has impelled the imputa- 
tion of an infinite dimension to the body. The motivation is that the body does 
not terminate sharply (for example, it has a jagged edge, or a blurred ending 
such as an ionized region) and the boundary surface is mathematically extended 
to infinity to avoid undesired end effects that would arise from a spuriously 
clean edge. The source of radiation and the observer have actually a well- 
defined position. With respect to this position, the boundary extended to 
infinity is receding from the radar. It is in the spirit of the infinite body 
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approximation that this distance variation be taken into account, in the same 
manner for the incident wave as for the scattered wave. Removing the plane- 
wave approximation for the incident wave necessitates the stipulation of the 
antenna pattern—which is reasonable. Thus, the limit process should be 
handled in a single step as per equation 3, with the incident field represented 
not by a plane wave but by a spherical wave with a specified angular variation 
of amplitude. The R — @ limit can then be consistently attained for a semi- 
infinite body, as it need only mean that the range is much longer than the 
wavelength. 

Formally, the problem is now set up in the standard way (for example, 
Kerr 1951) except for omission of the plane-wave approximation. The field 
scattered by a perfectly conducting body is given rigorously by 


(4) H* = —(1/47) J Gxpxo(rte*nas 


= unit vector normal to the surface S, 


where n 
H = total magnetic field at the surface, 
r 
k 


= range from. the radar to a point on the surface, 
= 27/X. 


If kyr > 1 forall ron S, 


(5) V(r-! e- #?) = —F ir! e- #F 

so that 

(6) H* = —(ik/4n) f *X(AXH) re *' dS. 
Ss 


Equation 6 is the starting point for the physical optics approximation. In the 
latter, the field at an illuminated point on the surface is set equal to twice the 
tangential component of the incident field 


(7) H = 2(H'—(H'-n) a] 

(as would be the case if the body were replaced by its tangent plane at the 
point), so that 

(8) ?X(”#XH) = —2H!' (7-n), 

while the field at a surface point in the shadow is set equal to zero. The depar- 
ture from the usual formulation comes in not proceeding now to the further 
approximations that r~ be factored out of the integral in equation 6 and that 
H' be taken as a plane wave. Instead, H' will be considered as a spherical wave 
with a specified angular amplitude diagram. As an example, the field due to an 
electric dipole oriented along the x axis will be used in the illustrations below: 
(9) H' = Hy (# XX) r— e- ®" 


leading to 


(10) H® = (ik/27) Ho f (#X*) (Fn) re P* dS. 








BRYSK: CROSS SECTION OF SEMI-INFINITE BODY 51 


II. BODIES OF REVOLUTION NOSE-ON 
As the most straightforward example of scattering from a semi-infinite body, 
these ideas will now be applied to backscattering from a body of revolution 
with the dipole source located on the axis of symmetry (taken as the z axis). In 
cylindrical co-ordinates, 2, p, %, 


(11) dS = p{1+(dp/dz)?]* do dz, 

(12) r = 2 (R+z)+Xp cos ¢+9p sin ¢, 
(13) r= (R+2)?+)’, 

(14) ? XX = 1 [9(R+z) —Zp sin 9], 


2x 


(15) do(?Xk) = 2nhr7' (R+z). 
0 


The slope of the tangent to the body is (dp/dz). The inclination with respect 
to the z axis of the ray to a point on the body is tan[p/(R+2z)]. The angle y 
between the incident ray and the tangent to the body is given by 


(16) tan y = tan {tan=! (dp/dz) —tan [p/(R+2)]} 
= [(R+2) (dp/dz) — p]/[R+2+p(de/dz)], 

from which it follows that 
(17) #-% = —sin y = —[(R+2)(dp/dz) —p)/{[(R+2)*+e"[1 + (de/dz)*]} * 

= —[(R+2)(dp/dz) —p]/r[1+(dp/dz)*]}. 
With these values, equation 10 becomes 
(18) HY = ~ik Hy f “de o(R+2)[(R+s)(dp/ds)—p] te ™, 
From equation 13, 
(19) rdr = [R+2+ (dp/dz)| dz 
so that 
(20) H*® = —ik Hos Jira e**" 5(R+2)[(R+2) (dp/dz) —p] 

x [R+2+p(dp/dz)}*. 


Repeated integration by parts of equation 20 yields 


(21) H* = —(Ho/2) & 25 > (ik) g”(R) 

where 

(22) g(r) = r-*(R+2)[(R+2)p(dp/d2) —p'[R+2+p(dp/dz)}“ 
and 


(23) g(r) = (d/dr)*g(r). 
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The terms evaluated at infinity must, of course, vanish. 
For a clearer view of the relative order of terms, it is helpful to change to 
the dimensionless quantities 


(24) t=r/R, u=2/R, v = p/R, K = 2kR, 
(25) G(t) = Rg(r) = 41+) [((L+u)v(dv/du) —v*][1 +u+v(dv/du]-, 
(26) G™(t) = (d/dt)"G(t) = R"(d/dr)"G(t). 
Thus 
(27) H® = —(/2R) * 5 GK)" G™ (1). 
n=0 


Provided the G(1) do not differ inordinately from each other, successive 
terms in the series are smaller by a factor of order (RR)~!. 

Inasmuch as a surface of revolution is most naturally described by specifying 
p as a function of z, it may be more convenient to use equation 19 to write 


the derivatives as 


(28) d/dt = (dz/dr)(d/du) = t{1+u+v(dv/du)|-'(d/du) 
so that equation 27 becomes 
(29) H* = —(1)/2R) e~*$ SS GK)" {t{1+-u+v(dv/du)}(d/du)}"G]yno, 
n=0 
where 
t = [(1+u)?+v*]} 
and the surface is defined by 
v = v(u). 


(4) Paraboloid 
The equation of the paraboloid is 


(32) p? = 2cz 


or, in terms of the dimensionless variables, 


(33) v? = 2cu/R. 

Hence, 

(34) G(1) = (c/R) [1+(¢/R)}" 

with all the higher-order terms finite. To lowest order, then, 
(35) H* = —9 (Hoc/2R?*) (1+ (c/R)]-! e-2#* 


and according to equation 3 
(36) o = mec’. 


This is Kerr’s answer (Kerr 1951) obtained directly, without resort to trunca- 
tion of the paraboloid and the associated Fresnel zone argument. 
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(B) Cone 

Use of the physical optics approximation for the infinite cone cannot be 
formally justified, since the tip does not satisfy the condition that the radius 
of curvature be greater than the wavelength. However, physical optics has 
often given good answers beyond its range of proved applicability, and it has 
accordingly been used extensively beyond this range, in particular for the 
cone. The equation of the curve is 


(37) p=ztana 

where a is the cone half-angle, or 

(38) v = utana. 

In this case, 

(39) G(1) =0 

(40) G'(1) = tan? a 

with all the higher-order terms finite. Hence, to lowest order, 
(41) H® = i (Ap tan? a/4kR?) 2? ** 

and 

(42) o = (r/4k?) tanta = (A2/167) tan‘ a. 


This answer was obtained by Kerr’s approach in Siegel and Alperin (1952). 


III. BROADSIDE SCATTERING 


The attractiveness of the concept of a cross section stems from the assump- 
tion that, if the source and observer are sufficiently distant from the scattering 
body, the scattered intensity will be independent of the range except for the 
geometrical factor R-*. There are configurations for which this is not true 
(i.e. for which the scattered intensity at large distances does not fall off as 
R~). For instance, certain bodies have an R~* dependence for broadside 
incidence; heuristically, this is identified with a dependence of the scattering 
on the size of the first Fresnel zone. For these cases, the formulation of the 
cross section as the infinite range limit does not apply; in fact, equation 3 
would give an infinite answer. Nonetheless, the term ‘‘cross section”’ is often 
used to denote the quantity that would have to be used in the radar range 
equation to describe the scattering at large distances, even though the ‘‘cross 
section’’ so defined differs from the more common variety in containing a 
dependence on the range (no matter how great the latter be). The extension of 
equation 3 to cover this contingency involves replacing the limit R > © by the 
more intricate statement that R be taken sufficiently large (albeit finite) that 
in an expansion in powers of R only the highest-order term need be retained. 
The normal situation then corresponds to this highest-order term being 
independent of R (in which case the two definitions are equivalent), while in 
the present instance the leading term will be linear in R. 
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(1) Circular Cylinder 
As an example, the “cross section”’ of an infinite circular cylinder for broad- 


side backscattering is computed from equation 10. The axis of the cylinder is 
taken as z axis, the normal from the cylinder axis to the source as y axis, and 
the polarization direction as x axis. Then 


(43) dS = pd¢ dz, 

(44) r = —3(R+p)+22+ 6p, 

where 

(45) Pp = Xcos¢+9 sind = 4%, 

so that 

(46) r = Xp cos ¢—J[R+p(1-sin $)]+2z, 
(47) r = 2+R?+20(R+ p)(1— sing), 
(48) XX = r{2[R+p(1— sin $)]+5z}, 
(49) 7-2 = —r—[Rsin ¢—p(1— sin ¢)]. 


At long ranges (R>> p), ¢ in the illuminated region extends from 0 to x. The z 
component of the integrand is symmetric in 2, the y component antisymmetric. 


Hence, 
(50) H*® = —(ikp/r)Hoz J d¢[R’ sin ¢—p(R+p)(1—sin ¢)"] f ee: 
0 0 


Subject to R >», r can be expanded: 


(51) or = (2° +R*)$+p(R+p)(2?+R?)- 1(1 —sin g) —2-1p°(R+p)?(2+R2) 32 
X(1—sing)?+.... 


The third term can be neglected in the phase if 
(52) (kp/R) K2r, 
which is the usual Fresnel far-zone requirement 
(53) R> (p?/A). 


Only the first term of equation 51 need be retained for the r~ factor. A couple 
of p's can also be consistently dropped in other factors, reducing equation 50 to 


(54) H* = —(ikR’p/x) Hod if de(2?-+R?)? exp{ —2ik[(z2°+R”)! 
+Ro('+R)4}} xf dafsin ¢— (o/R)(1-sin 9)" 


x exp[2ikRp(z°+R’) sin @]. 
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The first term in the ¢ integration is 


(55) f dg sin ¢ exp[2ikR p(s°-+R’)? sin ¢] 
0 


= j dosing cos[2kRp(s° +R’)? sing]+i f d¢ sing sin[2kRp(z’° +R’) ' sing] 
0 0 


wE,[2kRp(2° +R’) 4) +ixJ[2kRp(2° +R’) 4] 

where £; and J; are the Weber and Bessel functions of order one, respectively 
(Watson 1952). The second term yields the same functions of orders zero, one, 
and two. The contribution of the second term is thus smailer than that of the 
first by a factor of the order of magnitude of (p/R), and it will henceforth be 
neglected. Then 


(56) H*® = —ikR’pH3 f dz(2°+R*)~*{ Ey[2kRp(2’+R’) 4] 
0 


+iJ\[2kRp(2’ +R*)*]} Xexp { —2ik[ (2° +R’)! + Ro(2? +R’) ]}. 
Equation 56 can be evaluated by stationary phase. To highest order in R 
outside the phase factor, the result is 


(57) H® = —i2-w*R-*pHo3[Ei(2kp) +iJi(2kp)] exp{ —i[2k(R+p)+(x/4)]}, 
which leads to 
(58) ¢= mkRp*|E\(2kp) +151 (2p) |?. 


The polarization has been assumed above to be transverse to the cylinder 
axis. For longitudinal polarization, there would correspond to equation 48 


(59) 7X2 = —r{&[R+p(1—sin ¢)]+p cos ¢}. 


To leading terms in R, this would lead only to the change 2 — —£ in equation 
57 and leave equation 58 unchanged. Thus, the far-zone intensity of the field 
scattered by an infinite circular cylinder shows no polarization dependence in 
the physical optics approximation. 

In the course of this calculation, explicit use has been made of several 
requirements that the range exceed other characteristic dimensions. An 
additional condition that occurs in the derivation of the physical optics approxi- 
mation (although there has been much pragmatic success in using the method 
despite its violation) is that the wavelength be smaller than the radii of curva- 
ture of the surface, hence that 2kp > 1. If this last inequality be also satisfied, 
the Weber and Bessel functions can be replaced by their asymptotic values 
(Watson 1952): 


(60) E,(2kp) = — ¥i(2kp)+0[(2kp)-*] = —(rkp)—} sin [2kp—(3x/4)] 
+0[(2kp)-*”], 


(61) J1(2kp) = (rkp)—* cos[2kp — (3%/4)]+0[(2kp)-* I, 
(62) E,(2kp)+iJi(2kp) = i(wkp)—* exp{i[2kp — (3x/4)]}. 
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Equation 57 now reduces to 
(63) H® = —2-'R-3p 4g e-2#F 
and equation 58 to 

(64) o = tRp. 


It is interesting to note that the cylinder ‘‘cross section’’ obtained above 
(equations 58 and 64) is a factor of two smaller than would be obtained with 
the usual assumption of an incident plane wave. Thus, tampering with the limit 
processes can alter the answer. 

The limiting value of the “cross section”’ given by equation 64 is easy to 
grasp intuitively. The geometric optics cross section of a conducting body for 
specular backscattering is usually given by 


(60) iC - wR Rs 


where R, and k: are the principal radii of curvature at the specular point, and 
both are large compared with the wavelength. In the case of the cylinder, one 
of these radii of curvature is infinite (the longitudinal one) and the expression 
breaks down. A bit of care in tracing through the derivation of the geometric 
optics formula elicits the fact that it is the relative curvature of the body and 
wavefront which is pertinent. Ordinarily, the curvature of the wavefront at a 
large distance from the radar is negligible compared with that of the body, and 
equation 65 results. In the present case, the opposite is true along one axis and 
then the radius of curvature of the wavefront enters instead of the radius of 
curvature of the body, leading to equation 64. On the other hand, the assump- 
tion of plane-wave incidence takes this relative curvature into account going 
one way and not the other, thus underestimating the effective curvature 
(overestimating the effective radius of curvature). 
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MEASUREMENT OF RADIATIVE CAPTURE RESONANCE 
INTEGRALS IN A THERMAL REACTOR SPECTRUM, 
AND THE THERMAL CROSS SECTION OF Pu-240' 


W. H. WALKER, C. H. WEstTcotTt, AND T. K. ALEXANDER 


ABSTRACT 


An apparatus is described which detects y-rays emitted by a thin target 
placed in a well-defined neutron beam. It has been used to determine the Cd 
ratios of Au and Pu-240, from which the ratio of the resonance integral to the 
2200 m/s cross section for radiative neutron capture in Pu-240 has been deduced, 
using Au asa reference standard. Using this ratio and previously measured values 
of the resonance integral of Pu-240 and its effective cross section in two positions 
in the NRX reactor, three separate estimates of the 2200 m/s cross section of 
Pu-240 have been made. The mean value is 270+17 barns. 

In an auxiliary experiment to indicate the shape of the epithermal spectrum of 
the neutron beam, the activation Cd ratios of Mn and In were compared with 
that of Au. These results, combined with the known 2200 m/s capture cross 
sections of these nuclides, yield new values of the radiative capture resonance 
integrals for both Mn and In. 


1. INTRODUCTION 


When this experiment was started, the existence of a strong resonance in 
Pu-240 at about 1 ev was known, but accurate values of resonance parameters 
and thermal cross section were not available. Initially therefore, this investiga- 
tion was primarily directed at measuring the cadmium ratio for the capture of 
reactor neutrons in Pu-240, using an apparatus which was also intended to be 
used with other nuclides such as fission products. 

Although the results for Pu-240 have been given in an earlier paper (Westcott 
et al. 1958) a full description seems worth while, especially as the techniques 
used for observing neutron capture have not been widely used except in con- 
junction with time-of-flight work where the expected background counting 
rates are much lower than in our conditions. In this connection it may be 
interesting to note that the three Nal (Tl) crystals used as y-ray detectors 
were originally arranged close to the Pu target foil to provide maximum 
efficiency in case coincidence measurements were required. Heavy shielding 
was necessary to protect these counters from the radiation near the reactor and 
a considerable effort was required to develop a suitable collimator for the 
neutron beam. 

In the course of calibrations, incidental determinations were required of 
resonance integrals for two nuclei commonly used as standards. As these did 
not confirm the accepted values, details of the work and its conclusions are 
presented in addition to results on Pu-240. 

The relation between cadmium ratios and the ratio of resonance integrals 
to thermal cross sections is also more fully expoynded in the present paper 
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than has been done in earlier publications (Westcott 1956; Westcott, Walker, 
and Alexander 1958). It may be noted that the neutron spectrum in the 
moderator at the center of a large thermal reactor, such as the NRX reactor 
at Chalk River, is well suited for the application of this theory. 

Several recently published time-of-flight measurements (Fluharty 1957; 
Cété et al. 1958; Leonard et al. 1958) of the resonance parameters of Pu-240 
are in agreement, and predict a 2200 m/s cross section which has been con- 
firmed by a measurement with a crystal spectrometer (Pattenden 1958). By 
combining our value for the ratio of the resonance integral to the 2200 m/s 
cross section with determinations of the effective Pu-240 capture cross sections 
for positions in NRX, we obtain a value of 270 barns for the 2200 m/s cross 
section, in agreement with the spectrometer results. However, it should be 
noted that other recently published results on the resonance parameters 
(Evelstaff et al. 1958) and thermal cross section and resonance integral (Rose 
et al. 1958) differ appreciably from this value. 

2. BEAM COLLIMATION, SHIELDING, AND y-RAY DETECTION 


The essential components of the experimental arrangement used to detect 
capture y-rays are an intense, well-collimated beam of neutrons, a thin film of 
target material in the beam, and efficient y-ray detectors placed around the 
target position. These are shown in Fig. 1. The neutron beam is obtained from 
a graphite scatterer in an “instrument” hole which is tangential to the wall 
of the NRX calandria, as shown in the inset sketch. The 11-ft collimator unit 
extends outwards from a point 73 ft inside the pile shielding to the detector 
cavity inside the main experimental shield. The collimating hole is placed 
above the center line of the instrument hole so that a rotating steel cylinder 
located 2 ft inside the pile face can be used as a gate to shut off the beam. The 
main shielding is contained by walls of welded }-in. steel plate, mounted on 
two 4-wheeled carriages running on rails as shown in Fig. 1. The sections 
adjacent to the counter cavity and beam path are filled with lead, the remain- 
der with boron-loaded gelatin. The forward section, which carries the detectors, 
preamplifiers, and pump, normally remains in position against the pile face, 
and the rear carriage is moved back when it is necessary to gain access to 
counters or target. 

In its initial form the collimating hole was 1 in. in diameter over its entire 
length. Alternation of short sections of graphite, BsC powder in aluminum 
cans, and Pb served to slow down and capture neutrons, and to reduce y-ray 
intensities outside the beam. The background due to capture of scattered 
neutrons in the Nal (TI) detectors and the materials surrounding them 
proved to be greater than expected. This background was reduced to a low 
level by evacuating the beam path near the counters and using a collimator 
insert. The evacuated section extends from near the filter position to the ByC 
beam stopper while the insert starts near the fission chamber at the left-hand 
margin of Fig. 1 and ends at the end of the original collimator. The beam is 
reduced by the first (tapered) portion of the insert to a diameter of 0.5 in. at a 
point 2 ft in front of the target position, while from this point on, the insert 
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has an internal diameter of ? in., large enough to clear the neutron beam, but 
small enough to shield the scintillators from neutrons penetrating or scattering 
from the constricted }-in. section. Thus no neutron scattered only once from 
the walls of the collimator can enter the detectors. 

The target material is deposited on a .0003-in. thick Al foil mounted on a 
l-in. diameter Al ring. Both the ring and the holder used to place it in the 
target position are clear of the neutron beam. The aluminum-canned ByC 
shown in Fig. 2 surrounding the target tube helps reduce the background in 
the detectors due to scattered thermal neutrons. 

Capture y-rays from the target are detected in three Nal (Tl) scintillation 
counters arranged about the target position as shown in Fig. 2. Two of the 
counters are 4-in. diameter by 2-in. thick Nal (TI) crystals mounted with 
lucite light pipes on 2-in. RCA 5819 photomultiplier tubes. The section 
through one of these and its bismuth shielding is shown. 

Photomultipliers with a small photocathode area were chosen in preference 
to 5-in. diameter types because the latter had a poor gain-stability character- 
istic under changing counting rates. After counting at high rates, the gain 
returned to its low counting rate value exponentially with a half-life of about 
28 minutes, confusingly similar to the half-life of the [-128 activation induced 
by neutron capture in the iodine of the crystal. The use of a light pipe instead 
of direct coupling entailed some loss in resolution. 

The third counter consists of two independently sealed 2-in. diameter by 
2-in. long Nal (Tl) crystals placed end to end and viewéd separately by two 
RCA 5819 photomultiplier tubes. The collector outputs are added at the 
input to a cathode follower preamplifier. The gains were equalized by adjust- 
ment of the voltage across the dynode resistor chains until the counters gave 
optimum photopeak resolution when detecting Na-22 y-rays. 

To reduce the effect of stray magnetic fields on the gain of the photo- 
multipliers, they were surrounded by mu-metal shields. 


3. ELECTRONIC APPARATUS 

To record the y-rays emitted from the target, a counting system was re- 
quired in which the pulses from the scintillation counters, after pulse height 
discrimination, can operate either a 0.1-usec coincidence circuit or scalers. The 
design of the circuits ensured that the rise time of the pulses at the output of 
a linear amplifier and the time jitter of triggering of the discriminators were 
less than 0.1 usec. Hence it was only necessary to shape the output pulses of 
the discriminator, by clipping, into pulses with 0.1l-ysec half-width. These 
shaped pulses then operated a diode coincidence gate. 

Figure 3 shows a block diagram of the counting system and demonstrates 
one typical mode of operation. The current pulses from the photomultiplier 
tubes are integrated (time constant > 0.25 usec, the time constant of Nal (T1)) 
at the input of the cathode follower and shaped to pulses with approximately 
0.1-usec rise time and 0.5-usec decay time constant at the output of the 
cathode follower. The pulses are amplified with a maximum gain of 20 and 
passed directly to the input of a discriminator which triggers from negative 
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Fic. 3. Block diagram of the counting system showing a typical mode of operation. 


pulses in the range 0.1 to 1 volt. The discriminator produces a positive, 10-volt, 
fast-rising pulse when triggered. The dead time of the discriminator is 0.7 usec. 

The output pulses of the discriminator operate a scaler or, if coincidences 
are to be recorded, these pulses are clipped with a 0.05-usec shorted delay line 
to give pulses 0.1 usec wide. This width then determines the resolving time of 
the coincidence circuits. 

The coincidence unit, shown in Fig. 4, has four inputs A, B, C, and D. Four 
diode gates impulsed through ferrite-cored pulse transformers make up the 
coincidence gates to provide outputs A+B (—D), A+C (—D), B+C (—D), 
and A+B+C (—D), where + indicates coincidence and (—D) indicates 
optional anticoincidence. If all diodes in a given gate are cut off, the current 
flowing then charges the stray capacity on the input of the shaper and 
triggers it to produce a 15-volt output pulse. 

The resolving time was checked experimentally by measuring the chance 
coincidence rates when independent Co-60 y-ray sources irradiated each 
crystal separately. The measured value was 0.12 usec +5%. 
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Fic. 4. Circuit diagram of coincidence circuits. 


4, COUNTING CAPTURE Y-RAYS 


The capture rate in the target foil may be measured by either single or 
multiple counting of capture y-rays. The range of y-ray energies used in single 
counting is limited by the growth of I-128 in the scintillator due to the capture 
of stray thermal and resonance neutrons. Since this nuclide decays with a 
25-minute half-life, emitting a 2-Mev -particle in 93% of the disintegrations, 
and a 1.6-Mev £-particle in coincidence with a 0.4-Mev y-ray in the remainder, 
all counts below 2 Mev are usually biased out in single counting in order to 
avoid a time variation of the background. 

With multiple-coincidence counting, the rates are much lower, and the 
signal-to-background ratio higher. In double-coincidence counting, the low 
y-ray energy limit is again set by the I-128 decay; in this case, however, 
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pulses from all y-rays exceeding 0.4 Mev are acceptable. In triple-coincidence 
counting there is no lower limit to acceptable y-ray energies. 

The beam is monitored by the small U-235 fission chamber shown in Fig. 1, 
which counts neutrons scattered by a lucite ring at the inner end of the 
collimator insert. 


5. THEORETICAL INTERPRETATION OF CADMIUM RATIOS 


The treatment developed here is based on the convention for cross section 
and neutron flux first proposed by Westcott (1956) and adopted since by many 
laboratories (Pattenden 1958; Rose et al. 1958; Eastwood ef al. 1958); a 
similar scheme is used by Halperin and Stoughton (1958). It is applicable to 
large well-moderated reactors, and assumes that the neutron flux is composed 
of Maxwellian (or thermal) and epithermal components. The Maxwellian 
component is defined by a temperature 7, and the epithermal component is 
proportional to dE/E with a cutoff at wk7T, where uw = 5 for heavy-water- 
moderated reactors. The appendix contains a list of the symbols used in this 
paper with their definitions. 

The development of the definitions leading to the equation 


(1) & = oo(g+rs) 


can be found in an earlier paper (Westcott et a/. 1958). In this reference it is 
shown that, for a detector with cross section proportional to 1/v below about 
0.5 ev, the quantity s is proportional to W(7T/T>). In the treatment which 
follows the invariant quantity sx/(To/T) will be denoted by so. 

In the present paper we are concerned with the interpretation of cadmium 
ratios. For an infinitely thin detector in which the cross section is nearly equal 
to gm(vo/v) between wkT and Ega, the equation of Westcott ei a/. (1958) written 
in terms of sp and 7x/(7/T») becomes 


(2) ioe (, / ae 

ee To g/l = 
The parameter K is a function of the thickness of Cd, tabulated by Westcott 
et al. (1958) for both a normally incident neutron beam and for isotropic 
incidence on one side of a Cd foil, as required for a ‘“‘sandwich”’ detector inside 
a reactor. If a sharp cutoff at an energy Ecaq is assumed for the neutrons 
transmitted through the Cd, then 


oo Le! Seon 
4 Eo 
For application to the present case, equation (2) needs further refinement. 
The term {rx/(7/T»>)}-! mainly represents the contribution of thermal 
neutrons; if the Cd is thin enough a fraction of these can pass through the Cd, 
so that it is necessary to add to the denominator the term h{r/(T/T»)}~, 


where / is the transmitted fraction of the neutron density for the initially 
Maxwellian component. 
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The resonance capture contribution, so/g, is defined by the equation 

(3) Bing 2 fe ee 
g V © J uxr \2o0 E/JE 

ea | {t-4/ 2) e.2.£ (Scns a) dE 

a VT ukT \£%o E E V4 Ecoa \Z7o E Ye 
The second integral is, for our cases, essentially invariant with respect to T 
and the first term is a small correction, W say, arising in those cases where o 
is not proportional to 1/v. Values of W have been estimated for Au, In, Mn, 
and Pu-240, assuming that the cross section between wkT and Egg is all due 
to the first resonance; they are Way = 0.027, Wy, = 0.245, Wun = 0, and 
W pu-240 = 0.35. 

The last integral of equation (3) is reduced by resonance self-shielding in 
the target foils and by the attenuation of resonance neutrons in the Cd 
filter. The first of these factors can be calculated using the method of Roe 
(1945) for predominantly capturing resonances if their parameters are known. 
For the target thickness of Au, In, and Pu-240 used in these experiments self- 
shielding is important for only one resonance, while in the case of Mn, where 
the resonance is predominantly scattering and the method of Roe does not 
apply, the targets are sufficiently thin that self-shielding is negligible. 

The attenuation of the resonance neutrons in the Cd filter is a function of 
cadmium thickness (NV atoms/cm?), neutron energy, and angle of incidence. 
For the beam irradiations this reduction is given by the factor exp{ —ctotai 
X N(.030 in.)}, while for those in the pneumatic carrier it is given by the ex- 
ponential integral function Ee{ocapture N(.015 in.)}. In the latter case it is 
assumed that the foil sandwich approximates an infinite slab and the effect 
of potential scattering (~5 barns) is negligible. The slopes of these functions 
change sufficiently slowly with neutron energy that the value at the resonance 
peak energy (E,) may be used as the filter reduction factor. This assumption 
is valid for narrow resonances in the region of 1 ev such as Pu-240, but would 
not apply to those for which the width T was appreciably larger. 

It is therefore necessary to multiply the resonance capture contribution in 
the numerator (so/g) by Gi, and that in the denominator {(so/g)—W} by 
FG,, where ‘‘F’’ is the reduction factor due to the cadmium filter and ‘‘G,”’ is 
Roe’s self-shielding factor for the dominant absorption resonance of the target. 
Equation (2) then becomes 


0 nal FY saphrazes way Y 


using the approximation W = FG,W. 
Since all materials are to be compared with Au, a more useful equation is 


(5) Es Be, nf eis ONO Eat 
(FRea—1) au { (Gis0) /g}e+(1/Fe){(1/K)—Weth(rx/T/To) } ’ 
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where 





: 5 _(rVT/To)"(1/Fe){ A/—K) — Wet h(r/T/T 0) '} 
on * (PVT /T) = (1/Faa) { 1/K) — War thr T/T) 4} 
=1+W,rvV/T/T», 
since hK1, Wau K1, and K >rV/(T/T»). In these experiments ro/(T/To) 
is equal to 0.025 to the degree of accuracy required, so that Zyn = 1.0, 
Live = 1.006, and Z pu-240 = 1.009. 
6. CADMIUM RATIO MEASUREMENTS WITH GOLD 


In the measurements reported in this paper gold has been used as a refer- 
ence standard. The gold target used is a deposit of 1 mg/cm? evaporated on a 
0.3-mil Al foil. Table I lists recent measured values for the 2200 m/s capture 


TABLE I 
Gold cross sections 





Sd 
~ 


4.9-ev Other 
Reference oo o;(th) resonance resonances 





Carter, Palevsky, Myers, and 

Hughes (1953) 98.7+0.6 — 
Egelstaff (1954) 98.4+0.9 9+2 — — 
Allen, Stephenson, Stanford, 

and Bernstein (1954) ~- 8.7+40.9 — — 
Gould, Taylor, Rustad, 

Melkonian, and Havens (1957) 98.8+0.3 = -— 
Seidl et al. (1954) — — 
Wood, Landon, and Sailor (1955) — — 1470+37 
Argonne (1955) — — — 55+12 


42+15 


cross section, the thermal scattering cross section o,(th), and the capture 
resonance integral. The capture resonance integrals 2’ have been calculated 
from published resonance parameters since it can be shown that, for a single 
Breit—Wigner resonance, 


(7) 2’ = (x/2)o,T,/E, 


provided E,>>ukT, so that o is proportional to 1/v in the region of ukT. 

Values adopted in this paper are oo = 98.6 barns, o, = 9 barns, and 
2’ -= 1525) barnsso that:(So)aq = 1045. 

The validity of using capture y-ray counting to measure the Cd ratio of Au 
has been established by comparison with activation results. Agreement to 
1% and better was obtained between all three counters, singly and in co- 
incidence, taken over a range of biases. For single counts the range investigated 
was between 1.5 and 3.0 Mev. In addition Cd ratios obtained by activations 
in the beam agreed within the experimental error of 2% with values obtained 
in the same experiment by single counting of capture y-rays. No appreciable 
contaminant was detected in the activation experiment. 
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7. SELF-FILTERING WITH GOLD 


A further test of the reliability of the detection method used in these experi- 
ments is a measurement of the epi-Cd neutron transmission through Au filters, 
using a thin Au detector. If n; = o,f; and n, = o,f,, where ¢; and ¢, are the 
thicknesses of the filter and target foil respectively, in atoms/cm’, and 4g, is 
the resonance peak cross section in square centimeters, then it can be shown 
(Roe 1954, p. 32) that the transmission for this resonance is given by 


(8) T> (netne)Gi(nr tne) —m1Gilar) 
a mGi (ne) 


ointot oe 
— (Gi(m)) an ymGr(aet am) § 


Here ‘‘G,y’’ is the self-shielding factor referred to in Section 5. The approxima- 
tion in equation (8) introduces a negligible error in the cases of gold and 
Pu-240. 

If it is assumed that the signals received by the counters from the various 
resonances are in proportion to their resonance capture integrals, the over-all 
calculated transmission is the sum of the contributions from each resonance, i.e. 
(9) Teac = RPE ee: 


— 


where 2/ is the resonance integral contribution of the m™ resonance and 
7. sf of 
vo’ = 2,4+2,+.... 


In Table I] the results of measurements of epi-Cd neutron transmissions 
with their statistical errors are compared with calculated values. The error 


TABLE II 


Transmission of epi-Cd neutrons 


Au filter 





thickness, 
mg/cm? T calc (T meas — 030) /0.97 
4.67 0.816+0.010 0.810+0.005 0.803 
9.48 0.673+0.010 0.673+0.003 0.662 
14.15 0.573 +0.010 0.578+0.002 0.564 
22.93 0.434+0.010 0.459+0.002 0.442 
32.41 0.359+0.010 0.380+0.001 0.360 
1490 0.050 +0 .005 0.079+0.004 0.050 





assigned the latter is estimated and is due primarily to errors in reading the 
curves. The last measurement (very thick filter) was done separately, in 
conjunction with the temperature measurement reported in the next section. 

The two sets of transmissions show a monotonic divergence which is well 
outside the limits of error for the thicker filters. This can be accounted for by 
assuming a constant signal in the counters from the target, equal to 3% of the 
epi-Cd signal with no Au filter, as indicated in the last column. Such a signal 
could arise from resonant capture in unidentified high-energy resonances, 
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from capture in the detectors of potential-scattered neutrons, or from capture 
in a contaminant in the target. The first possibility is improbable since it would 
require a contribution to 2’ equal to that from all the identified resonances 
other than the one at 4.9 ev. A scattering signal of 3% is not consistent with 
the agreement to better than 2% between in-beam and activation comparisons 
of the Cd ratio of Au, as described in Section 6. Since the presence of a con- 
taminant cannot be similarly ruled out nor its effect on the Cd ratio calculated, 
it is only possible to say that the Cd ratio of Au determined by capture y-ray 
counting is correct to 2% or better. 


8. THE “TEMPERATURE” OF THE BEAM NEUTRONS 

To determine the temperature of the Maxwellian component of the beam 
neutrons a measurement was made of the transmission through a 1.49 g/cm? 
gold filter, of neutrons having energies less than the Cd cutoff, using the 
standard gold foil as detector, with single counting. The measured transmission 
was 0.645. 

Values for transmission using 1/v filters and a thin 1/v detector have been 
calculated as a function of filter thickness and neutron temperature. Correc- 
tions were made for scattering in gold, for the departure of the gold cross 
section from 1/yv (g = 1.01 at 130°C) (Westcott 1958) and for 1/v capture 
due to nitrogen (a9 = 1.88 barns) (Hughes and Harvey 1955) in the flight path. 

The temperature corresponding to the corrected transmission was 130+10° C 
where the assigned error in the measurement is due chiefly to the uncertainty 
in the scattering cross section of Au. 

This result is in good agreement with thermocouple measurements in the 
graphite outside the calandria, and with other measurements of neutron 
temperature in the graphite at two self-serve positions (Bigham 1958). 

The effect of the nitrogen in the 5-meter flight path between the graphite 
scatterer and the filter position is made up of two parts, a scattering effect 
causing a loss of about 20% for neutrons of all energies, and the capture of 
neutrons in the (”,a) reaction. At a neutron temperature of 130° C the latter 
effect is estimated to be a signal reduction of 33% for Maxwellian neutrons, 
and of 1% for epi-Cd neutrons. No correction is required if the results are 
treated using equation (5). 


9. THE EPITHERMAL SPECTRUM OF THE BEAM NEUTRONS 


In the derivation of equation (5) it is assumed that the epithermal spectrum 
is proportional to dE/E, the distribution function predicted by neutron 
diffusion theory for epithermal neutrons in the moderator near the center of a 
well-moderated reactor with a large resonance escape probability such as 
NRX. The irradiation station of the NRX pneumatic carrier is located in an 
empty lattice position in the median plane of the reactor near the central axis. 
Unfortunately, the absence of fuel in this position depresses the neutron 
spectrum at higher energies below the 1/E distribution which occurs elsewhere 
in the moderator. A calculation of the flux at this location in the region of the 
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lowest Mn resonance (337 ev) and near the gold resonance (5 ev), assuming 
Fermi age distributions about line sources of 3-Mev neutrons, indicates that 
the former is 14% smaller than the 1/E extrapolation of the latter. Although 
the higher Mn resonance will be affected to a greater extent, corrections for 
the various simplifying approximations made and for the resonance capture in 
U-238 in the fuel will tend to reduce this discrepancy, so that 14% may be 
considered an upper limit. The deviation from 1/E predicted between 1 ev 
and 5 ev was less than 1%. 

The epithermal spectrum of the beam neutrons was compared with that in 
the pneumatic carrier position in the moderator by measuring the Cd ratios 
of Au, In, and Mn. All measurements were made by activating the target 
materials, with the exception of those for Au in the beam. In this case capture 
y-rays were counted before and after the beam activations of Mn and In. 
Pertinent information on the activation targets is given in Table III. 





TABLE III 
Activation targets 
Target Cd Irradiation period, min 

Irradiation Target thickness, thickness, ——————— —— Activity 
position element pug /cm? inches No Cd Cd measured 
Pneumatic carrier Mn 240 .015 1.00 10.0 Mn-56 
(2.6 hr) 

Pneumatic carrier In 60 O15 1.00 3.00 In-116 
(54 min) 

Pneumatic carrier Au 80 .015 In same Au-198 
capsules (2.7 days) 

as Mnand In 

Beam Mn 3500* .030 5.0 60.0 Mn-56 
In 840* .030 10.0 29.0 In-116 


*Mounted on 1-in. diameter Al ring and irradiated in regular target position (see Fig. 2). 


The activation counting was done in the same geometry in all cases. The 
detector was a 2-in. X2-in. Nal (Tl) crystal coupled to an RCA 6655 phototube, 
and shielded by 2 inches of Pb. To give the largest signal-to-background ratio, 
only the region about the photopeak was counted, using a single-channel 
analyzer. Corrections for the activity of the aluminum backings were measured 
by counting foils which were irradiated with the targets. 

The relative masses of the target materials used in the pneumatic carrier 
irradiations were measured by activation counting following subsequent 
irradiations in which the pairs of foils were placed face to face, separated only 
by a thin foil of Al. The relative masses of the beam targets were measured by 
comparing the capture y-rays from each pair. The values obtained for Rea in 
these measurements are listed in Table IV. 

From equation (5) it can be seen that {(FRca — 1)2}/{(FRea — l)au} 
will be sensibly independent of the epithermal flux density if # = 0 and so 
remains constant, i.e. if the shape of the epithermal flux spectrum remains 
constant. Auxiliary conditions are that G,, g(T °C), 1/K, and F (on the left- 
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TABLE IV 
Measured values of Rca. Corrected values of {(FRcea-1)z}/{(FRea-1) au} 


ete = a 


Shielding 

factors {(FRca-1)2} / 
Irradiation Target Rea ——————— Normalization {(FRca-1)au} 
position element (measd.) F r factor (normalized) 


1.034 15.80+0.5 


OQ 


Beam Mn 39.4+0.7 0.98 
Au 3.53+0.07 0.98 


SS 


com 
oo 


1.034 15.70+0.5 


Beam Mn 37.3+0.55 0.98 
Au 3.41+0.07 0.98 


SS 


oco- 
oo 


Pneumatic carrier Mn 31.1+0.4 1.00 1.060 15.64+40.2 


Au 3.04+0.01 1.00 


Beam In 3.20+0.05 0.96 
3.32+0.07 0.98 


ss 


0.895+0.035 


oo 
Oo 
“ID 
0 & 


Pneumatic carrier 2.955+0.018 0.95 0.886+0.015 


3.04+0.01 1.00 


ss 


hand side of the equation) remain the same. In the present experiment these 
factors vary and ‘‘h’’ is not always zero, but since their values can be calculated 
it is possible to normalize the measured ratios to a standard condition. The 
shielding factors F and G, are listed in Table IV; the thermal neutron penetra- 
tion factor, ‘‘h’’, is 0.0025 for the thickness of Cd used in the pneumatic carrier 
irradiation and zero for the beam; 1/K = 0.522 for the thickness of Cd and 
direction of neutron incidence in the beam, 0.502 in the pneumatic carrier; 
for the beam, with T = 130°C, gau = 1.012, gun = 1.000, and gin = 1.041; 
for the pneumatic carrier, with T = 38°C, gay = 1.006, gman = 1.000, and 
Zin = 1.022 (Westcott 1958). The standard conditions chosen are: h = 0; 
F = G, = 1.00 for every case; 1/K = 0.522; and T = 20° C. The normaliza- 
tion factors are listed in Table IV, as are the normalized values for {(FRca 
—1)2}/{(FRea — lau}. 

The similarities of the values in Table IV for both Mn and In for the two 
locations indicate that the epithermal spectra of the beam and in the vacant 
lattice position have substantially the same form within the energy range 
1 ev to 1000 ev. 


10. THE ACTIVATION RESONANCE INTEGRALS OF Mn-55 AND In-115 


The values of so for Mn-55 and In-115 can be obtained directly from the 
weighted mean values of {(FRcea—1)2}/{(FRea—l)au} of Table IV using 
equation (5) and the values of g, F, G1, and 1/K given above for the standard 
conditions. The resonance integrals, 2’, can also be obtained since the thermal 
capture cross sections of both Mn-55 and In-115 are known (Hughes and 
Harvey 1955; Rose et al. 1958). Values of so, oo, and 2’ are given in Table V. 
The value of 7.2+0.2 obtained for Mn on this basis has been adjusted to the 
value given in the table to allow for the depletion of the flux at 337 ev relative 
to that at 5 ev discussed in the previous section, the error quoted having been 
increased to allow for uncertainties in this correction. 
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TABLE V 
Values of so and 3’ for Mn, In 
Target g, 20°C So oo, barns >’, barns 
Au 1.005 17.45 
Mn 1.00 0.615+0.015 13.2+0.4 7.8+0.8 
In 1.017 20.2 198+2 3530+ 100 





Values previously reported for the Mn resonance integral include one of 
5 barns quoted by Macklin and Pomerance (1956), 4.5+2.5 barns (Rose, 
Cooper, and Tattersall 1958), and 8.2 barns (Dancoff et al. 1947). For the 
In-115 resonance integral our value is to be compared with the value of 2640 
barns quoted by Macklin and Pomerance (1956) and one of 3500+250 barns 
obtained for naturally occurring In by Rose et al. (1958) using a pile oscillator. 
Using resonance parameters from the 1958 edition of BNL-325 the calculated 
value of D’ is 3050+150 barns. 


11. PREPARATION AND CALIBRATION OF Pu SOURCES 


Neutron capture in Pu-240, as in Au-197, is dominated by a single large 
resonance. The method of determining the Cd ratio by counting capture 
y-rays should therefore be applicable to Pu-240. The isotopic composition of 
the Pu samples used in the measurement of the Cd ratio of Pu-240 are given 
in Table VI. The Pu-239/241 sample was especially prepared to have the same 
ratio of Pu-239 to Pu-241 atoms as the Pu-240-enriched sample. 


TABLE VI 


Isotopic abundances of Pu target material 
(in atom per cent) 


Isotopes 




















Sample Ratio, 
designation 239 240 241 242 Pu-239/Pu-241 
Pu-240 11.46 88.01 0.49 0.034 23.4 
Pu-239/241 92.13 3.86 3.95 0.06 23.3 
TABLE VII 
Alpha counting of plutonium targets 
Specific Target 
Target a-Emitting Geometric efficiency Counting rate, activity, thickness, 
designation area, cm? of chamber c.p.m. d.p.m./mg  yg/cm? 
Pu-240(A) 1.267 8.306 X10-5 11,750 4.66 X 108 240 
2.85 8.306 X 10-5 26,590 4.66108 241 
Pu-240(B) 1.267 8.460 X 1075 12,160 4.66 X 108 243 
2.85 8.460 X 1075 27,320 4.66 X 108 243 
Pu-239/241(A) 1.267 8.757 X1075 597 1.50108 35.9 
2.85 8.757 X10-5 1380 1.50X 108 36.9 
Pu-239/241(B) 1.267 8.757 X10-5 405 1.50108 24.3 
2.85 8.757 X1075 910 1.50108 24.3 
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Two targets of each sample were used. The Pu-240 targets were scanned 
with a “pinhole’’ a-probe and found to be uniform within +15%. Alpha 
counting was carried out in a low-geometry ionization chamber to determine 
the absolute and relative amounts of the deposits. As an additional check on 
deposit uniformity the emission areas were restricted by disks with circular 
apertures of 3-in. diameter and 3-in. diameter. The results are shown in 
Table VII. 

The a-spectra of both Pu samples were investigated to determine the 
proportion of Pu-238 present. Concentrations were negligible in both cases, 
of the order of .01 at.%. 


12. DETERMINATION OF THE Cd RATIO OF Pu-240 


Several measurements were made using single counting. In addition to the 
Pu targets and the standard Au target a .00031-in. foil of Al was inserted to 
measure background counting rates. The Pu counting rates had to be corrected 
for the contributions of Pu-239 and Pu-241, as well as those of the background 
and the Al mounting foil. This contribution can be estimated from the net 
Pu-239/241 target counting rates (i.e. less Al foil and background rates) by 
multiplying by a factor equal to the ratio of the relative abundance of the 
Pu-239 and Pu-241 in the Pu-239/241 target to that in the Pu-240 target 
multiplied by the ratio of the Pu thicknesses (ug/cm?) in the two targets. This 
factor is of the order of unity for both Pu-240(A) and Pu-240(B). The results 
of all runs are shown in Table VIII. 


TABLE VIII 
Cadmium ratios of Pu-240 and Au-197 


Pu targets used Cadmium ratios 


Pu-240 = Pu-239/241 Au-197 Pu-240 (So) Pu-240 





A 3.38+0 .067 2.31+40.035 36.9+1.3 
A, B A, B 3.3340 .067 2.31+40.030 36.3+1.2 
B 3.41+0.068 2.34+0.040 36.74+1.4 


A,B A 3.31+0.066 2.28+0.025 36.9+1.2 
B 


Equation (5) is used to calculate the values of so for Pu-240 shown in the 
last column. For T = 130+10°C the tabulated value of ‘‘g’’ for Pu-240 
(Westcott 1958) is 1.060+0.003; it is dependent only on the location of the 
large resonance (E, = 1.056 ev). For the thicknesses of Pu-240 used in this 
experiment the self-shielding factor ‘‘G,’’ is approximately 1 and is not too 
sensitive to the resonance parameters. The value used in calculating so was 
(Gi) pu-2co = 0.983. The factor calculated for shielding by the .030-in. cadmium 
filter at the Pu-240 resonance energy (F) is 0.94. The mean value of So is 
36.7+0.7. 

If neutron capture takes place predominantly through a single energy state, 
as in Au, the energy distribution of the y-rays should not change when the Cd 
filter is inserted. For gold this was verified by comparing the Cd ratios deter- 
mined by capture y-ray counting at various biases and by activation; for 
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Pu-240 only a check at different biases is possible. The Cd ratios were found 
to be constant within experimental error for single y-ray counting at several 
biases ranging from 2 Mev upwards. (The mean value of these measurements 
is shown as run No. 4 in Table VIII.) 


13. THE THERMAL CROSS SECTION, RESONANCE INTEGRAL, AND 
RESONANCE PARAMETERS OF Pu-240 
Our value of so for Pu-240 can be combined with other measurements 
involving irradiations in NRX to yield a value of oo, the 2200 m/s cross section 
for Pu-240. The epi-Cd resonance integral, 2¢q, has been measured by Cornish 
and Lounsbury (1956) as 8700+800 barns. For the thickness of Cd used in 
their experiment 2’ = Yeq—0.45e0, so that 
_ 2 (Sca—180) 
ft So 
In addition, two measurements of the effective capture cross section for 
Pu-240 have been made; one in an empty fuel rod position (Fields et a/. 1956) 
(¢ = 530+50 barns), the second inside a fuel rod (Craig et a/. 1958) (¢ = 1000 
+50 barns). The values of rx/(7'/T») for these two positions are known to be 
0.026+0.002 and 0.070+0.006 respectively. The thermal neutron temperature 
is taken to be that of the moderator, 38° C, for which g = 1.03 (Westcott 1958). 
The two values of oo, calculated using equation (1), are listed in the last two 
lines of column 5, Table IX.? 


= 264 barns. 


i) 


TABLE IX 
Pu-240 parameters 





Tres, 70, >’, 

Source I, milli-ev kilobarns barns kilobarns >'/oo 
Fluharty (1957) 34.543 176* 295t 8.5 28 .9T 
Coté et al. (1958) 34.543 173 290t 8.4t 29 .9T 
Leonard, Seppi, and 

Friesen (1958) 32.341.5 — —_ — 30. 8t 
Egelstaff et al. (1958) 40+3 130+5 300t 7.45f 24.8t 
Ioffe and Okun (1956) 42 190* 480t 11 .3f 23 . 6T 
Krupchinsky (1957) -- -- 460 +45 10.0+2.8 21.7 
Rose, Cooper, and 
Tattersall (1958) - — 370+40 11.3+1.0 30.5 
Cornish and Lounsbury 

(1956) —_ — _ 8.5,+0.8 -- 
Pattenden (1958) — — 273248 — _ 
Halperin and Stoughton 

(1958) ae — 285 +25 —_— — 
This paper 30.6+0.6* — 264+30f — 32.5+0.6 

278+30§ 
268 +30)| 





*Calculated from other parameters. 

tCalculated from resonance parameters, including contributions from higher resonances (see text). 

{Derived from measured value of 2’/eo, using 2’ from Cornish and Lounsbury (1956). 

§Derived from 2’/oo and @ (Pu-240) = 1000 +50 barns in an NRX fuel rod (r VT/To = 0.070+ .006). See Craig 
et al. (1958). 

||Derived from 2’/oo and ¢ (Pu-240 = 530+50 in an NRX empty lattice position (r WT/To = .026+ .002). See 
Fields et al. (1956). 


This table supersedes that given by Westcott et al. (1958); some of the corrections (notably 
F) described in Section 5 above had not then been applied to the results. 
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The mean value for the 2200 m/s capture cross section is 







oo = 270+17 barns 






and for the resonance capture integral is 







>’ = 32.5X270 = 8.78+0.55 kilobarns. 





If a single resonance dominates capture in both the thermal and epithermal 
region it can be shown that 











o FT, 





ie 70 FEY EY"{1— (Bo/Eo) 





Combining equations (10) and (7) for the 1.056-ev resonance, 





aa se ee pene Bey 
so = VT Jo L r fe , 
(So) 1.056 = ae 


It is known, however, that other resonances contribute to LY’ and go. Using 
published parameters for these resonances (Fluharty 1957) their contributions 
were calculated to be 150+50 barns and | barn respectively. Thus 


‘ ys wv 
2 (2 1.056-F —~ rem 


59 = 


) = (So) 1.0561 (so) rem 


VJ wT oo 


so that (so)1.056 = 36.7 —1/Wx = 36.1 and Fr = 30.6+0.6 milli-ev for the 1.056-ev 
resonance. 

In Table IX our values for 2’/o 9 and the derived parameters and cross 
sections are compared with recently published values. The first five lines give 
results of parameter measurements. In estimating 2’ and 2’/oo from these 
parameters the 150-barn contribution from higher resonances has been added. 
We only obtain full agreement with the time-of-flight measurements of the 
three groups of American workers (Fluharty 1957; Cété et a/. 1958; Leonard 
et al. 1958), and with the values of oo obtained by Pattenden (1958) and 
Halperin and Stoughton (1958). 


14. SELF-FILTERING WITH Pu-240 

A portion of the Pu-240-enriched material (Table VII) was alloyed with 
pure aluminum and from this stock two disks $ inch in diameter were fabricated 
for use as filters. Chemical analysis showed that the alloy was 12.6% Pu by 
weight, and that the thickness of the two disks corresponded to 13.0 mg/cm? 
and 27.4 mg/cm? of Pu-240. 

A number of measurements were made of the epi-Cd counting rates with 
either or both of these filters inserted at the filter position, and with Al filters 
having a similar aluminum content. The self-filtering transmissions for Pu-240 
resulting from these measurements are compared in Table X with values 









' 





“Jj 
qo 
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TABLE X 
Self-filtering in Pu-240 








Transmissions 
Filter thickness, ._——————— -——————— —_—— 
mg/cm? of Pu-240 Peis Teste 
13.0 0.247+0.015 0.26;+0.01 
27.4 0.165+0.010 0.16 +0.01 
40.4 0.134+0.008 0.13;+0.01 





calculated using equations (8) and (9), taking into account contributions from 
all known resonances (Fluharty 1957). For the 1.056-ev peak the values 
adopted were I, = 0.032 ev and I’, = 0.00245 ev. 

The agreement between 7 hneas and T eae gives added support to the parameter 
measurements of the three groups of American workers. It also indicates that 
there is no scattering or contaminant contribution of appreciable magnitude, 


as was the case with Au. 
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APPENDIX—SYMBOLS USED AND THEIR DEFINITIONS 


Symbol Definition 

Vo neutron velocity of 2200 m/s. 

Eo energy of a neutron with velocity v (.0253 ev). 

Eca the ‘‘effective’’ cutoff energy for a cadmium filter. 

E the energy at the peak of a Breit-Wigner resonance. 

o(E) absorption cross section for neutrons with energy E. 

To absorption cross section for a neutron with velocity vo (the 2200 
m/s cross section). 

c & = oo (g+rs), effective cross section ina flux having both Maxwel- 
lian and epithermal components. 

om the effective cross section in a pure Maxwellian flux distribution for 


a given temperature. 
Or the total resonance cross section (i.e. excluding potential scattering) 
at the peak of a Breit-Wigner resonance. 


> ee Bf. 2 4 
g f= bu/oo= = if ae Siow (4) exo( E)a kT) ° 





> >= 7 o(E) gE the resonance absorption integral. 
wkT E 
>’ LD’ =z — goo r- Eo dE , the resonance absorption integral 
J kr ee tts 


in excess of the effective Maxwellian cross section. 


E 
, Loa = E- foe) ae , the resonance absorption integral under 
p E 





a cadmium filter. 
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r = f (au/4) is the epithermal density index, where f is the frac- 
tion of the total neutron density in the epithermal distribution. 
the ratio of the reaction rate in an unfiltered neutron flux to the 

reaction rate in a Cd-filtered neutron flux. 


Ty 


alo do 


oe / To _ / 4 

So=S - ™ ae ° 

To = Eo/k = 293.6° K (20.44° C), “‘temperature”’ of a neutron with 
velocity vo (& is Boltzmann’s constant). 

K = Rea rV(T/T>) for a 1/v detector (g = 1, 59 = 0) (see text). 

the full-width of the resonance peak o(total) = o,/2. 

r,, Tl, the partial widths for radiative capture and resonance scattering 

({! = f,+T, if no other processes are involved). 





ANALYSIS AND SYNTHESIS OF RADIATION PATTERNS 
FROM CIRCULAR APERTURES'! 


A. IsHIMARU AND G. HELD 


ABSTRACT 

Part I considers the problem of determining the source distribution over a 
circular aperture required to produce a prescribed radiation pattern. In particular, 
the problem of optimizing the narrow broadside pattern from a circular aperture is 
discussed in detail and an improved design method over Taylor’s for line source is 
devised. Numerical examples are given. 

Part II deals with the analysis of the radiation pattern from a circular aperture 
from 7 to r2 with the traveling wave type source functions. Expressions suitable 
to the analysis and the synthesis are obtained and the narrow-beam and shaped- 
beam synthesis are discussed. 

INTRODUCTION 

The problem of determining the source field over an aperture in order to 
produce a required radiation pattern has been considered for the rectangular 
aperture by Woodward in 1946. Taylor (1955) extended the method to include 
a design which approaches the optimum pattern. However, the synthesis of 
circular apertures has not yet been discussed except for the method outlined 
by Johnk (1954), who employed the conventional source function of the form 
(1—r*)? and as a consequence, had to use the Schmidt orthogonalization 
process, which is rather difficult to apply to an actual synthesis problem. 

In this paper we investigate the synthesis problem for circular apertures by 
proposing to represent the source function by a series of Bessel or Hankel 
functions. This representation yields the radiation pattern in a form similar 
to the one encountered in the rectangular aperture case, and a technique 
similar to Woodward's and Taylor's may be employed for the circular aperture 
case. For broadside antenna we devised a method to improve Taylor’s design 
of optimum pattern. For radiation at an arbitrary angle, we employ the 
traveling wave type sources and discuss the required modifications. 

In Part 1, we discuss the constant phase broadside pattern by representing 
the source by a series of Bessel’s functions and we clarify the relationship 
between the radiation field and the circular aperture source and develop the 
general form of the radiation function. We then define the optimum pattern 
for aperture antenna, formulate this optimum condition, and devise a method 
to realize this optimum pattern. Numerical examples are given and compared 
with the conventional design and Taylor’s design. 

In Part II, we discuss the radiation from a source composed of traveling 
wave type modes. The source is expressed by a series of Hankel functions 
over the range from 7; to 72. The radiation field from this source is investi- 
gated in terms of sidelobe level, beamwidth, and the choice of 7; and r2. The 
synthesis of narrow beam and shaped beam is also discussed. 


‘Manuscript received August 6, 1959. 
Contribution from the Department of Electrical Engineering, University of Washington, 
Seattle 5, Washington. 


Can. J. Phys. Vol. 38 (1960) 
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PART I. CONSTANT PHASE 


1. Source Function and Radiation Pattern 
The radiation pattern from a circular aperture is given (Silver 1949) by 


(1) g(u) = J flr) Jour) rér 


for a source independent of ¢’, where (see Fig. 1) 


Fic. 1. Circular aperture source (p, ¢’) and the field point (R, 9, ¢). 


g(u) = radiation pattern, 


f(r) = aperture field distribution, 


= — on 6, 


a = radius of the aperture. 


The source function which has been employed in the past (Silver 1949; 
Johnk 1954) is of the form 


f(r) = -r’)?, p=0,1,2. 


This function yields the radiation pattern of the form 


J 
g(u) = 2p) 2) 


From the viewpoint of the synthesis, this pattern is quite inconvenient 
since this is a non-orthogonal function and the only way of utilizing this is 
through the Schmidt orthogonalization process which is hardly practical 
(Johnk 1954). Also this offers no means of optimizing the pattern which is 
the main concern of this Part I. 
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For the purpose of synthesis, we use the source function of the form 





N 

(3) f(r) = oe AnJ (Un?) for0 <r <1, 
n=0 

= 0 forr > 1. 


Each term of (3) is chosen to satisfy the homogeneous boundary condition 
at r = 1. Thus 


(4) UnJ (Un) thJo(un) = 0 


where h is a constant. It may be seen that if h — infinite, each term of (3) 
approaches zero at r = 1 and f(1) = 0. If & = 0, each term of (3) has maxima 
or minima at r = 1. As will be shown later, the optimum pattern requires 
h = 0, and this case will be discussed in detail. However, physical structure 
of an antenna sometimes requires the vanishing source field at the edge and 
in this case, # — infinite should be used. 

Substituting (3) into (1) and employing the Lommel integral formula 
(Jahnke and Emde 1945), 


J Four) Folur) rd = yy [uw o(tgt) Je ur) — thy Ti(war) Jo(ur)] 


we obtain 


(5) g(u) = Dan) fu Ji(u) +h Jo(u)- 


The value of a, to produce g(u) can be determined by the value of g(u) at 
u = u, in the following manner: 
If u approaches one of u,, say ux, then the function 


uJ6(u)+hJo(u) 
2 


u.—Uur 
is zero except when uy = u,. When u, = u,, this function becomes 


2, \ tun 
Film) SEP. 


Therefore, a, is determined by g(u,) by the following equation: 


Que 
(6) Om = Gta) Tn) BM) 


In order to utilize eq. (6) and to determine the proper a, to produce a 
desired pattern g(u), it is necessary to investigate in what form g(u) should 
be represented. 
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For this purpose let us write (5) in the following form: 








(7a) glu) = lu So(u) +hJo(u)) De ag 22, 
(7b) = [uJ5(u)-+hJo(u)] WE naa | 
\? 
My, [-(¢) | 
[uJo(u)+hJo(u)] Il | 1-( roe s 
(7c) sie cece n=! Gall 


i -()) 


Py(u?’) in (76) is a polynomial in uv? of order N and can be determined by the 
location of its N zeros. The N zeros of Py are expressed by u,+t,, » = 1, 
2,... Nin (7c). The constant is chosen to normalize g(u) such that g(0) = 1. 


2. Optimum Pattern 

In order to obtain an optimum pattern, we list the following statements: 

(a) t, Must Be Real for an Optimum Pattern 

It is well known (Silver 1949) that the best pattern can be obtained only 
when all the zeros lie on the real axis, and that any complex zeros tend to 
produce higher sidelobe levels and broader beamwidth. Therefore, the zeros 
of Py must be real, and since u, is real, t, is also real. 

(6) The Constant h Must Be Zero for an Optimum Pattern 

In (7c), we note that for the range u > wy4;, the zeros of g(u) are Uy, 
n = N+1, N+2,... in eq. (4). In the range u < uy41, we wish to place NV 
zeros U,+t, so as to produce a desired pattern. Now, the optimum pattern 
requires the minimum beamwidth for a given sidelobe level. This can be 
achieved only when the zeros are placed as near to the origin as possible 
while keeping the sidelobes at a given level. 

It can easily be seen from eq. (4) that as h becomes smaller, the zeros u, 
for n > N+1 move towards the origin. But the negative values of h can 
not be used because the first zero is lost and the mth zero for h < 0 is farther 
from the origin than the mth zero for h > 0. Thus, the optimum condition 
requires h = 0 and the first zero uo coincides with the origin. 

The radiation pattern (7c) must be written then as 


«Th [1-(4) ] 


where u, is a positive root of J;(u) = 0. 
It may be noted from eq. (3) that the source field is finite at the edge 
r = a for the case h # ~. If h-+@ may be used, the source field always 


approaches zero as r > a. 





(7d) g(u) = 
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(c) Definition of Optimum Pattern 

Optimum pattern is defined as the one having a minimum beamwidth for 
a given sidelobe ratio (or a maximum sidelobe ratio for a given beamwidth) 
for a given N where N —1 is the number of terms used to describe the source 
function as in (3). The significant difference in the definition of the optimum 
pattern between the linear array and the aperture lies in this additional 
condition concerning N for the latter. For a given N, we have a freedom 
of choosing N zeros or u,+#, in (7d) to make the pattern optimum. Since 
we have only W zeros to adjust, the best we can hope to obtain is the equal 
level for the first NV sidelobes. The rest of the sidelobes are tapered down 
according to the equation (7d) and we have no control over those sidelobes 


(see Fig. 2). 





U,+t, : Ut ta: Uy+ts, Us Us Us 
' ‘ 
M, Me Ms 


Fic 2. Optimum pattern (NV = 3). Equal level for the first N sidelobes. Tapered sidelobes 
for u > uy+. 


(d) N Must Be Finite 

It may be obvious from Fig. 2 that the larger the .V, the closer the zeros 
Un +t, can be chosen towards the origin for a given sidelobe level. It seems then 
that the infinite NV or the infinite series for the source would achieve the truly 
optimum pattern. That this infinite series is physically impossible was shown 
by Van Der Maas (1954) and Taylor (1955) for the line source. It is shown 
in the Appendix that a similar argument can be made for the circular aperture. 
If the source function is an infinite series, the equal sidelobes are extended 
in all the region u,; < u < © and this requires the value of the source function 
f(r) at the edge r = 1 to become infinite. Therefore in order for f(r) to be 
finite, NV must be finite. 

(e) Minimum and Maximum N, Supergain 

In order to find the practical range of N, let us consider a pattern as shown 
in Fig. 3. Let us consider the design for 30-db sidelobe level. If we employ 
N = 2, we can choose the first two zeros to get 30-db level for the first two 
sidelobes. However, the third sidelobe cannot be controlled and though the 
level may be slightly lower than 28 db, the third lobe may remain above 
30 db. Therefore N must be at least 3, for then the fourth lobe has 31 db or 
lower and the situation that all the sidelobes are equal to or below 30 db is 
established. In general, the minimum N must be chosen such that the (V+1)th 
sidelobe of 2/,(u)/u has the level below the desired sidelobe level. 





i ie 
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PAG) 
2-0 





(Cc) N=3 
Fic. 3. Choice of N. For 30 db, N = 2 is not enough, N = 3 is minimum. 


The maximum N is determined mainly by the supergain consideration. If 
N is chosen so large that u, lies in the invisible region, then g(u) in the region 
(2ma/X) < u < Uy41 must be increased to attain the prescribed sidelobe 
level. Thus, this leads to larger reactive power and a supergain antenna. 
Therefore the maximum N must be such that uy is in the visible region. 

In fact if wy is in the invisible region, we can place all N zeros in the visible 
region, thus obtaining a smaller beamwidth than the optimum design in this 
paper for the same sidelobe level. However, since this requires excessively 
large reactive power in the region (24a/A) < u < uyy; this antenna is not 
practical. 

Between the above two limits, V must be determined by practical considera- 
tions such as the amount of computational labor (later it will be shown 
that the solution of N linear equations is required) and the physical character- 
istics of the source, which may not be able to produce too high an harmonic. 


38. Optimum Condition and Solution 

Our next task is then to find a method to determine the zeros u,+#, in (7d) 
in order to realize the optimum pattern. To this end, we formulate the optimum 
condition in the following manner. It is seen from the previous section and 
Fig. 2 that there are N sidelobes in the region u,; < u < uy4;. Assume that 
the peaks of those sidelobes occur at u = M,, m = 1,2,3,...N (see Fig. 2). 
Then. the optimum condition is expressed by the following N equations: 


(8) lg(Mi)| = |g(Ms)| = |g(Ms)| =... = le(My)| = 6 


where 3 is the sidelobe level and g(u) is shown in (7d). M, are functions of 
ty, t2,...t,. Therefore, this is a system of N equations for N unknowns 


bey Basins «Bue 
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The solution of (8) seems rather difficult, but the following procedure is 
found quite effective in solving a problem of this type. Let us write equation 
(7d) in the following manner: 


I+ pt hy to = 
2Ii(u) EE _ a ‘| 


S Lie | 


Let us for a moment assume that ¢, are small compared with u,, and that 
#2? terms are negligible. Now in eq. (8), M, are located at about the middle 
point between u, and u,41. Therefore, for small ¢,, we approximate the product 
in (9) by a series, 





(9) g(u) = 





‘ x Qu 
2Ji(Mn)| 1+ 2) =F tn 








(10) g(Mn) = eee 
N 9 
Mal 1+ z= i 
n=l Un 
Therefore, eq. (8) becomes 
Pe 
(11) ya Qnatn = CC. m= DD x athe 
n=1 
2un M,,b 2 
(i) Onn = 7M |2T,(Ma) 
Cm = = =~]. 











The equation (11) seems to be N linear equations for N unknowns. However, 
there are two points which must be clarified in order for (11) to be treated 
as a linear equation. 

(a) t, Must Be Small 

The equation (10) holds only for small t,. However, this is not an essential 
difficulty, because this process can be repeated and after the repeated process, 
t, becomes so small that equation (10) holds almost exactly. The repeated 
process is as follows: 

Suppose we solved (11) and found ¢,; g(u) in (7d) is not yet optimum 
because ¢, may not be small. To find the pattern closer to optimum, let the 
zeros be given by the u,+/,+#}. Then 


Hh L- Geta) ] 
I r ()'] 





(13) g(u) = gi(u) 
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TM L-(¢) | 


gi(u) is obviously an improvement over 2/J,(u)/u. And if we compare (13) 
with (7d), we see that (13) has the same form except that [2/,(u)]/u and u, 
in (7d) are replaced by gi(u) and u,+¢, respectively in (13). Therefore, to 
get the better pattern we use 


where 


(14) gi(u) = 


N 
(15) z. a ee m= eae cole 
n=1 


where @m, and cp, are shown in (12) with u, and 2J;(M,»,)/M,» replaced by 
Unt+t, and gi(M,,). 
Solving (15), we obtain ¢} and we get the improved pattern 


2n(u) TI Sle Sea 
(16) g2(u) = net bath ttn 


ly 
u 1—\— 
I] Un 
The further improved pattern may be obtained by letting the zeros be given 
by u,+t,+ti+t)! and solving the equation similar to (15) for #2, and theoreti- 


cally this process can be repeated indefinitely to obtain the truly optimum 
pattern. However, for practical purpose, ¢} usually gives satisfactory results. 
(b) Choice of Mn 
The above discussion concerning the linear equation for ¢, is valid only 
when M,, is independent of ¢,. But M,, is actually a function of ¢,. This diffi- 
culty can be avoided by assuming M,, independent of t,. One such choice is 





_ Umt+Um+1 
(17) Mn = ore 


This assumes that the peak point of the sidelobe is at the midpoint between 
zeros. This causes an error but the error for the sidelobe of the pattern 
[2J,(u)]/u is less than 0.4 db. Considering the great simplification made 
possible by the assumption (17), this unavoidable error of 0.4 db is certainly 
tolerable. 


4. Numerical Example 

As an example, let us consider a design with a sidelobe level of 25 db. As 
is seen from Fig. 3, the minimum JN is 2. However, since three linear equations 
can be easily solved, we take N = 3. Obviously, the higher V will give the 
better pattern. Now we employ the following constants for eqs. (11) and (12). 
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= zeros of J;(u) M,, = maximum point of J;(u)/u 
= 3.83 M, = 5.33 
= 7.016 M2 = 8.53 
us = 10.173 M; = 11.70 
ee 
$8 


Solving the linear eq. (11) we get gi(u) in (14) 


ano] (5) L1G) L- (aos) | 
1-(s45) L1-G4s) -Gta) | 


This gi(u) exhibits the obvious improvement as shown in Fig. 4, but its first 
sidelobe is still higher than the desired 25 db. This is because of the relatively 


b (25 db down) 


(19) gi(u) = 





s—= sind 


Fic. 4. Example of synthesis for 25-db sidelobe level. 


large ¢, for which the equations (10) and (11) hold only approximately. The 
next attempt in (15) must improve this pattern considerably, because ¢} is 


now very small. 
After computing #2 in (15) we obtain the improved pattern 


41-Gi5) L1-Gés) L-Gia) ] 


In Fig. 4, the improved pattern g,(u) and the further improved pattern g2(u) 


are shown together with the pattern 2/,(u)/u. The half-power beamwidth 
can be approximately given by 


(20) ge(u) = 
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Awmttht+ti 
oa radians. 


D 
In Fig. 5, the synthesized pattern is shown with the conventional 24.6-db 
pattern 8/2(u)/u. Both have the same sidelobe level, but the marked improve- 


(21) 6 < 1.02 






Fic. 5. Comparison between the optimum and the conventional pattern for the same 
sidelobe level of 25 db. 










ment of the synthesized pattern is obvious in the reduction of the beamwidth. 
The source function can be found by (6). Thus 








(22) f(r) = 1+0.912 Jo(3.83 7), 
+0.595 Jo(7.0167), 

— 0.987 Jo (10.173 r) for O:<S-r-< Ft, 

= 0 for r > 1. 










5. Comparison with Taylor's Design 

Taylor's design procedure for a line source can be adapted to the circular 
aperture in the following manner: 

For u > uy, the zeros are roots of J;(u). 

For u < uy41, we employ the zeros which are similar to that of the ideal 
pattern 









cos o/(u?—A”*) 
cosh A 





: cosh A = sidelobe ratio. 


(23) g(u) = 
Thus 


(24) 
where 











Un = ovV/[A?+(n—1/2)*x?] for mn = 1,2,...N, 









Un+1 


7° VIA +(W+D—1/20 8] * 
For the example of the previous section, 
N =3, A = 3.57, o¢ = 1.11. 


(25) 
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So, we obtain 

u, = 4.34, 
(26) u2 = 6:7. 



























The pattern is given by 


C7) gee) = —— : 
Uu Uu Uu 
| 1-(s%5) [-Gs) [1-Gets) | 


The first zero is at 4.34 while the optimum design yields 4.2885 which gives 
a narrower beam and this improvement of the optimum design over Taylor’s 
design is due to the fact that Taylor’s design does not insure the equal side- 
lobe level for the range “4; <u < uUyy41. 


PART II. TRAVELING WAVE TYPE MODES 
1. Formulation of the Problem 

Let us consider a circular aperture in a perfectly conducting plane (Fig. 6). 
The electric field is assumed to lie in a plane of aperture and to be directed 
in ¢o direction. Furthermore, we assume no variation of field in @9 and the 





Field point iz 


(REO) Arig 
Os 









Source point ie 
(er ¥ #) oS 
Fic. 6. Source and field points. 





fields are functions of radius only. This choice of field configuration is made 
on the basis of a recent experimental investigation on a radial surface wave 
antenna that was performed in our Antenna Laboratory. However, the method 
to be described in this paper can be applied equally well to other configurations 
with little modifications. 

The Huygens’ principle of the electromagnetic vector problem gives the 
necessary relationship between the source electric field and the radiation 
fields. The electric field produced by a given electric field over an aperture 
in a conducting plane is given by (Lewin 1951) 

—jkR | 


ds. 


é 


R 





1 
(1) E= tux (xz) 










: 
i 
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The far-zone radiation field is then given by 


(2) Es = _ c= J (—nXE)-i, e*” '"'8 ds 
3 ha = oe —IkRo i, pdr i-ig 
8) i? eae (nXE) -ig e” 8 ds 


where Nn is a unit vector outward normal to the surface S, and i,, ig, is, i, are 
unit vectors as shown in Fig. 6. Since E = Ei,,, we have 





(4) Ey, = se eo f E sin (¢— $0) e”*” °°" ds 

- jk —jkRo jkr cosy 
(5) Ey = ane’ E cos 6 cos (¢@— go) e€ ds 
where 

(6) cos y = sin 6 cos (¢—¢o) 


ds = r dr dd. 


Since E is a function of r only, the integration with respect to ¢o can be 
readily performed to yield 





E, =0 
(7) E, = a e 7**° cos 6 g(6) 
0 Gal é 
r2 
(8) g(6) = f f(r) Ji (Rr sin 6) rdr 
r1 


where f(r) = E(r) is used to conform with the notation in Part I. It is seen 
that the radiation pattern is a product of cos @ and the function g(@). Our 
problem of synthesis is then that of finding the source function f(r) required 
to produce a prescribed radiation function. 


2. Source Functions 
Two different forms of the source functions may be considered and each 


one demands a different approach to our synthesis problems. One is a source 
function with constant phase and varying amplitude. The other is the function 
with varying phase and varying amplitude. Both functions are capable of 
producing the desired beam at an arbitrary angle, and some general aspects 
of these problems will be discussed here. 

(a) Constant Phase 

The source function f(r) takes a form 


(9) f(r) = SY anJ1 (cnr) forO<r<a 


n 


= (0 forr >a 


where a, are real. Each term is similar to the different modes in a circular 
guide. In general, the choice of a, may be made such that 


(10) Ont] { (amt) +hI (ana) = 0. 
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The magnitude of the source function at the edge depends on / as in Part I. 
This problem will not be treated in this paper since the situation is similar 
to the Part I. 

(b) Varying Phase 

It is often desired that the source is composed of the traveling wave type 
modes. Our particular choice is 


(11) f(r) 











Dd apNnH{” (anr) forri <r < ro, 







0 forr <7, 





PS ts 


The Hankel function H{” is chosen because of our choice of time variation 
e', so that H\” represents a traveling wave in a positive r direction. N, is 
a complex number which gives the phase shift between each term and also 
serves to normalize the radiation pattern; a, is a real number representing 
the strength of field contained in each term. 

In the subsequent sections, an analysis and a synthesis procedure for the 
sources represented by (11) will be explained in detail. 










3. Analysis of Radiation Pattern 
We wish to examine the radiation function 





(12) g(0) = WOT sin 6) rdr 





for the sources given by eq. (11). First, we take one term of the series in 
eq. (11). Then, we get 





(13) g(0) = aN J H (ar) J,(kr sin @) + dr. 





At this point let us choose 7; or a, such that 





(14) nr, > 1. 


This does not restrict the applicability of our solution, because as will be 
seen later, H‘? (a,r) yields the main beam at @, where & sin 6, = an, and we 
can always choose r; to satisfy (14). Under this condition (14), the Hankel 
function in (13) can be approximated by its asymptotic form whenever it is 
convenient. 

The integration in (13) can be easily performed by the Lommel integral 
formula. However, this integration does not give us a form suitable to our 
analysis and synthesis problems for the useful range of 6. Thus, we consider 
two cases, one dealing with the useful range of 6, 0g < @ < 2/2, and the other 
dealing with the angle @ nearly equal to zero, 0 < 6 < &. 

(a) 04 <0< 4/2 

In this range of 6, we employ the form 


(15) (8) = 6 [g, (0) +e # g0(8)], 
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(16) e””9,(0) = Salts J HH (cogr)H® (kr sin 6) r dr, 
2 i 
, ere 
(17) 7? 9,(9) = OnNn J H}” (aar)H}” (kr sin 0) r dr, 
2 a 


where we employed the relation 


(18) Ji(x) = 3[AY(x)+H1?(x)] 


and the phase functions p(@) and ¢(@) are inserted to simplify the expressions 
later and are independent of m. The integration in eqs. (16) and (17) can be 
readily performed, and if we choose 6, such that 


(19) kri sin 0, > 1 


then all the Hankel functions can be approximated by their asymptotic forms. 
The results can be written as 


= My sin (u— tn) 
(20) g1(u) = dn of * ecu 
) 
- Up SIN (U+Un) 

(21) g2(u) = a, V 2. ade. 
(22) p(6) = aT 4 dae 
(23) (6) = 4% v5 sin 0-42 

qg eo ro sin 9° 
where 


ro = 1/2(r1 +12), 
r= 1/2(re—n1), 


(24) “= = ra sin 6, 
Un = Alay 


and the normalizing constant N, in (16) and (17) has been chosen as 


ai Ant janto 
(25) Nn Papin ec. 
The visible region of u is then 


2 Ta 


r 





(26) O0<u< 


It is important to note that there is a one-to-one correspondence between 
the intensity a, and the value of gi(u) at u = u,: 


(27) £1 (Un) = An. 
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The physical interpretations of eqs. (20) and (21) may be made by considering 
the radiation pattern produced by a traveling wave on a rectangular aperture. 
It is well known that the aperture field in a rectangular aperture expressed by 


(28) frlx) = € 7", —x—g < x < x 
produces the radiation pattern 


sin (u—Un) 


at gi(u) = u—uU 


where 


“= — X4 sin 6, Un = OnXq. 

We can see that eq. (20) has exactly the same form as (29) except / (u,/u). 
Therefore, gi(u) is similar to the radiation from the source wave traveling 
in a direction of the observation point. Similarly, g2(u) represents the radiation 
from the source wave traveling in the opposite direction. 

The radiation from the circular aperture is a combination of these two 
gi(u) and ge(u) as in eq. (15). It is understandable that since the field over 
the aperture travels in all directions around the axis r = 0, g(u) must be a 
sum of two waves traveling in opposite directions. (Fig. 7). 

gi(u) behaves like a sin x/x shifted by u, except for ~/(u,/u), and has a 
main beam at up; g2(u) is relatively small in the rangeO < u < u, (see Fig. 8 (a)) 


fi Id oint 
a ae 


Field point 


173, (u) 


(b) 


Fic. 7. Similarity between two rectangular apertures and circular aperture. 

(a) Two rectangular apertures. The one on the right side carries the source wave /i(x) 
traveling to the right. The aperture on the left carries f2(x) traveling to the left. 

(b) Circular aperture produces the radiation pattern similar to gi(6) and g2(@) from two 
rectangular apertures. The phase difference is (24/A)(2ro) + (4/2). 





ISHIMARU AND HELD: RADIATION PATTERNS FROM CIRCULAR APERTURES 93 


Therefore, the first approximation for g(u) may be given by neglecting g2(u): 


(30) \g(u)| = |gi(u)]. 


The maximum error introduced by (30) may be given by the maximum 
value of g2(u) which is 


{ u, 1 
(31) Ge u U+U,’ 


For u near u, or peak point, this error will be approximately 


(32) 





93Cu) 


—-—- ENVELOPE 


a A(uy 


30 4nT St) oT 


Fic. 8. (a) Radiation functions g:(@) and go(6). 
n=, ro = 7.92 X, ra = 3.46 A, 
6, = 60°, un = 6n, visible range of u < 6.92 7. 
(b) Radiation pattern g(@) of eq. (34) and envelope. 
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The detailed shape of the pattern g(w) may be found by considering 
(33) lg(u)| = {lgr(w)|? + |g2(w)|?+2lgrge cos ¥|}# 


where y is the phase difference between g; and ge. This phase difference is a 
function of u so that y increases by 2m as u increases by the amount (rq/ro)m 
as is seen from eqs. (15) and (23). Since rz < 70, as u varies by 7, ¥ undergoes 
at least a complete cycle of 27 or more. 

Therefore, the radiation function g(#) consists mainly of gi(u) and the 
effect of g2(w) is the small high-frequency components superimposed on g; (x), 
and g(u) varies between |gi/+/go| and ||gi/—|gel| (see Fig. 8 (d)). 


(b) 0<0< 6, 

In this range of angle, the asymptotic form for Hankel functions no longer 
holds and we must investigate the original equation (13). We note that for 
this range, gi(u) and g2(u) arrive at the observation point with almost the 
same magnitudes and little phase difference and it is no longer advantageous 
to treat the radiation field as composed of two waves traveling in the opposite 
directions. 

The integration in (13) yields after some manipulation 


(34) o(0) = Se%s, "a Hal ¢-slenn-D (4 — B) 
B —Qy v 
(35) A = Js (8r2)—j* Ju(Grs) 
(36) B= y/ 74 gien("2-1 J, (Br) —janJo(Br1)] 


where 6 = &sin 6. 

In eq. (35), B/a, is very small and A is approximately equal to J;(8r2), 
which has a peak value of 0.582 at Br. = 1.84. In eq. (36), the second term 
is also small and B is approximately given by 


(37) IB] = 4/ 2 J(6n). 


The approximate maximum value of |A —B] will be then 


3/2 
(38) |4—Blmax ~ 0.582-+0.92 (2) é 
2 


Therefore, the approximate maximum value of g(u) near @ = 0 is given by 


aan car 3/2 
(39) On / i / a | 0.582-+0.92 («) | 
Un Ya Ye 


(c) Dividing Line 
The angle 6, can be found by equating the peak value of (15) and (39). 


Thus 
(40) =| : | |- 4/ = V/2[ 0 582+-0.92 (=)""] 
os u utun! J Ge Qun Ya - ’ 1s f 


1 
U— Up, 








+ 
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From this, we obtain approximately 





(41) Ug=127 for 71 K ro. 


"1 
itz: 
(d) Analysis of Pattern 
It may be interesting to investigate the detailed shape of the pattern in 
view of the above analysis. In Fig. 8(a), gi(@) and g2(@) are plotted for 
r, = \ and ro = 7.92 and 6, = 60°. In Fig. 8(), the exact pattern com- 
puted from eq. (34) is shown together with gi(@) and the envelope given by 


Un 1 1 
/* l tata | for 0 > 64 


and eq. (39) for 0 < 6 < 6. 

It may be seen that the actual pattern consists mainly of gi(@) and the 
effect of g2(6) is the high-frequency components of small amplitude. The total 
pattern never exceeds the envelope given above. 





4. Source Function which Vanishes at the Edges 

Since we are considering the aperture on a conducting plane, it is required 
that the electric field at 7; and re be zero. In this section, we wish to examine 
such a source distribution and the radiation pattern it produces. 

Each term of the source function given by (11) does not vanish at 7r;, and 
ro, but a suitable combination of the terms will produce this vanishing field. 

Let us approximate eq. (11) by using the asymptotic form of Hankel 
function 





ll 


1 ts a 3 
ia AAS eden  forrs <r < 4a, 


n fe—F 5 r 
= 0 forO KR r<r, ra <r. 


(42) f(r) 


We combine the mth and the mth terms and letting 





An tam _ 
+o 9 = On 
(43) 
Qm — An, 
-_— 


~ 


and choosing a, and a, in such a manner that 


OnV On = OnV am = AV ap 





(44) 1 
ar = (2k+1)5——, k =0,1,2,.. 
d 
we get 
ae ge nr—Po 
(45) f(r) =K 3 p> a, cos (2k+1) ; 


where 
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(46) K = constant, 


a 2 V/2 Tay eltrrot si r 
fo—?1 F 


The series in eq. (45) is symmetric about r = ro and vanishes at 7; and 
ro. If a, is real, eq. (45) represents a traveling wave with a constant velocity 


and its amplitude varies approximately as 1/1/r. 
The radiation field arising from such a source distribution can be easily 


obtained: 


— 2(2k+1) 5 (—1)** cos (u—upy) 
(47) gi(u) = » ay 2 pie 


(u—u)*—| (2k+1) = 


— 2(2k+1) 5 (—1)*** cos (w+u) 
go(u) = DO ao /# ; senaeen 


(u+-u)*-| (2k-+1) z) 





(48) 


2 


There is a one-to-one correspondence between a, and the value of g;(u) at 
U = Up (2k+1)x/2. 


. / u 
(49) gf uy 2k+1) =| =a, / et 7 
/ up (2R-+1) 5 
Also, Vu g:(u) is an even function of (u—w,). 
5. Narrow Beam 
If we desire a narrow beam directed at an angle @,, it is most convenient 


to choose 


. 2n . 
(50) a, = - sin 6,, Up = Pa. 


The general form of the radiation pattern gi(u) can then be written as 


. a a Pryl (u~tp)’] 
(51) gi(u) = / * cos (u —Up) gee ae 
fh fi-(o*, 
(2k+1) 5 


where Py[(u—u,)*] is a polynomial in (u—u,)* of order NV. This can also be 
written using the locations of zeros as follows: 





gi(u) = 4/ cos saci IT [(-"=*)'] 


(52) = : ui 
N ar 
& Dp 
IT 1 Sain Tv 
(2k--1) 9 


where », are zeros of gi(u). 
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Since the magnitudes of g2(u) are relatively small as compared with g,(w), 
we can neglect g2(u) as a first approximation. Then equation (52) is similar 
to the radiation pattern (7d) of Part I, and the same technique may be used 
to optimize the radiation pattern. It is also possible to include the effect of 


















nd 2(u) by considering the radiation pattern |gi(#)|+-|go(%#)|max in such a manner 

g Pp g g 
ty that the actual sidelobe level is always equal to or lower than the prescribed 
| level. The detail of the treatment of this point is given in the previous report 
ly by Ishimaru and Held (1959). 

(a) Sidelobe Level 
If the source function consists of one term for k = 0 in eq. (45) we have 

(53) 

This is the simplest case and the source field is 1/+/r times the half sine 

wave. 

The first sidelobe level is approximately at «=u,+2z and the maximum 

magnitude is given by 

. Up 1 1 
; (54) les) Lea(0)lnas = 00 4/22 | ——! _ wigs aie ail go 
L 

2 Tv 
| (w—u)*-(2) (u+u,)*-(Z) | 
If u > u,, then this becomes 
4 1 

55 a aio8 Jt, 

(55) . Iba? due 

The sidelobe ratio is 

" ] 1 ey 

56 cnet ae. ae © 
t We 15) dae 2 





As u, becomes large, this ratio approaches 1/15, which is 23.52 db. 

(6) Beamwidth 

When 6, is not near +/2 and the beam is sufficiently narrow, cos @ in eq. 
(7) is practically constant. We obtain from the equation 






(57) cos (u—4uU,) 4; 45 
Lach Seca F 
Tr 


= 
“ 









the beamwidth in terms of u, 
(58) 


The beamwidth in terms of @ is 







(59) 
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If @ is near r/2, it becomes more difficult to obtain the sharp beam. The 
beamwidth must be computed from 


cos@ cos(u—u,) 1 
(60) ed ae 


cos 4, 1 et 
v 


L 2 
(c) Choice of r; and rq 
r,; must be sufficiently large so that eq. (14) holds 


(61) kr, sin @, > 1. 
rz must be chosen so as to give a desired narrow beamwidth. From (59) 
(62) ra 3.74 


\ ~ 2m cos 6,:A0 


where A@ is a desired beamwidth. 


6. Shaped Beam 

As stated in Section 4, »/u g:(u) is an even function of (u—u,). Therefore, 
by the source function described in Section 4, it is not possible to obtain an 
arbitrary function of u. For example, if we wish to obtain the cosecant pattern 


(63) | cos 6g(@)| = cosec (=~<) for 0, < 6 < 02, 


< 


nil 0<0<6,and <4<~, 


we cannot choose u, arbitrarily. Instead, we must choos* u, so as to produce 
this desired pattern. It is possible to choose u, such that 
2ar 
(64) a ee 
r 
Then, all the visible region of u is on one side of u, and we could assign any 
prescribed arbitrary pattern in this region. Of course, the other side of u, is 
in the invisible region of u and this would tend to increase the reactive power. 
If we neglect ge(u) in (63), we can find the source function to approximate 
the radiation pattern in exactly the same manner as the Woodward method 
(1946). 





CONCLUSIONS 


In Part I, the synthesis of circular aperture antenna is considered with 
particular emphasis on the optimum pattern. Taylor’s narrow-beam synthesis 
technique for a line source is applied to a circular aperture and the improve- 
ment over Taylor’s design is shown together with the definition and the 
limitations of the optimum pattern. 

In Part II, the circular aperture with sources composed of traveling wave 
modes is analyzed and its similarity with rectangular apertures is discussed 
in detail. Some considerations are given to narrow-beam and shaped-beam 
synthesis. 
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APPENDIX 









N must be finite. 
Let us consider the radiation pattern for circular and rectangular apertures: 


2J,(u) I Li-(45) | 


(1) wfL-@ 
smu Th L1-Ge) | 
oT [-( yi 


If N becomes infinite, both saianinia approach 


fi [-(a)'I. 


Therefore, we have the same ideal (but unrealizable) optimum pattern 


_ cos V(u"—A’*) 
(3) Bs 6) = Cosh A 


where cosh A is a sidelobe ratio. 
If this g,(u) is to be realized, the source function must be given by 


















(2) gr(u) = 




















Fey bil): 


Jo(ux 
For large k, 






_(-)* 
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uy, = ka+— and i) > i 
? 4 Be(u) V2 cosh A 
The source function is a series of a,Jo(u,7) and this will become at 7 = 1 
u 
eon) LED. 










As u, — ©, the term becomes infinite and f(r) at r = 1 is divergent. 
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THE FREEZING POINTS OF HIGH PURITY METALS AS 
PRECISION TEMPERATURE STANDARDS 


V. THERMAL ANALYSES ON 10 SAMPLES OF TIN WITH PURITIES 
GREATER THAN 99.99+%! 


E. H. McLAREN AND E. G. Murbock 


ABSTRACT 

Extensive thermal analyses have been made on 10 samples (suppliers’ 
analyzed impurity contents <0.2 to <100 p.p.m.) of high purity tin, including 
zone-relined metal; liquidus points have been intercompared with a precision of 
about 0.0002° C and alloy melting ranges have been examined following different 
types of freezing with and without overnight anneals near the solidus tempera- 
ture. Samples of nominal 99.9999% purity tin were found to have such narrow 
alloy melting ranges that any ambiguity, arising from unknown impurity con- 
centrations, in specifying the liquidus point of pure tin is well inside 0.001° C; a 
value of 231.913° C (Int. 1948) was found for the standard liquidus point of pure 
tin. An account is given of the supercooling that was observed on the bulk 
samples and of anomalous structures that were found on melting curves. An 
appendix gives the results of long-term intercomparisons of the temperatures 
realized in four water triple point cells. 


INTRODUCTION 

Earlier work* (Part II, 1957b) has shown that, with appropriate experi- 
mental conditions, it is possible to obtain highly reproducible plateaux that 
have nearly constant (within 0.0001° C) temperatures for long durations on 
the cooling curves of high purity samples of tin. The purpose of the present 
investigation was to increase our knowledge of the tin point by examining 
samples of the highest purity tin available from several refiners and to deter- 
mine a value for the liquidus point of pure tin on the International Temperature 
Scale. 


MATERIAL 


Table I lists the origin, grade, and form of the tin metal as received and 
Table Il shows the chemical analyses supplied by the sources. Attempts to 
estimate the expected lowerings of the liquidus from that of the pure element 
and the alloy melting ranges of the samples using these chemical analyses and 
known phase diagrams were abandoned due to the uncertainties in locating 
the liquidus for the tin-iron system at the tin-rich end and to the 
usual difficulties in locating the solidus for very dilute alloys. Campbell et al. 
(1949) report a eutectic transformation at about 0.003% Fe with the eutectic 
temperature not more than 0.1° C lower than the freezing temperature of the 
pure element; precision thermal analyses on a series of tin-iron alloys in this 
concentration range should yield an accurate determination of this liquidus. 


1Manuscript received October 1, 1959. 

Contribution from the Division of Applied Physics, Heat and Solid State Physics, National 
Research Council, Ottawa, Ontario. 

Issued as N.R.C. No. 5468. 

*Reference to preceding papers in this series is made by part number and date of publication. 


Can. J. Phys. Vol. 38 (1960) 
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TABLE I 


Identification of the tin samples 











Origin Grade Form 


Vulcan Detinning Co. Spectrographic, VS-52SP(1-5) Bar, 1.20.7 X7.5 cm 
Sewaren, N.J., (zone refined) 

U.S.A. Spectrographic, VS-52SP(51-55) Bar, 1.20.7 X7.5 cm 
(zone refined) 


Vulcan Detinning Co. Extra Pure, VS-65 Bar, 2.2X1.5X18.5 cm 
Sewaren, N.J., (electrolytic) 

U.S.A. Extra Pure, VS-24 Bar, 2.2X1.5X18.5 cm 
(electrolytic) 


The Consolidated Special Research (34-37, 42) Bar, 250 g 
Mining and (electrolytic) 
Smelting Co. 
Trail, B.C.,Canada High Purity, MR139 Bar, 2.5 X2.5 X25 cm 
(electrolytic) 


Johnson, Matthey Specpure. J.M. 530 Rod, 0.6 diam. X16 cm 
and Co. Ltd., 
London, E.C.1, 
England 


National Bureau of Standard Sample 42d Cast ingot, 
Standards, 3.8 diam. X8.0 cm 
Washington 25, 

D.C., UA 


Sn.12¢ ‘Philipp Bros. Zone Refined Rod, 1.6 diam. X22 cm 
(Canada) Ltd., (received Feb. 1959) 
Montreal 25, Que., 
Sn.10f Canada Zone Refined Rod, 1.6 diam. X22 cm 
(received Nov. 1958) 


*See Scribner (1952). 
tJapanese zone-refined tin. 


TABLE II 
Impurity contents of the tin samples supplied by the sources 





Weight, p.p.m. 





Pb Fe Cu i i In Ca As Mg Al 


<Q. n.d. 
<0. n.d. 
0.2 
<0.2 


5 nd. <1 
<50 <56 6 
tr. v.w. VW. 
Sn.10* tr. V.Ww. VW. 


3 
3 
2 

2 


<0. 


Note: n.d., not detected; tr., trace; v.w., very weak. 
*Qualitative spectrochemical analysis only. 


EXPERIMENTAL 


In the earlier investigation (Part II, 1957b) thermal analyses were made on 
only two samples of tin, Sn.2 and Sn.3; the liquidus point of Sn.3 was found 
to be about 0.008° C higher than Sn.2 and alloy melting ranges were found to 
depend both on the type of freeze and the duration of anneal prior to remelting. 
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(A) Thermal Survey 

The apparatus and techniques used in the thermal analyses have been 
described previously (Parts I-IV, 1957a, 6; 1958a, b); dry nitrogen atmos- 
pheres are now always used over the 1100-g tin melts which are contained in 
high purity graphite crucibles (manufactured and purified after fabrication by 
United Carbon Products Co., Bay City, Michigan, U.S.A.). 

Resistance thermometry for this survey was based on measurements of the 
2-ma resistance (Rj,) of the thermometer at 15cm mid-point immersion in 
the liquid metal corrected to 1 atmosphere nitrogen pressure (for pressure 
corrections see McLaren 1959) and single daily determinations of the zero 
current resistance (R%.,) at the triple point of water. For plotting and corre- 
lating results in this type of thermal survey it is a convenience to work with 
temperatures in arbitrary ratio units defined as Rj,/R%.,; for the Meyers 
thermometer, S155, used in this study 37 X 107° ratio units ~ 10-*° C. 

Figure 1 shows typical cooling curves (Sn.3) that are obtained when a 
sample of high purity tin transforms under different rates and types of freezing. 


MINUTES FROM END OF SUPERCOOL 
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Fic. 1. Typical freezing curves obtained on high purity tin. 


The very fast normal freezes are made by extracting the pyrex tubes holding 
the crucibles of molten tin from the furnace and allowing solidification to take 
place in room air. Slow normal freezes are obtained when undisturbed freezing 
takes place in a slowly cooling block furnace; cooling curves from such freezes 
are characterized by well-rounded contours, a consequence of the large amount 
of supercooling found on freezing high purity tin, which rule out their use for 
precision-fixed temperature points. The bottom curve in Fig. 1 illustrates the 
type of cooling curve that is obtained by nucleating (technique described in 
Part II, 19575) the melt in the throat of the furnace and then immediately 
lowering it into a slowly cooling furnace block; this outside nucleated slow 
freeze results in a plateau of nearly constant temperature on the cooling curve 
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which serves as an excellent index for a temperature standard and, as an 
additionally desirable feature, this type of freezing has very short recalescent 
periods. 

In this survey, estimates of the melting ranges of the samples were deter- 
mined from heating curves after very fast freezes and after outside nucleated 
slow freezes, and in addition the ranges of freezing temperatures were estimated 
for all samples from the cooling curves obtained on outside nucleated slow 
freezes; only freezes, and melts after freezes, that had a supercooling of > 4° C 
are used in this survey. Table III gives the approximate durations and trans- 











TABLE III 


Transformation conditions for the thermal survey 








Freeze 
Transformation Type Designation Figure duration 














Freeze Very fast normal = 1 10 min 
Freeze Outside nucleated a 1, 2(a), 2 hr 
slow 3 
Prior freeze Solid Melting 
duration anneal duration 
Melt After outside b 2(b;), 4 2 hr 2 hr 2 hr 
nucleated slow 
freeze 
Melt After outside bo 2(b2), 5 2hr 17 hr 2hr 
nucleated slow 
freeze with 
overnight anneal 
Melt After very fast a 2(c), 6 10 min 2hr 2 hr 
normal freeze 
Melt After very fast Ce 2(c2), 7 10 min 17 hr 2 hr 






normal freeze with 
overnight anneal 












formation conditions for several series of cooling and heating curves from 
which estimates of the alloy ranges of transformation temperatures were 
made for the 10 samples of high purity tin; samples were annealed at 215+5° C. 

Resistance balances were estimated to the nearest micro-ohm and times 
were recorded to the nearest second. From direct time-temperature (in ratio 
units) plots of these measurements the percentage of the total transition time 
spent in equal temperature intervals was derived* and the information was 
then used in two ways to make effective comparisons of the alloy transition 
ranges of the 10 samples. 

First, histograms of the percentage time vs. temperature were plotted from 
which the temperature widths at the 1% level of the total transition time 
were estimated for all samples. Figure 2 shows a set of these histograms that 
were obtained for sample Sn.9 from Vulcan Detinning Spectrographic grade 
tin. Since this material was found to be the purestf tin that was examined in 


*This operation approximates the temperature derivative of the direct time-temperature 
curve; a plot of these time increments against temperature is referred to as the derived inverse 
(or inverse-rate) transition curve (see Hume-Rothery et al. 1952). 

+Purest sample is used in the thermal sense; narrow alloy ranges are taken to show that the 
associated samples contain dissolved impurities in amount and kind that place the least strain 
on the solvent lattice and hence their solid solutions most closely approximate an ideal pure tin. 
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Fic. 2. Inverse transition curves derived for Sn.9: a, outside nucleated slow freeze; b;, bz 
melts after slow freezing; ci, cz melts after fast freezing; subscripts 1, 2 refer to melting following 
2-hour or overnight solid anneals respectively. 
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this investigation, these inverse curves were derived from percentage time in 
temperature intervals of ~0.00013° C (i.e. 50X10-° ratio units) and hence 
show alloy transition ranges of the purest tin under the highest temperature 
resolution that was used in this survey. In Fig. 2 temperatures are plotted 
directly from arbitrary ratio units; the plateau temperature on the outside 
nucleated slow freeze was not selected as the arbitrary temperature zero on 
these plots because temperature differences from histogram to histogram may 
be uncertain to 0.0005° C due to the shortened temperature measurement 
plan and the approximate differentiation operation used in the survey. 

For the wider examination of all the samples durations in temperature 
intervals of ~0.00027° C (i.e. 100 X10-* ratio units) were used in plotting the 
histograms from which the values of the widths at the 1% level of the total 
transformation time listed in Table IV were obtained. The estimates (Table IV) 


TABLE IV 
Width at the 1% level of the total transition time (°C) 


Sample a by bs a Ce 
Sn.9 0.0013 0.0023 0.0022 0.0015 0.0015 
Sn.8 0.0013 0.0023 0.0016 0.0015 0.0015 
Sn.13 0.0013 0.0030 0.0011 0.0047 0.0026 
Sn.6 0.003 0.004 0.004 0.003 0.003 
Sn.3 0.003 0.008 0.005 0.004 0.002 
Sn.7 0.003 0.006 0.006 0.005 0.006 
Sn.2 0.004 0.008 0.012 0.007 0.008 
Sn.11 0.006 0.009 0.014 0.006 0.006 
Sn.12 0.006 0.011 0.019 0.006 0.005 


Sn.10 0.008 ~0.02 ~0.03 ~0.01 ~0.02 


Note: Alloy transformation ranges: a, outside nucleated slow freeze; bi, b2 melts after slow 
freezing; ¢1, ¢2 melts after fast freezing; subscripts 1, 2 refer to transformations following 
2-hour and overnight solid anneals respectively. 


of the widths at the 1% level, with the exception of samples Sn.9, Sn.8, and 
Sn.13, are given only to the nearest millidegree, since variations of this amount 
were found for this melting range parameter on repeated transformations of 
the same type on the samples of lower purity tin. 

The temperature widths at the 1% level of the total melting time from 
similar histograms have been used previously as useful indexes for comparing 
the alloy transformation ranges of temperature among several samples of zinc 
and cadmium but they must be used with some reservations with tin. The 
large and variable amount of supercooling found with high purity tin makes it 
very difficult to arrange nearly identical freezing environments for different tin 
samples. Although the furnaces may be adjusted to give the same rates of 
cooling on the samples prior to extraction into room air or into the furnace 
throat, depending on the type of freeze, it was not possible to have identical 
initial crystallization conditions exist in the different melts on the same types 
of freezes, or, indeed, in a single sample on repeated freezes of the same kind. 
A consequence of this failure to duplicate the thermal environment at the 
moment of primary crystallization results in a diversity of growth conditions 
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which causes variations in the initial segregations of soluble impurities even 
on the same type of freeze. These variable segregations on freezing in turn give 
rise to variable secondary arrests on melting that last well over 1% of the total 
melting time and make the 1% level melting range parameter much too 
sensitive for comparing other than very pure tin samples (i.e. samples with 
alloy ranges of <0.002° C). 

For these reasons the derived information ‘‘% duration vs. temperature inter- 
val”’ obtained from direct time-temperature curves was used in an alternative 
manner to make reliable comparisons of the transition ranges of all the samples. 
The percentage of duration was incremented* to 95% of the total transition 
time beginning at and incrementing in both directions from the temperature 
interval having the maximum percentage of duration (i.e. the temperature 


100 
80 
60 


40 


OF TOTAL FREEZING TIME 


to 


20 





° 
09,0001 0.001 0.01 
FREEZING RANGE (°C) 


Fic. 3. Fraction of the total freezing time spent in given freezing ranges on outside nucleated 
slow freezes: derived from ‘a’ curves. 


range in which the greatest amount of the sample transformed). The summa- 
tions were cut off at 95% of the total transition time for two reasons: the 
estimations of the total transition times may be in error by a few per cent and 
also the initial and final few per cent of the total time of a melting or freezing 
curve may reflect curvatures, caused by thermal lag in the furnace block, that 
might not be quite identical from melt to melt because the final adjustments 
on furnace powers are made only 10 or 15 minutes before the arrest. 

The total transition time is arbitrarily estimated by selecting the duration 
that appears to encompass the complete transformation on the direct time- 
temperature plot. With the exception of samples Sn.9 and Sn.8, difficulty was 
sometimes experienced in accurately identifying the time for the liquidus break 

*This operation approximates an integration of the inverse transition curve between specified 
limits; by summing in both directions from the intervals having maximum durations on all 


samples, melting range comparisons are based on the heart of the melting range information 
and variable end effects are minimized. 
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especially on 6; type (see Table III) melting curves; the temperature would 
quickly break upward for a few tenths of a degree from the liquidus point but 
during the next few minutes it would intermittently rise, drop, or remain 
nearly steady before positively following an all liquid heating curve. 

This erratic behavior is likely caused by the movement of remaining small 
amounts of solid still dispersed in the liquid melt. It is probably significant 
that this phenomenon is most pronounced on melting shortly after the slow 
freezing which appears to cause the greatest amount of segregation of the 
dissolved solutes in the ingot. 

Figures 3-7 inclusive show the results of plotting the ‘““% duration” sum- 
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Fic. 4. Fraction of the total melting time spent in given melting ranges following outside 
nucleated slow freezes with 2-hour solid anneals: derived from ‘b,’ curves. 

Fic. 5. Fraction of the total melting time spent in given melting ranges following outside 
nucleated slow freezes with 17-hour solid anneals derived from ‘b,’ curves. 
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mated in the above manner against transition range for the several series of 
freezes and melts on the 10 samples of high purity tin; trend lines instead of 
proper histograms are drawn on these plots to reduce any confusion in following 
particular samples. 

Several observations may be made on the results shown in Figs. 3-7. 
Narrowest* alloy transition ranges for all the samples were found on the outside 
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Fic. 6. Fraction of the total melting time spent in given melting ranges following very fast 
normal freezes with 2-hour solid anneals: derived from ‘ce,’ curves. 

Fic. 7. Fraction of the total melting time spent in given melting ranges following very fast 
normal freezes with 17-hour solid anneals: derived from ‘cs’ curves. 


*Narrow and wide alloy ranges are used here to denote temperature ranges over which the 
bulk of the individual samples transformed; undoubtedly there are regions in each ingot (e.g. 
near grain boundaries and the small volume immediately adjacent to the thermometer well) 
where segregated solutes (see Weinberg 1959; Kishkin and Bokstein 1955) become so concen- 
trated that these regions have solidus temperatures that are very much lower than the tem- 
perature at which general melting of the bulk of the particular sample begins. 
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nucleated slow freezes (Fig. 3); widest alloy ranges were found on remelting 
after these freezes allowing only short anneals near the solidus temperature 
(Fig. 4). By inspection of these curves it is easy to grade the samples on the 
basis of the width of their alloy transition ranges; the fact that their relative 
positions in an array based on increasing alloy range vary following different 
types of freezing and amounts of annealing is probably the result of variations 
in the distributions of different impurity solutes between grain volume and 
grain boundary regions. In order to illustrate these differences and to obtain 
numbers for specific comparisons among the samples, the alloy ranges corre- 
sponding to 50% of the total transition times were obtained from Figs. 3-7; 
these estimations are given in Table V, the final column of which anticipates 
results from the intercomparison of plateau freezing temperatures. 
TABLE V 
Temperature range (°C) for 50% duration 
Liquidus point 


rel. to Sn.9, 


a bh be ‘ 2 





0.0001 0.0003 0.0001 0.0 
(—7) (—20) (—7) 
0.0001 0.0003 0.0002 0.0 
(—17) (—17) (—17) 
0.0008 0.0003 .0004 
(—14) (—7) ) 
0.0012 0.0004 0.0005 0009 
(—12) (—10) (—10) 
0.0020 0.0008 0.0004 0.0002 .0012 
(—4) (—6) (—6) (—13) 
Sn.7 .000: 0.0022 0.0015 0.0015 0.0009 0010 
(—6) (28) (—8) (—1}) 
Sn.2 0.0008 0.004 0.005 0.004 0.003 —0.0090 
(—7) (—7) (—5) (—6) (—9) 
Sn.11 0.002 0.007 0.007 0.004 0.002 —0.0025 
(—6) (-8) (—6) (-—9) (-8) 
Sn.12 0.002 0.008 0.009 0.003 0.0015 —0.0048 
(—5) (—11) (—5) (—8) (—5) 
Sa.10 0.005 0.04 0.04 0.015 0.003 —0.0132 
(A) (—6) (—4) (—9) (—9) 


Note: Temperature ranges (°C) in which the samples spent 50% of their total transition times: 
a, outside nucleated slow freeze; b:, b2 melts after slow freezing; ci, c2 melts after very fast freezing; 
subscripts 1, 2 refer to transformations following 2-hour or overnight solid anneals respectively. The 
amounts (°C) of supercooling prior to solidification are given in parentheses beneath the respective 
transition ranges; for ready reference, a summary of the liquidus point intercomparison is given in 


the final column. 


The Vulcan Detinning Spectrographic tin (Sn.9 and Sn.8) has the narrowest 
alloy ranges of all the samples examined; further, the transition range of this 
metal is not greatly affected by the type of transformation or amount of solid 
anneal. Tables IV and V show that the samples which have liquidus points 
within 0.001° C of that of Sn.9 have alloy ranges which are both wider and 
more sensitive to transformation conditions than are transition ranges of 
Vulcan’s zone-refined material. The relatively poor showing of the Japanese* 


*We have examined a sample of nominal 99.9999% zone-refined lead from this source and 
found it to compare quite favorably with lead (not zone-refined) of this purity from other 
refiners. 
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zone-refined tin (Sn.12 and Sn.10) would suggest that these samples contain 
dissolved solutes that are not efficiently carried out with a moving zone. From 
this investigation it appears that for gauging alloy transformation ranges for 
tin in the purity range 99.99% to 99.9999%, the estimates (Table V) of the 
alloy ranges corresponding to 50% of the total transition times are probably 
more useful than the temperature widths (Table IV) at the 1% level from 
inverse curves because the latter parameter is too sensitive to minor variations 
in segregations of dissolved solutes on solidification following the large and 
variable supercooling. However, the 1% level parameter has the necessary 
sensitivity to detect significant differences in the alloy ranges of samples of 
the highest purity (i.e. samples with alloy ranges of <0.002° C). 

(a) Supercooling 

During the course of these thermal analyses some knowledge was gained of 
factors affecting the supercooling of high purity bulk tin melts (mass, 1.1 kg; 
size, 3.5 diameter X19 cm) contained in machined graphite crucibles. 

A high degree of supercooling is particularly important in the attainment of 
the plateau temperatures on outside nucleated slow freezes; it is essential that, 
when a melt strikes, the rates of nucleation and grain growth are sufficiently 
great to provide the release of enough latent heat of fusion to immediately 
raise the sample temperature to its liquidus point. For the 10 samples used in 
our experimental arrangement a supercooling of 4° C or greater provided large 
enough rates of nucleation and growth to assure attainment of the plateau 


temperature} for Vulcan Spectrographic grade tin the plateau temperatures 
could be obtained with only 2° C supercooling. 
The degree of supercooling attained before melts strike appears to depend 


onacombination of many variables of which the following general observations 


were made: 

(i) Least pure samples, Sn(2, 11, 12, 10), tended to supercool <10° C. 

(ii) It is best to handle high purity tin under an inert atmosphere and to 
omit topping the melts with the usual graphite powder, as an oxide inhibitor. 
Inadequate (<1° C) supercooling* was experienced, even with Vulcan Detin- 
nings’ Spectrographic grade tin, when graphite powder (200 mesh) became 
distributed throughout a sample; on loading Sn.8 for the second time some of 
this powder from the first loading became dispersed through the liquid metal 
and gave rise to difficulties in attaining the desired 4° C supercool on outside 
nucleated slow freezes. Thereafter graphite powder topping was omitted when 
the samples were loaded. 

A sample that has been contaminated with graphite powder in the above 
manner may be cleaned by heating the sample ~50° C above its liquidus 
point and allowing time for the powder to float to the surface where it may be 
spooned off with a pyrex ladle; cycling the sample through a few freezes and 
melts also appears to aid in clearing out the external nucleating agents. 

(iii) The extent of the supercool appeared to depend on the amount and 

*These experiences on bulk samples appear to be contrary to that reported by Pound and 


LaMer (1952); they found that powdered graphite did not catalyze the nucleation rate in 
finely divided tin droplets (diameters, 1-30 microns). 
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duration of the previous overheat; melts that were held 10° C above their 
liquidus points overnight nearly always supercooled more than the same melts 
when refrozen <2 hours after melting. Further, when melts were held <2° C 
above the liquidus points for <2 hours and then refrozen the observed super- 
cooling was usually <4° C. These trends became accentuated when using the 
least pure samples, Sn.2, 10, 11, 12, to such an extent that we conclude that 
the supercooling of samples having alloy melting ranges of, say >0.01°C 
definitely exhibits a dependence on the thermal history of the liquid sample. 
In many respects, these observations are similar to the effect of thermal 
history on the supercooling phenomena on tin observed by Webster (1933), 
the existence of which, Turnbull (1950, 1953) attributes to the presence of 
crystal nucleation catalyst particles in the liquid. Webster does not give either 
the mass or the purity of the bulk tin used in his experiments: his tin was 
contained in quartz crucibles open to the air. 

(iv) On many freezes the rate of cooling was observed to slow up and indeed 
even stop and reverse slightly, in the temperature range up to 4° C undercool, 
before increasing rapidly to a final momentary arrest before general nucleation. 
At this temporary arrest* in the region up to 4° C supercool significant nuclea- 
tion must be occurring to release enough latent heat to stop the cooling of the 
1.1-kg melts that was observed. When the samples were in this rather indecisive 
but sensitive state, general nucleation could apparently be triggered by a loud 
noise near the furnace. In contrast, when the melts had safely cooled below 
this sensitive region vibration and noise appeared to have no further effect on 
nucleation; nor did 10 r/min y-rays from a Co® beam therapy unit. Whether 
the nucleating agents present in the 0—4° C range are different from the ones 
active in the general nucleation is not known. 


OCCURRENCE 
@ 


So 


0 5 10 15 20 25 
AMOUNT OF SUPERCOOLING (°C) 


Fic. 8. Histogram of the amount of supercooling attained before general nucleation on 
freezing all samples. 


*These observations refer only to supercooling in freezes when the melts are extracted into 
the furnace throat or into room air (i.e. relatively rapid supercooling >1° C/min). When this 
type of nucleation occurs on a normal slow freeze (i.e. rate of supercooling <0.1° C/min) the 
melt temperature has been observed to rise and fall in this 4° C range several times before 

eneral nucleation begins; similar indecision in the recalescent state was noted on normalfslow 
uous in high purity zinc (see Fig. 5, McLaren 19575). 
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Webster (1933) found that his samples supercooled either 3 to 4° C or 20 to 
25° C and he reported only one observation between these two ranges; the 
information given in Fig. 8 includes supercooling data obtained for all samples 
using the freezes that were made in both the thermal survey and in the. pre- 
cision intercomparison of the liquidus points. 

(b) Anomalous Structures on the Front End of Melting Curves 

In our earlier work* examples were presented of anomalous structures that 
were found on the front ends of many melting curves on high purity tin. These 
were sudden discontinuities on the melting contours in which the melt tempera- 
ture abruptly rises a few millidegrees and may or may not then decrease 
0.002° C before completing the transformation with either a long alloy slope of 
steadily increasing temperature or a nearly constant plateau-like feature, pre- 
sumably dependent on the concentration ranges of the associated dissolved 
solute in the solid that is being melted. Figure 9 shows a selection of similar 
front end structures (usually found during the first quarter of the melting 
time) that were observed in this investigation. Neither the prominence nor the 
form of the structure appeared to correlate with the amount of supercool or 
the type of the prior freeze; on the other hand, perhaps the most interesting 
structures were found on melting curves made on ingots following overnight 
solid anneals. 

The most significant observation on these structures made in the present 
investigation was that no marked evidence of front end structures was found 
on the melting curves on Sn.9 and Sn.8. Since these samples also had the 
narrowest alloy melting ranges we conclude that the front end structures are 
probably associated with impurity segregationsf in the ingots. The lack of a 
correlation of these effects with prior supercool and type of freeze seems to 
militate against a strain anneal explanation whereas impurity diffusion might 
account for the greater frequency of these structures on melting curves follow- 
ing overnight anneals. We speculate that a possible cause of the anomalous 
superheating that is observed might be a decrease in internal stress in the 
solid following a rapid diffusion of impurity atoms at the liquidus temperature; 
further, because these effects were not observed on samples Sn.9 and Sn.8 we 
consider that instrumental causes similar to those discussed by Smit (1957) are 
unlikely in our experimental arrangement. The latter observation also appears 
to rule out the possibility of changes in internal stress (and hence variations 
in the melting temperatures) with a release of trapped gas since all samples 
were solidified under dry nitrogen. 

A final observation was that the maximum temperature attained on a front 
end anomaly never exceeded the temperature of the liquidus break on the 
melting curve; this factor clearly differentiates this type of apparent over- 
heating from the superheating of 1 to 2° C found by Haykin and Bénet (1939) 
on melting large single crystals of high purity tin. After completion of the 

*Part II (19576): In Fig. 13 of this reference the curves labelled (d) and (b) should have read 
(6) and (d) respectively. 

tDoping Vulcan's Spectrographic grade tin with radioactive solutes and correlating auto- 
radiographs of solidified ingots with their corresponding thermal analyses should test this 


conclusion and provide information on the thermal effects and distributions of specific 
impurities. 
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Fic. 9. Portions of melting curves showing anomalous front end structures. Types of 
melting curves: Sn.9(c2), Sn.13(d1), Sn.7(c2), Sn.6(b2), Sn.13(b2). 
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thermal survey a melt (Vulcan Extra Pure from Lot VS-65) was decanted 
30 minutes after striking on a slow outside nucleated freeze; prior supercool 
was 11° C and the furnace was adjusted for a ~2-hour freeze. Figure 10 shows 
a photograph of the sectioned shell that remained in the crucible on decanting; 
a smooth interface, containing many dendrites, appears to have advanced in a 
nearly uniform fashion (grain size, irregular 5-20 mm) from the walls of the 
crucible. The dendrites* that project ahead of the smooth interface undoubtedly 
grew very quickly in the minute or so after general nucleation as the melt 
recovered from the large supercool (see Weinberg and Chalmers 1951, 1952). 
Cylindrical shells have previously been found on decanting zinc melts that 
were partly solidified in the block furnace; they provide an ideal environment 
for the sensing coil of the resistance thermometer by both enclosing and 
shielding it with the freezing interface. 


(B) Intercomparison of Liquidus Points 

A careful intercomparison of the plateau freezing temperatures attained on 
outside nucleated slow (~3 hours) freezes on the 10 samples of tin was made; 
the measuring sequence was similar to that used in the intercomparison of 
samples of high purity zinc (Part III, 1958a). For such measurements it is 
normal to compare the liquidus points of two samples at nearly the same time 
with a single thermometer in order to minimize the effects of any short-term 
systematic variations in bridge, thermometer, or water triple cell temperature 
during the comparison. The fully corrected zero current resistances at the 
liquidus points of a single Meyers thermometer, S155, at 15cm mid-point 
immersion corrected to 1 atmosphere dry nitrogen pressure over the melts 
were determined; corresponding resistances at the triple point of water (triple 
cell 29) were used to baseline all measurements. Figure 11 gives the results of 
comparing samples Sn.8 and Sn.13 against Sn.9, eight determinations being 
made on each sample; Table VI summarizes the intercomparisons that were 
made on all the samples. 

Sn.9 and Sn.8, Vulcan-Detinning’s zone-refined tin, had the narrowest alloy 
transition ranges and for this reason are considered to be the samples which 
have purities nearest to an ideal ‘pure’ tin; and since the bulk of these samples 
transformed in <0.002° C under a variety of transition conditions (Figs. 3-7) 
it is very probable that their liquidus points are <0.0005° C from that of pure 
tin. Further, since a value of 231.9119° C (Int. 1948) for the liquidus point 
of Sn.3 under a total pressure of 1 atmosphere was determined previously 

*We do not discount the melting of the material in these dendrites as a possible cause of the 
anomalous structures observed on the front end of melting curves; since they grow during a 
period of very rapid freezing their average impurity concentration is thought to be considerably 
greater than that of material that solidifies at the plateau freezing temperatures. 

+Plateau freezing temperatures are used synonymously with the liquidus points of the alloys; 
plateau freezing temperatures are invariably closely associated with the temperature of the 
liquidus breaks on melting curves. It is thought that, on the appropriate phase diagrams for 
these samples, the liquidus for the range of concentrations that solidify at the plateau freezing 
temperature is so close to horizontal that the resistance thermometer is unable to detect its 
slope; and in addition, the material freezing out at the plateau freezing temperature has a nearly 


constant composition which closely approximates the equilibrium concentration in the solid at 
the liquidus temperature. 
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Fic. 11. Intercomparison of plateau temperatures of Sn.8 and Sn.13 with Sn.9; Rrp and 
Rgn, zero current fully corrected resistances of the thermometer at the triple point of water and 
the tin point respectively. 


TABLE VI 
Intercomparison of plateau temperatures 





Ratio Rga /Rrp Reproducibility 
——_—_—_—————————_ (equivalent 10-*° C) Plateau 
Mean No. of 1ex =——————_ temp. rel. to 
date Sample freezes Mean std.dev. Std.dev. Spread Sn.9, °C 


.89252821 +51 +1.4 —0.00123 
3605 +20 +0.5 


.89252877 +33 +0.9 
3593 =+37 +1.0 


.89253525 +54 +1.5 
3205 +33 +0.9 


.89253227 +24 +0.6 
3220 +31 +0.8 


.89249855 +20 +0.5 
53188 +24 +0.6 


.89248301 +215 +5.8 
53212 +29 +0.8 


.89252282 +10 +0.3 
3196 +17 +0.5 


.89251428 444 +1.2 
3206 +25 +0.7 
.89253034 +35 +0.9 
3189 +22 +0.6 


29.4.59 Sn3 .89252819 +46 +1.2 
*Sn.9 3262 +30 +0.8 


*The average Rgn/Rrp for the 56 freezes on Vulcan Detinnings’ Spectrographic grade tin (Sn.9, Sn.8) was 
1.89253215 with a standard deviation equivalent to +0.00009;° C and a spread of 0.00044° C; this reproducibility 
is comparable with that of any fixed temperature point on the International Temperature Scale. 
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(McLaren 1959) from measurements with 11 standard thermometers, it is 
deduced, from the additional information given in Table VI, that the standard 
liquidus points of samples Sn.9 and Sn.8 have a value 231.9131° C (Int. 1948). 


DISCUSSION 


The selection of samples of highest purity metal on the basis of narrowest 
alloy transition ranges was used to establish precision thermometric fixed 
points at the standard liquidus points of zinc, cadmium, and indium; in the 
particular samples that were examined in the zinc and cadmium investigations 
the ones with the narrowest alloy melting ranges also happened to have the 
highest liquidus temperatures. An interesting difference has been found in this 
study on high purity tin: Sn.6 was found to have a liquidus point that is 
0.0009° C higher than that of the metal with the narrowest transition range. 
The explanation of this difference undoubtedly lies in the nature of the specific 
impurities contained in the respective samples. 

On solution, solutes which are more soluble in liquid than solid alloy lower 
the liquidus points of their alloys from that of the pure solvent; conversely 
solutes more soluble in solid than liquid alloy raise liquidus temperatures on 
solution (for a discussion on these points see Pfann (1958)*). Therefore to 
explain the value of the liquidus point that was found for sample Sn.6 we look 
for a solute that raises the liquidus temperature on solution in tin; antimony 
appears to be the most likely impurity on two counts. Firstly, Vulcan Detin- 
ning’st Extra Pure grade tin is prepared by further electrolytic refining of 
selected samples of Vulcan 99.99% Electrolytic tin which contains about 
20 p.p.m. antimony. The Vulcan Spectrographic grade tin is prepared by zone- 
refining selected lots of the Extra Pure grade tin which further effects a partial 
removal of any remaining antimony and substantially reduces the level of 
other impurity solutes. On analysis the limit of detection of antimony in the 
final products is thought to be about 5 p.p.m. Secondly, since a reasonable 
suspicion of antimony in Sn.6 has been established it is also significant that 
6 p.p.m. antimony in tin would raise the liquidus point by the observed 
0.0009° C (phase diagram, Metals Handbook 1948). 


CONCLUSIONS 


This investigation on 10 samples of the highest purity tin that could be 
obtained has provided the broad evidence required to confirm the results of 
the earlier investigation, namely, that the liquidus point of high purity tin is a 
thoroughly reliable thermometric fixed point of the highest reproducibility and 
experimentally may be simply and quickly realized in any well-equipped 
laboratory. This study has shown that alloy melting range parameters are 
sufficiently sensitive to enable the purest samples of tin to be easily selected. 

*Pfann (1958) also gives references to the purification of tin by zone refining; Aleksandrov 
et al. (1956) give a thorough account of the purification of tin by zone-melting techniques and 
in addition ny on the effectiveness of continuous heating under high vacuum for the 


elimination of highly volatile impurities. 
{We are indebted to Robert H. Taylor of the Vulcan Detinning Co., Sewaren, N.J., U.S.A., 


for this information. 
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A value of 231.913° C (Int. 1948) was found for the standard liquidus point 
of pure tin; samples that melt over a range of <0.002° C will have liquidus 
points within 0.001° C of this value. 

The three fixed points, triple point of water and the standard liquidus points 
of high purity tin and zinc, provide a ready means for any industrial, univer- 
sity, or institutional laboratory to quickly and reliably carry out high precision 
calibrations (for interpolation considerations see McLaren (1959)) of resistance 
thermometers for subsequent use in the range 0 to 630° C (Int.). 
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APPENDIX 


During March 1959 an intercomparison of the temperatures realized in six 
water triple point cells was carried out with the same measuring scheme that 
was used in earlier intercomparisons (Part I, 1957a). This laboratory has 
retained four of the National Research Council cells, namely triple cells 29, 
24, 21, 31, that were included in the former surveys and recently has added 
two more cells, 203 and 204, prepared by the J. and J. Instrument Co., Silver 
Spring, Maryland, U.S.A.; this group of cells now constitutes the grid for 
assessing our realization of the triple point of water and for gauging changes 
in particular cells. Table VII gives the results of the present intercomparison; 
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TABLE VII 


AR‘p of cells relative to T.C. 29. The limiting uncertainty in a 
difference measurement is of the order of 3 wQ : 10 uQ=0.0001° C 


March 1959, pQ 


Nov. Dec. /1954, Apr. May/1956, 
pQ 


Ist week 2nd week 


0 
+3 
“Pd 
—16 


each entry was determined from the average of 15 independent determinations 
of the annealed* zero current resistance at the triple point of water, Rrp, made 
on each cell during a 4-day interval. Table VII shows that triple cells 29, 24, 
and 21 have maintained their relative triple point temperatures during the 
past 3 years while that of triple cell 31, already considered substandard, has 
decreased another 0.0002° C. As in the earlier investigations no systematic 
measurable (i.e. >3u2) change of cell temperature with age of the ice mantles 
was found from the second to the fifth day. 


*All thermometers were annealed to stability at 450° C prior to each 4-day intercomparison. 









THE RADIATION CHARACTERISTICS OF A SINUATE ANTENNA! 


S. C. Lox? 








ABSTRACT 


Rigorous expressions for the radiation field of an antenna of sinusoidal shape 
are derived on the assumption of a travelling-wave type of current distribution 
along the conductor. In order to illustrate the derived results an endfire antenna 
is designed and some calculations of radiation field are presented. 







INTRODUCTION 


The zigzag antenna described first by Cumming (1955) is characterized by 
two significant properties: namely, the moderate bandwidth in the V.H.F. 
and U.H.F. ranges, and the strong axial beam of radiation with very low side 
lobes. In order to derive the far-field expressions, Cumming treated the 
structure as an array of V-shaped elements. Assuming a travelling-wave type 
of current distribution along the conductor the two electric field components 
are determined for one turn. Considering the single turn as one element in an 
array, the antenna radiation pattern is obtained by multiplying the unit 
pattern by an array factor appropriate to the number of turns and the turn 
spacing of the antenna. The radiation properties of a single zigzag antenna 
were further investigated and reported in a recent article by Sengupta (1958). 

A different and somewhat more rigorous method of solving the radiation 
problem of the zigzag antenna has been presented by Wong and Loh (1959). 
Assuming a travelling-wave type current distribution, the field components 
are formulated in terms of a vector potential. The total vector potential, and 
consequently the radiation field of the antenna, are determined by integrating 
the appropriate integrals over the length of the antenna. 

In the present paper, the radiation problem of an antenna of sinusoidal 
shape (see Fig. 1), which is very similar to the zigzag antenna mentioned above, 
is investigated. Expressions for the radiation field are derived by using the 
method similar to the one employed by Wong and the author. In order to 
illustrate the derived results, an endfire antenna of six turns is designed for a 
frequency of 2650 Mc/s and some calculations of the radiation pattern are 
presented. The theoretical values calculated from the derived expressions agree 
quite well with the measured pattern. 






























DERIVATION OF THE RADIATION FIELD 


The study was made on an antenna having the configuration shown in 
Fig. 1. The co-ordinate system is chosen in such a way that the antenna lies 
entirely on the x-z plane, the axis of the antenna being along the 2 axis. 


1Manuscript received October 1, 1959. 

Contribution from the Radio and Electrical Engineering Division, National Research 
Council, Ottawa, Canada. 

Issued as N.R.C. No. 5486. 

2National Research Council Postdoctorate Fellow. Present address: Department of Engi- 
neering, University of Hong Kong, Hong Kong. 
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Fic. 1, The sinuate and the zigzag antennae. 


Neglecting the effect of the ground screen backing the antenna, it was assumed 
that the antenna would carry a current distribution of the travelling-wave 
type with constant amplitude, and that the reflections from corners and the 
unterminated end of the wire would have no effect on the radiation pattern. 

The problem will be formulated in terms of a Hertzian vector potential. In 
order to facilitate the analysis, the vector potential will be determined in 
rectangular components, and by applying the appropriate rectangular- 
spherical co-ordinate transformation equations, the field is expressed in terms 
of a single spherical component. In the following analysis, the rationalized 
M.K.S. system of units is employed, and the time variation, exp(iwt), is 


understood. 
We shall first write down the parametric equations for the antenna: 


x =acost 
(1) y=0 0<t< bo, 
z= (b/2n) t 


where a is the amplitude of the cosine curve and b the turn spacing. 
The element of arc length, ds, along the conductor is given by 


(2) ds = /(dx’+dz") 
= (a/m)\/(1—m’ cos’ t) dt 


where 
Bot US sed 
Vii+ @/42'07)] 


Thus the arc length s is 
a 
(3) s= a (m,t) 


where 


K(m,t) = fVaem cos’ t) dt 


is an elliptic integral of the second kind. 
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The function K(m,t) can easily be expressed in an approximate form given by 


(4) K (m,t) = c, t—ce2 sin 2t 


= 1-(848 +g Be m4. “a 
a ee ee ee 


2 
mm‘ 
Stet: 

To determine the radiation field of an arbitrarily shaped radiator, we shall 
employ an expression for the Hertzian vector potential. Let us examine the 


vector potential at point P, distant R from the origin of the co-ordinates. 


where 


C2 = 


Y 


Fic. 2. Sinuate antenna located at the origin of the co-ordinate system. 


With reference to Fig. 2, the Hertzian vector potential is 


(5) Il = pees fi I(s) exp(ik? +p) ds, 


inn Tel 
where 


« is the permittivity of free space, 

k = (2m/X), where X is the wavelength in free space, 

I(s) is the current in an element of length ds, 
located at a point denoted by the vector 7, 

Bp is a unit vector in the direction of R. 


Writing in terms of the co-ordinate system shown in Fig. 2, we can expand 
eq. (5) as follows: 


(6) II = eg ees fi I(s) exp[ik(x cos ¢ sin 6+ sin ¢ sin @ 


ims melw 


+z cos 6)] ds. 
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Let us assume the current distribution along the conductor is a travelling 
wave of constant amplitude and is given by 


(7) I(s) = Tuexo( —é#s) 


where J» is the amplitude and p is the phase velocity of the current along the 
conductor relative to the velocity of light in free space. 

Referring now to the rectangular components of the Hertzian vector 
potential II, it is seen that the expression for II, is given by 


a —alo x(k) fe 


4 relw 


so(-i8 = s) exp[ik(x cos ¢ sin 6+2 cos 6)] 
X sin ¢ dt. 


(8) IL 


In terms of the parametric equations (1), one can obtain upon substitution 
in eq. (8) the following result: 


ome ay aia to 
(9) See =—ale exp(—tkR) J sin Hhexp if “ K(m,t)+a cost cos ¢ sin 0 
0 


+5. t cos alta 


4 Trelw 


Simplifying eq. (9), we get 
—aly exp(—tkR) 


4 relw 


(10) Il, = 


where we introduce 


and 
= ka sin 6 cos @¢. 


In order to obtain a series expansion for eq. (10), we shall use the following 
well known Fourier—Bessel expansions, 


co 


(11) exp(ix sin 0) = pay J, (x) exp(iné), 


exp(ix sin 20) = Jn(x) exp (126). 


Utilizing the above expressions, eq. (10) becomes 
—alIy exp(—tkR) 
4 relw R 


The above integral can easily be evaluated by expressing sin ¢ in terms of its 
exponential equivalent; therefore 


(12) II, = » » J, (B) J. mf sin t exp[tat+72mt+ int]dt. 
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_alo e apt peat tel} 
8) B= 8 rein X x Jm(B) In(7) atn+2m+1 
exp[i(a-+n+2m— 1)to]— ty. 
atn+2m—1 

The z component of the Hertzian vector potential can be derived in exactly 
the same way. It is obvious that there is no y component. The expression for 
II, is given by 

_ =bIo exp(—ikR) pecellietest alte 

It is usually desirable to express the Hertzian vector potential in spherical 
co-ordinates. This can be easily achieved by utilizing the following trans- 
formation equations: 

II, 

(15) I] 6 
II, 







4- 










(14) Il, 







II, sin 8 cos @ + II, sin @ sin ¢ + II, cos 6, 
II, cos 8 cos @ + II, cos 6 sin @ — TI, sin 6, 
— II, sin @ + II, cos ¢. 














Also, by using the following vector relationship, the electric field components 
can be derived from the Hertzian vector potential. 


(16) E= w well ++¥ (V IT) 


where u is the permeability of free space. 

The expressions in eqs. (13) and (14) for the radiation field in the x and z 
directions are calculated for the sinuate antenna of six turns, for different 
values of p. It is found that the best-fitting pattern is obtained when p = 0.95. 










ENDFIRE OPERATION 
In order to illustrate the derived results, a sinuate antenna of endfire 
operation has been designed and tested. The condition which should be 
satisfied in order to have maximum radiation along the @ = 0° direction is as 
follows (see Fig. 2): 







—> 
art 


Nico 


; arc length AB— 





An experimental six-cycle array, backed by a circular ground plane, was 
constructed on the basis of the foregoing analysis. A design frequency of 
2650 Mc/s was chosen for convenience in pattern taking; the E- and H-plane 
patterns were measured. Both Cumming (1955) and Sengupta (1958) investi- 
gated only the radiation characteristics of a six-element array, and their 
results differ somewhat in beamwidth and bandwidth. Cumming found that 
with the zigzag antenna a bandwidth of 50% and a beamwidth of 35° can be 
readily achieved. However, Sengupta observed that a sharper beam is possible 
and found the half-power beamwidth of 22° in E plane and 28° in H plane. 

In the present paper the author has also studied the radiation characteristics 
of a six-cycle array. It is found that the measured E-plane pattern of the 
sinuate antenna shown in Fig. 3 agrees quite well with the calculated pattern 
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Fic. 3. E-plane patterns at 2650 Mc/s. 
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Fic. 7. Measured H-plane pattern at 2700 Mc/s. 


shown in the dotted curve. A number of radiation patterns for different 
frequencies are also presented in Figs. 4, 5, 6, and 7, and it can be easily seen 
that the half-power beamwidth in the £ plane only varies from 28° at 2660 


Mc/s to 32° at 2700 Mc/s. 


CONCLUSIONS 


The general radiation characteristics of a sinuate antenna are very similar 
to those of a zigzag, but the former has a narrower bandwidth. 

The analysis presented in this paper seems quite satisfactory when compared 
with the measured radiation pattern of a six-cycle sinuate antenna. However, 
it has been observed that the radiation characteristics of this type of antenna 
are more sensitive to the axial length of the radiator than those of a zigzag. The 
endfire properties disappear as the number of turns is increased to seven and 
eight, but it regains them at nine turns. A fully satisfactory explanation will 
only be possible when the mechanism of wave propagation along the conductor 


is thoroughly understood. 
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DIFFRACTION BY A SPHEROID! 


BERTRAM R. LEvy AND JOSEPH B. KELLER 


ABSTRACT 


The diffraction of a spherical scalar wave by a hard or soft spheroid is investi- 
gated theoretically. First the diffracted field is determined by the geometrical 
theory of diffraction. Then for comparison the corresponding boundary value 
problem is solved exactly in terms of a series of products of spheroidal functions. 
The series involves the “radial” eigenfunctions which correspond to appropriate 
complex eigenvalues. Asymptotic expansions are derived for these functions for 
large values of the variable and the parameter. When used in the series solution, 
these expansions yield the asymptotic form of the diffracted field for incident 
wavelengths small compared to the spheroid dimensions. This result coincides 
precisely with that given by the geometrical theory. This agreement provides 
another verification of that. theory. The expression for the field is used to 
calculate the backscattering and the field on the spheroid. The electromagnetic 
backscattering is finally computed with the aid of a theorem which relates it to 
the two scalar results. 


1. INTRODUCTION 


A spheroid is a solid of revolution generated by rotating an ellipse about 
one of its axes. We consider the problem of determining the field which results 
when a wave is incident upon a spheroid. The customary manner of dealing 
with such a problem involves formulating it as a boundary value problem 
for the reduced wave equation. This problem is then solved by means of a 
series of products of spheroidal wave functions. The resulting series has been 
extensively studied in the literature. A compilation of the properties of the 
spheroidal functions contained in this series together with a complete biblio- 
graphy of previous work on diffraction by a spheroid is contained in the 
book by Flammer (1957). However, this well-known solution is practically 
useless when the wavelength of the incident field is small compared to the 
dimensions of the spheroid. It is just this case which we treat here, and 
therefore we must use a different procedure. In fact, we use two different 
procedures, both of which yield the desired result. 

Our first result, given in Section 2, is based upon the geometrical theory 
of diffraction (Keller 1958). This is a theory which describes the field scattered 
by any object directly in terms of geometrical properties of the object and 
certain coefficients determined from a few canonical problems. We have 
previously applied this theory to diffraction by convex cylinders (Keller 1956) 
and convex solids (Levy and Keller 1959). Whenever possible we compared 
the results of that theory with the asymptotic form (for wavelength small 
compared to object dimensions) of the solution of the corresponding boundary 
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value problem. In every case the agreement was perfect, including that of 
the elliptic cylinder (Levy 1958). Although it has not yet been proved that 
this agreement will always occur, we are convinced that this is the case. 

In Section 3 we formulate and solve the associated boundary value problem. 
Then by deriving asymptotic expansions of the spheroidal functions the 
solution is asymptotically expanded. Our results closely parallel those of Ritt 
and Kazarinoff (1959), who evaluated the field on the surface of the spheroid. 
We find that for points off as well as on the spheroid the asymptotic expansion 
of the diffracted field agrees with geometric theory given in Section 2. How- 
ever, to recognize the identity of the two results for points not on the surface 
of the spheroid it is necessary to utilize certain addition formulas for elliptic 
integrals. 

As applications of our result, we compute the field scattered back to the 
source, and then find the backscattering cross section. We also evaluate the 
field on the surface of the spheroid. 

All of our results pertain to a prolate spheroid. The incident field is pro- 
duced by a point source on the axis of revolution. The field is a scalar which 
satisfies the reduced wave equation and the spheroid is either hard or soft. 
In the case of backscattering, the cross section for an electromagnetic wave 
hitting a perfectly conducting spheroid is also given. This calculation is based 
upon a theorem proved in Appendix IJ. It relates the electromagnetic back- 
scattering for axial incidence, on a body of revolution, to the corresponding 
scalar results for hard and soft bodies. 


2. GEOMETRIC CONSTRUCTION OF THE DIFFRACTED FIELD 


Let us consider a prolate spheroid in a homogeneous medium, having 
propagation speed c. Suppose that a point source of radiation, located on the 
axis of revolution of the spheroid, emits a wave of angular frequency w. Then 
k = w/c is the propagation constant of the medium for such a wave. When 
the incident spherical wave strikes the spheroid, reflected and diffracted 
fields are produced. The reflected field u, can be constructed by the method 
of geometrical optics, as described, for example, by Keller, Lewis, and Seckler 
(1956). The diffracted field ug can be constructed by means of the geometrical 
theory of diffraction, as has been done by Levy and Keller (1959) for diffrac- 
tion by any smooth convex scattering object. When their result is specialized 
to the case of a point source on the axis of a body of revolution it becomes 


__ ae ‘ rn .'" Dm(P1)Dm(Q1) 
(1) ua(P) tel 47Ro stiles) m - T 

1+exp( er f oir) 

0 


S oP J Tope ue 
* exp tht si) -J. andr} "Usere(se— po pa)) 


+ —— Pal E Pale) exp it (+s) £ aia 
1+exp( iar ff oar) 
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In (1) P denotes the point at which the field is evaluated. The plane con- 
taining P and the axis we will call the xy plane. Then s; and se denote the 
lengths of the tangents from P to the spheroid in this plane and P;, P, 
denote the corresponding points of tangency (see Fig. 1). The points Q; and 


Fic. 1. A cross section of the ellipsoid showing the source Q and a field point P. Two rays 
from Q are tangent to the ellipsoid at Q; and Q2 where they produce diffracted surface rays. 
The ray from Q; travels a distance t; along the ellipsoid to the point P; from which it sheds a 
diffracted ray which passes through P. In addition it continues to encircle the ellipsoid and 
sheds another ray through P each time it passes P;. The ray from Q2 sheds rays through P each 
time it passes through P2. These rays pass through the axial caustic at Pj before reaching P. 


Q»2 are the points of tangency of the tangents from the source point Q to the 
spheroid in this plane. The distances from the axis to P;, P2, and Q; or Qe 
are denoted by 71, r2, and 79 respectively. The total arc length of the generating 
ellipse is denoted by 7, while t; is the arc length along the ellipse from Q; to 
P, and ¢2 is the arc length from Q2 to P2. The length of the tangent from P; 
to the point Pat which it intersects the axis is denoted by pi, while p2 denotes 


te om 


Fic. 2. A cross section of the ellipsoid showing the points Q; and Q2 at which two incident 
rays are tangent to it. The inc\dent field is a plane wave coming from the right. The tangent 
rays produce diffracted rays wich travel distances 4; and ¢, to a point P on the surface. 


the distance from P2 to Pj. The integrals are with respect to arc length + 
along the generating ellipse. The decay coefficient a(7) is defined in terms 


of b(r), the radius of curvature of the ellipse, by 
(2) am = e~ */6(R/6)1/8g,,6-2/3, 


The constant gm is defined in Appendix I. For the soft spheroid, gn = 
while for the hard spheroid gm = g°. When it becomes necessary to dis- 
tinguish between these two cases the coefficient a, must be given the same 
superscript as gm. The diffraction coefficient D,, depends upon the radius of 
curvature at the point at which D,, is evaluated. For the soft spheroid 
D, = D® is given by . 


(1) 
m 
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(3) DY ee wl gu 6728 ye {A (qe) } "eS ye 
For the hard spheroid D, = D® is given by 


(4) D® = a! 9-4 67s pos gow A ey eft /24 pe 


The Airy function A is defined in Appendix I. 

The result (1) can be rewritten in another form by introducing prolate 
spheroidal co-ordinates £, 7, ¢. These co-ordinates are most easily defined by 
relating them to the Cartesian co-ordinates x, y, z. Let the axis of symmetry 
of the spheroid be the x axis and let the origin be at the center of the spheroid. 
Finally, let 2h denote the interfocal distance of the generating ellipse. Then 
the prolate spheroidal co-ordinates may be defined by 


x = hcosh é cos 7, 

y = hsinh £ sin cos ¢, 

z = hsinh £ sin nsin ¢. 
The angular variables 7 and ¢ range from 0 to mw and 0 to 2m respectively, 
while 0 < — < ~. The spheroid may be defined by the equation = a. The 
point source Q will be assumed to lie at 7 = 0, § = & > a. The n co-ordinates 


of the points of tangency Q:, Pi, Qe, and P» will be denoted by m, 2, 3, and 
na respectively. Their other co-ordinates are § =a and either ¢ =0 or 


The geometrical quantities in (1) can be expressed in terms of these co- 
ordinates. Most of them have previously been calculated by Levy (1958) and 
the remaining ones can be easily found. When these expressions are inserted 
into (1), the result for the field diffracted by the soft spheroid becomes 


(5) ug(P) = 2? 27-97 6-4 ef? R-MS A (cosh a sinh a)" 


X (sinh £ sinh & sin 7)~”” (cosh” &) —cosh’a)~ 2 {fa(a, m1, 02) 


2x 
—if,(a, —03, —04)} LA"Gg?)} 7 A-tenp seh , (cosh’a—cos® n)'”* dn 
Sy 


2x 
+ir\” (kh)"*(cosh a sinh a)*”” J (cosh? a—cos” 9)~”” ant | 
0 


In (5), {2 = 6-/8e%/3g™ and } = 2r—n3, while nf = 4 if the ¢ co-ordinate 
of Pz is 0 and nf = 2r—n, if the ¢ co-ordinate of P: is r. We assume that 
the ¢ co-ordinate of P is 0. The function f, is given by 


(6) fr = (cosh? a—cos” n2)** { (cosh £ cos 7—cosh a cos 2) 
+ (sinh £ sin »—sinh a sin n°} exp| ik tanh £o(cosh” )—cosh? a)” 
1/2 


+ikh{ (cosh & cos n»—cosh a cos 72)’+ (sinh £ sin 7—sinh a sin n2)"} 


"2 
+ikh f (cosh? a —cos” 9)? dn+irs” (kh)** (cosh a sinh a)*”* 
rT 


"2 
x f (cosh? a—cos? »)~*” an| , 
m1 
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For the hard spheroid (5) applies if we replace 


nm” by 7h” and {A'(gh”)}-? by 3{gi?47(gh?)} 


The preceding result (1) or (5) is not valid on the surface of the spheroid 
nor on the axis of symmetry because both of these places are caustics of the 
diffracted rays. However, it can easily be modified to be valid in these places 
by applying the appropriate caustic corrections (Levy and Keller 1959). To 
obtain the field on the surface, we must multiply the mth term in the first 
sum in (1) by the factor 


(7) Rm = 276 RMS st?.4 (g) exp(—iks,—ix/12)b-*(P;), 
multiply the mth term of the second sum in (1) by R» with s; replaced by 
Sg and P, replaced by Pe, and then let s; and sz 0. Upon doing this and 


noting that P,, P2, and P become identical, while 7; = rz and P; = Pe, we 
find that (1) becomes on the surface of the hard spheroid 


(8) ug(P) si eho (4 wRo)'2 is ake ade! sade’ Mihalis (OVA 2 ae Bak ft 


Dm(P)A (Qn . 
cone ee are {D.(0) exp( its, f str) 
) Qa 


m Y 
1texp( ier f Omar 
0 
r 
—1Dm(Q2) exp(itts— f ms ‘ 
Qe 


For later use it is convenient to specialize (8) to the case of plane wave 
incidence. If the incident plane wave is given by u, = e~* then the diffracted 
field is given by (8) with e"°(4rRo)—'! replaced by 1. The points Q, and Q2 
become the ends of the minor axis of the generating ellipse of the spheroid. 
Thus, when we write (8) in terms of spheroidal co-ordinates for plane wave 
incidence, we obtain 
1/2 


- m(cosh a) 
(9) ug(P) = (sin »)'(cosh” a—cos’ )'4 


cos 
x 7 ex {G(1/2, n)}—iex {G(n, 34 /2)} 
<1 [1+exp{G(0,21)}lam A(Qm) 





Here we have introduced the function G(a,8) which is defined by 


8 
(10) G(a,8) = ikh f (cosh? a—cos” n)"dn+i(kh) "7. (sinh a cosh a)*” 


8 
x f (cosh* a—cos’ 9)~/"dn. 


To calculate the backscattering cross section of the spheroid, we again 
specialize (5) to plane wave incidence with wine = e~*?. We then apply the 
method of Levy and Keller (1959) to correct uq on the axial caustic. When 
this is done we find that the leading term in the expansion of uq(P) for a 
point on the axis becomes, for large values of x, 
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3n/2 


(11) ug(P) = —be™*"* exp ir(h)"” (cosh a sinh a) ns (cosh? a 


/2 
—cos’ »)~"”” an a’ kh (sinh a cosh a)” (xq Aes 


In (11), r = 67'%e*/8g, The values of 5, g, and 4 or A’ are given below for 
the two cases: 


6=1 
(12) gq = 1.469354 | for the hard spheroid; 
A = 1.16680 | 
6 = gA*3-'(A’)-? 
(13) = 3.372134 | for the soft spheroid. 
A’ = —1.059053 


The total backscattered field u is the sum of the diffracted field ug(P) and 
the geometrically reflected field u,(P). To calculate u,(P) we use the methods 
of Keller et al. (1956). We then obtain for the leading term in u,(P) the 
result 

hcosh’ a 


(14) u,(P) = + 5 exp{ik(x—2h cosh a)}. 


Here, the positive sign holds for the hard spheroid and the negative sign for 
the soft spheroid. Upon combining (11) and (14), we obtain 


4? cosh’ a 
<3 inh a 


E _25n*(kh)"*(sinh a)*” 
gA*6'3(cosh a)? 





(15) u(P) = exp{ik(x—2h cosh a)} 


exp|2ikh cosh a+ix/3 


3n/2 


+ir(kh)'(cosh a sinh a)*”* f (cosh? a—cos? aan} | . 


«/2 


Here, the plus sign and (12) hold for the hard spheroid and the negative 
sign and (13) hold for the soft spheroid. 


3. ANALYTIC SOLUTION 


In terms of our previously defined spheroidal co-ordinates the reduced wave 
equation corresponding to a point source located at & = &, » = 0 is 


5(n)6(E—&o) 


2 2 2 es 
(1)  ugettyyt+ue coth +n, cot n+ (kh)” (cosh” ——cos’ n)u = Fshunh tein 


In (1) we have made use of the fact that u is independent of ¢. The boundary 
value problem to be solved is to find a solution u(é,n) of (1) which satisfies 
the Sommerfeld radiation condition at § = ©. On the surface of the spheroid 
u must satisfy either of the following boundary conditions: 
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(2) u(a,n) = 0 (acoustically soft spheroid) 
or 

(3) ur(a,n) = 0 (acoustically hard spheroid). 


If we replace the right-hand side of (1) by zero, the resulting equation 
has product solutions of the form ®(é)¥(n). The functions ® and ¥W satisfy 
the equations 


(4) = +coth ‘ (kh)?(b’—cosh’ £)@ = 0, 
ly dv 
(5) ag? toot 9 Gt (eh)*(b—cos" ot wt: 


In (4) and (5) 6 is an arbitrary separation constant. We now note that for 
certain complex values 6, of 6, (4) has outgoing solutions V (é), which satisfy 
either of the conditions V‘?(a) = 0 or V‘?’(a) = 0.* Consequently, we may 
try to solve the boundary value problem by assuming that uw has the form 


@ 


(6) u(t) = X Wa(n) Vn (E). 


If we substitute (6) into (1) and make use of the orthogonality of the V{”, 
we find that W,, satisfies the equation 





dW. 5(n) 
(7) “Or 3" +cot sé + (eh) (b;—cos’ n)W, = Cron - 
In (7) the constant C, has the value 
eco —1 
(8) C, = V2(¢)| aah J {VS Ce} sinh & a | 


The orthogonality condition used in obtaining (7) is 


(9) J Vi? (€) Vi (é) sinh & dE = Spm J ”¢V(e)}? sinh € de. 


This condition follows from (4) and the boundary conditions satisfied by the 
V, An evaluation similar to the one carried out by Levy (1958) shows that 
for the soft spheroid 


(10) J {Vn(€)}*sinh & dé = {2(kh)6,°}~'Va"(@) sinh a Via). 
For the hard spheroid we obtain instead 


(11) ‘i {Va"(€)}*sinh & dE = —{2(kh)"b,”}~"V,"(@) sinh a ALO. 


*When it becomes necessary to distinguish between the sets of eigenvalues we will use b(t) 
to denote eigenvalues corresponding to the boundary condition (2) and 5°2) to denote eigen- 
values corresponding to the boundary condition (3). 
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In order to solve (7) we let W{?(n) be a solution of the homogeneous 
equation obtained from (7) by setting C, = 0. We require W‘(n) to be 
regular at n= 0. A simple calculation then shows that W{(y) is an even 
function of » and that W‘)(xr—») is also a solution of the homogeneous 
equation. Since 7 = 0 is a singular point of (7) the subsequent analysis will 
be simplified by shifting the delta function to 7 = mo and, finally, letting 
no — 0. Then the conditions determining W,(n) are that it be regular at 
n = 0 and = 7, that it be continuous across 7 = m0, and that W/(n) have a 
jump of magnitude C,(sin 90)~! across 7 = mo. By applying these conditions 
we find that W,(n) is given by 


Waln) = —Cy{sin nown(no)} "Wa (10) Wa (2-7), mind a; 
(12) W,(n) = —C,{sin nown(no)}~*Wa(e—n0)Wa(n), 0<1< 
Here w,(n) is defined by 
(13) won(n) = Wa?(9—n) Wan) + Wa (0) Wa —2). 


We note that w, is the Wronskian. of two linearly independent solutions of 
the homogeneous equation. Therefore, a simple calculation shows that it has 
the value 


(14) won(n) = 2222) © oy (2/2) WP e/2). 


sin 
Upon using (14) in (12) and letting 79 > 0 we obtain 
(15)  Waln) = —C,{2Wy?(4/2)Wa?"(e/2)}"Wa?(0) Wa? (2-2). 


We can obtain another expression for W,(n) by using (13) directly in 
(12). Upon doing this and using the fact that W{(n) has a logarithmic 
singularity at 7 = 7,* we obtain 


(16) Wa(n) = —Cy{lim sin noWa?"(e—10)} "Wa? 4—n). 
790 
In Appendix I, equation (18), we have evaluated the limit appearing in (16). 


We use this result, together with (16), (10), and (8), in (6). Then, in the case 
of the soft spheroid, (6) becomes 


Mi OS by Vn’ (Eo) Vn (€)Wa (=n) 


~ 2sinha 77 a : 
cose (61—cos's)""an} Va" (a) a5 Ve? (a) 
0 





(17) un) = 


For the case of the hard spheroid (17) still applies if the right side is multi- 
plied by —1 and if V{?(a) and V”’(@) are interchanged. 

We shall now obtain the asymptotic form of u for large values of kh by 
using the asymptotic expansions of the spheroidal functions given in Appendix 
I. We find, in the case of the soft spheroid, that (17) has the asymptotic form 


*See Appendix I. 
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(18) u(E, 9) ~ — 0 727876480 R- ME" (sinh — sinh £ sin )”” 
X (coshasinha)** (cosh? £o— cosh? a)~"/*(cosh? —cosh’a)~ (cosh*a —cos” n)~ 


x es {A’ ay exp{iF(t) +4F (y)} SPL = tenet ae 


In (18) we have introduced F() and G(n) which are defined by 


1/4 


g 
(19) F() = kh - (cosh? ¢—cosh? a)/?—i(kh)“*z" (sinh a cosh a)? 


é 
x J (cosh? ¢—cosh’ a)~ "dé, 


ll 


on 
(20) G(n) = kh J (cosh? a—cos” 9) "dn+i(kh)?7<” (sinh a cosh a)?” 


” 
x { (cosh? a—cos’ )~/"dn. 
0 


In the case of the hard spheroid (18)-(20) apply if 7°” is replaced by 7° 
and {.1’(g‘)}~? is replaced by 3{q{?.4?(g)}-!. 

It isa = i matter to show that the result for u given in (5) of Section 2, 
obtained by the geometric theory of diffraction, agrees precisely with (18) 
obtained by asymptotically expanding the exact solution. The necessary 
calculations are contained in Levy’s work (1958). The identity of the ampli- 
tudes follow from equation (I-1) of Levy’s paper, while the identity of the 
phases follows from (54), (55), (56), and (57) of that paper. 

Let us now calculate the diffracted field at a point P(a,n) on the surface 
of a hafd spheroid due to an incident plane wave e~**. To do this, we write 
(17) for a hard spheroid, set £ = a, multiply by 47h cosh & exp(—7kh cos £0) 
and let & — © to obtain 


peitsOy 11/6 1/82/18 oo b, Wi.” (4-2) exp{ikh f(b,)} a 


9 rw 
(21) u(a,n) cea 2, 


f i. 2 2 \1/2 \ 9 yay 
cos} k : (b,-—cos 7) ‘“dn ab {V, (a)} 
Here 


. to 
(22) f (bn) = lim \ i és (cosh? ¢— 63) "dé—cosh nat. 


§0-00 osh *bn 


Using the methods of Appendix I we find 


(23) 2 via) ~ — 9 (eh) *2""(sinh a)~"(cosh a)"e-?**/"gqA (ga). 


ab 


The evaluation of f(5,) has been carried out in a slightly different notation 
by Ritt (1959) and yields the result 


(24) f(b) = — "cos! n) "dn. 








| 
| 
| 
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Using (23) and (24) in (21), together with the asymptotic expansion of W{” 
contained in Appendix I, we obtain 


(cosh a) © eta __ 5 ttatn)— i@(n) outa 
25) u(an) ~ — OS oon > 755: 
(25) 4(@.) ~ Ging)" (cosh a—coe’ a)” 2» (1 eg? (g,7) © 


Here G(n) is defined by (20). This result (25) is identical with equation (9) 
of Section 2, obtained by the geometric theory. 


4. CONCLUSION AND ELECTROMAGNETIC BACKSCATTERING 


We have determined the scalar field diffracted by a hard or soft spheroid 
of dimensions large compared to the incident wavelength. This field is given 
by (1) or (5) at points off the spheroid surface and off the axis. On the spheroid 
surface it is given by (8) for the hard spheroid and by (9) if the incident field 
is a plane wave. The backscattered field on the axis is given by (12) for an 
incident plane wave. All these results have been derived by using the geometric 
theory of diffraction and, independently, by asymptotically expanding the 
solution of the appropriate boundary value problem. They are in a form which 
is suitable for numerical evaluation. 

As an application of our results we shall now calculate the backscattered 
field from a perfectly conducting spheroid when a plane electromagnetic 
wave is incident along the axis. To do this we could apply the geometric 
theory of diffraction directly. However, it is simpler to utilize our previous 
results by means of the theorem proved in Appendix II. According to that 


theorem 
(1) E= (Ug — Up) E,. 


Here E is the backscattered electromagnetic field due to the incident field 
E,e—“*. The scalar backscattered fields are ug and uy for the soft and hard 
spheroids respectively, with the field e~“? incident. Upon using (15) of 


Section 2 for ug and uy, (1) yields 


_ —E,hcosh*a ( a (kh)"*(sinh a)** 
(2) E= papa rg exp{ik(x—2h cosh a)} \ 1— 6" (cosh a)! aye 


X exp{2zkh cosh a+in/3} [3-4 - exp} ini (kh) (cosh a sinh a)*/* 
3n/2 \ 
x J (cosh’ a—cos* )~"”* dn ¢+g2'A~* exp \ir(th) _ 


/2 
3n/2 " 
X (cosh a sinh a)”” f (cosh? a—cos” 9)” ant |) ’ 


w/2 


Here A and gq are given by (12) of Section 2, while A’ and q; are given by 
(13) of that section. The last term in the brace can usually be neglected. 
From (2), the electromagnetic backscattering cross section can be computed. 
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APPENDIX I 
ASYMPTOTIC EXPANSIONS OF THE SPHEROIDAL FUNCTIONS 


In this appendix we study the asymptotic expansions as kh + © of the 
various solutions of equations (4) and (5) in Section 3. It will be convenient 
to denote kh by k. To analyze equation (4) of Section 3, we set 


(1) #(¢) = (sinh &)-!/1(E). 
Then we find that [ satisfies the equation 
(2) I’ —k?(6?—cosh? £+2-!k-?—4—!k-? coth? £)T = 0. 


Since — > a > 0, we may apply the methods of Olver (1954) to obtain uniform 
asymptotic expansions of I'(é). The calculation is similar to the one carried 
out by Levy (1958). It yields in the present case for the outgoing solution, 
V®(£), of equation (4) of Section 38, 


(3) VO (€) ~es 31/3q-1 (52 — cosh? £)-!/4(sinh £)-1/24 (31/3e— te /3p2/3e) | 
In (3) ¢ is defined by 


g 
(4) 2 a = i‘ (b°—cosh? x) "dx. 
3 cosh! 


The Airy function 4 (¢) is defined by 


(5) At) = J “cos (2° —tz)dz. 


In the case when £ > cosh~! 6 we can expand the Airy functions in (3) to 
obtain 


(6) v(e) pa erie k!89-*-"? (cosh? t—b*)~"*(sinh ao 
g 
x exp ik f ,, (cosh e— 0) ax} 
cosh” *b 
The calculation of the eigenvalues, },, is exactly the same as in Levy (1958) 
and yields 


(7) b, ~ cosh a+7,k-?/3(sinh a)?/3(cosh a)~!/3, 
Here 7, is defined by 
(8) Tr = G7'8etr/3q,,, 


In (7) and (8) b:, tz, and g, are to have superscripts 1 for the soft spheroid 
and superscripts 2 for the hard spheroid where g‘” is defined by A (gq) = 0 
and g is defined by A’(g®) = 0. 

To obtain the asymptotic expansion of W{”(n) we must proceed in a 
manner different from that which we followed in the case of the corresponding 
elliptic cylinder function. This is due to the singularities of equation (5) of 
Section 3 at » = 0 and n = z. Accordingly, we apply a boundary layer type 
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of analysis to obtain the asymptotic expansion. However, we do not find it 
necessary to introduce a stretched independent variable because we consider 
only the leading term in the expansion. The use of such a variable would 
permit us to obtain further terms. In the neighborhood of 7 = 0, equation 
(5) of Section 3 becomes 


1 ° 
(9) er: Y'+p'y = 0, 
where 
(10) p? = k?(b?—1). 


Thus the solution regular at » = 0 is proportional to Jo(py). We may normalize 
it to have the value 1 at 7 = 0. Then 


(11) v ~ Jo(pn) for 7» <1. 
Similarly, in the neighborhood of » = 7, equation (5) of Section 3 has the 
form (9) with » replaced by r—n. Hence, we may conclude that 

(12) ¥ ~AJolp(r—n)} +BNolp(x—7n)} form—n <1. 


In (12) A and B are as yet undetermined. Away from 7 = 0 and 7 = 7, we 
may apply standard W.K.B. arguments to obtain 


, ( n. 
C cos {ef (b°—cos” 0)""adnh 4D sin iL f (6° —cos” 2)'"an 
- ala eae ae 


his... NY 
'““(b°—cos’ 9)" 


(18) ¥~ (sin 7) 


The constants A, B, C, and D are obtained by matching (11) and (12) to 
(13) in their respective common regions of validity. To obtain C and D we 
set » = ek~!/*. Then the argument of Jo in (11) is large as R—@ and we 
may expand Jo asymptotically for large argument. Upon doing this and 
requiring that the resulting expression be identical with (13) we obtain 


(14) C = D = (xk), 


To obtain A and B we set 7 = r—ek~'/*, Then upon expanding (12) asymp- 
totically and comparing with (13) we obtain 


sin {ef (b?—cos” n)"*an} ; 
0 


(15) A= 
and 
(16) B= ~cos {& f (6'—cos" 9) aa 


Using (14) in (13) we obtain the following result, which is valid except 
in the neighborhoods of 7 = 0 and 7 = =: 


(17) W(n) ~ (2/mk sin n)/?(b’—cos? 9)" cos} J "teauil n)"én—a/4} 











§ 
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A simple calculation, using (12) and (16), shows that 


(18) lim sin n ¥(4—7) = 29 "cos ye f (b°—cos” n)""én¥ 
0 


9-0 


We note that for 6 = b,, ¥(n) = WO? (n). 

In order to expand equation (17) of Section 3 asymptotically, we evaluate 
the expansions contained in this appendix at the eigenvalues 6,. When (7) 
is substituted into (17), we find that 


(19) W2?(n) ~ (2/ mk sin 7)? (cosh? a—cos’ 9) "4 


wn 
x cos} J (cosh? a—cos” n)"dn+k*/*r, (sinh a cosh a)*? 


" 
x J (cosh? a—cos* ay ay «/4} 
Similarly, the substitution of (7) into (6) yields 


(20) Vir (Ee) wee R92"? (cosh? ¢—cosh’ a)"* (sinh ¢)”” 


{ g 
X exp \ik J (cosh’ ——cosh’ a) "dt —ik"/*z, (sinh a cosh a)*/* 


a 


t 
x f (cosh” ¢—cosh’ a) "ag ‘ 


Finally, we shall obtain the expansion of 


ant 
To do so, we note that (3) may be differentiated with respect to both ¢ and 


6 to yield the results 


(22) Va) ~ 9 'g43?9 (27,) 42" (sinh a)~2*(cosh a)" "24" (qq) a 





(23) - = VP (a) ~ we h83"8 (274) 4R* (sinh a)" (cosh a) 77 *"/3.4 (qn) Shs 


Here ¢, is the value of ¢ at § = a, 6 = b,. Upon differentiating (4), we find 


(24) = ~ = $5"? (rq) !2-"*(6inh @)*(cosh a)". 
We also make use of the following identities: 

(25) 27M = tH, 

sas Ie _ Ay Ot 


0b «At, 0b * 





3 
y 
4 
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From (25) and (7), we find that 
Wn 1/37, —2/3 1/3 
(27) a 2°"(sinh a) “(cosh a)". 


Upon making use of (27), (22), and (23), we obtain the desired asymptotic 
expansion 


-1 
(28) vi" s va} tn me 3-49-38 (cosh ss 


X (sinh a)**{A'(qp)}~°. 


APPENDIX II 
A THEOREM ON BACKSCATTERING 


Let E be the electromagnetic field scattered back from a perfectly con- 
ducting body of revolution. The incident field E,e~* is assumed to be a 
plane wave incident along the z axis, which is the symmetry axis of the body. 
We shall show that E can be expressed simply in terms of ug and uy, the 
scalar fields scattered back from the same body with the soft and hard boun- 
dary conditions respectively. In these cases the incident field is the plane wave 
e~***. The result which we shall deduce is the backscattering theorem 


(1) E = }(ug—uy)E,. 


Our proof of this relation is based upon the geometrical theory of diffraction. 
Therefore its range of validity is the same as that of this theory. 

Let us first consider the field on the reflected ray. In both scalar cases 
and in the electromagnetic case this field is a product of four factors. They 
are the value of the incident field at the point of reflection, a reflection co- 
efficient, a geometrical or amplitude factor, and a phase factor. The latter 
two factors are the same in all three cases. The reflection coefficient is —1 
in the soft case, +1 in the hard case, and —1 in the electromagnetic case. The 
value of this last coefficient depends upon the fact that the incident field is 
tangential to the surface at the point of reflection because the incident ray 
is normal there. From the preceding statement and the relation —1 = 3{—1 
—(+1)}, we see that (1) holds for the reflected field. 

To prove that (1) also holds for the diffracted field, we first treat the case 
of a smooth body. According to Levy and Keller (1959), the diffracted field 
E, at any point P off the axis of the body is given by 


(2) E,(P) = x { (E,- do) dé v,(P) + (E,-d4o) diu,(P)}. 


The unit vectors dj, and dj are normal to the body at the points where the 
jth diffracted ray through P is initiated, and where it leaves the surface, 
respectively. The unit vectors d3, and dj are tangent to the body and normal 
to the surface ray at these same points. The scalar functions v; and v2 are 
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constructed according to the geometric theory of diffraction with the con- 
stants appropriate to a hard body, while 4, and wu: contain those for a soft 
body. 

As the point P recedes from the body in a direction parallel to the axis, 
the points at which the rays through P leave the body tend to limiting 
positions. The points of initiation of these rays also tend to limiting positions. 
All these limiting positions lie on that cross section of the body which is 
normal to the axis and which has the largest radius. Let us call the plane of 
this section the xy plane. It is shown in Fig. 3. We have chosen the y axis 





Fic. 3. A cross section of the ellipsoid through its center normal to its axis of revolution. 
An incident ray tangent to the surface at A travels half way around the ellipsoid to B, from 
which point it sheds a ray which returns to the source. The vectors d}, and dj, are respectively 
normal and tangential to the surface at A, and both are normal to the incident ray. The vectors 
d3, and d3, are associated with the ray incident at B, which shedsa ray back to the source after 
travelling around to A. This figure is relevant when the source is far from the ellipsoid. 


to be parallel to E;. The angle between the xz plane and the plane which 
contains P and the symmetry axis we denote by ¢. Then the limiting positions 
of the vectors in (2) are 


(3) dio = —d3) = —d} = d} = (cos ¢, sin ¢, 0), 
(4) dio = —dj) = d} = —dj = (—sin ¢, cos ¢, 0). 
When (3) and (4) are inserted into (2), it becomes 


(5) Ea(P) - — (E; -d}o)d20(v1-+02) + (Ey -d30)d30(w1+u2) 


a { — (uituetvi+02) sin 29, ui+U2—01—V2+ (Ui +u2+01+02) 
X cos 2¢, 0}. 


As P tends to the axis of revolution, the expression for E(P) on the right 
side of (5) becomes infinite. This is because the axis is a caustic of the diffracted 
rays. In Keller (1957) there is a method for modifying fields which become 
infinite on an axial caustic. It applies to scalar wave functions which are 
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proportional to sin m¢ or cos n@. Therefore it can be applied to each component 
of (5), provided the y component is separated into a term proportional to 
cos 2¢ and another independent of ¢. To apply the method a term in sin n@ 
or cos v@ is first multiplied by the factor 


(6) Cy, = 3(2mrkp cos 5)!/? sec{ kp cos 5— (n+4)x/2}J, (Rp cos 5). 


In (6) p is the distance of P from the axis and 7/2—6 is the angle which the 
rays make with the axis. The correct finite value of the field on the axis is 
obtained by letting p tend to zero after multiplication by C,. The result is 
finite because the factor (p cos 6)!/? in C, cancels the singular factors in 1, 
Uo, V1, and vo. It follows that on the axis the terms with ” ¥ 0 vanish because 
J,(0) = 0 for n # 0. Thus only the terms with » = 0 in (5) remain when 
P is on the axis, and (5) then yields 

(7) E,(P) A \o, lim Co(ui+u2—v1—22), of = a {0, us—us, 0} 


p+ 0 


= 4(u$—un) E,. 
The final form of (7) is the desired result (1) for the diffracted field. In obtain- 
ing it we made use of the fact that uf and uf, the backscattered diffracted 
fields for soft and hard bodies, are given by 


(8) ug = lim Co(ui+u2), 
p+0 

(9) us = lim Co(v;+02). 
p+ 0 


Numerical calculations of the backscattering cross section of a prolate 
spheroid have been carried out by Siegel et al. (1956). However, these calcula- 
tions were limited to small values of kh and hence are not directly comparable 
with the values given by (7). 

The same result (7) holds when the body has an edge. The proof is similar 
to that given above so we shall just sketch it. For simplicity we shall confine 
ourselves to the case of a circular disk. If Fig. 3 now represents the disk, then 
the field at P is again given by (5). Of course the functions 1, 2, v1, and v2 
are now those appropriate to scalar diffraction by a disk. The negative sign 
in the first term on the right in (5) arises because the backscattered ray travels 
in the direction opposite that of the incident ray. Since v; and v2 describe the 
diffraction of the tangential component of magnetic field H, the opposite 
relations between E and H on the two rays introduce the negative sign. From 
(5) we obtain (7) as before, and thus (1) is proved. 

We may also prove (7) for a body with an edge like that of a wedge, at 
which the two tangent planes meet at a non-zero angle. For that purpose 
the wedge diffraction coefficients given in equation (A10) of Keller's work 
(1957) must be used. The right side of that equation is missing a factor 
— (2/n) sin (x/n) due to a typographical error. Also in Keller’s paper (1957) 
the second columns of the diffraction matrices in equations (Al2) and (A13) 
must be multiplied by —1. 
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NOTES 


A COMPLEX NON-ANALYTIC MANIFOLD AND 
CONFORMAL MINKOWSKI SPACE-TIME 


D. K. SEN 


COMPLEX MANIFOLDS IN GENERAL 
A real 2n-dimensional manifold V2, is said to admit a complex structure if 
it is possible to introduce co-ordinate neighborhoods in the manifold with 
complex co-ordinates of dimension n. 
Let the co-ordinates of Vo, be (x') and let there be two co-ordinate neigh- 
borhoods connected by the co-ordinate transformation 


vs if 
(1) eg seg (er) x= x(x’) 


with a non-vanishing Jacobian, where the Latin indices take the values 
Beye <i ae oe et 


te 


(2 Ja” = x*+ ix" 
lz 


with its complex conjugate e = xix" 
where the Greek indices take the values 1, 2,...m. There is then a one-one 
correspondence, (2%, 2) = (x'), and (2) may be considered as the co-ordinates 


fi 


of the associated n-dimensional complex manifold. Equation (1) may then be 


writtel 
(3 g = sg (z,2); z = 3 (s,2). 
unplex manifolds which have so far been studied by the mathema- 
re complex analytic ones which are defined as follows. 
Let .. and N’ be two complex co-ordinate neighborhoods and P be a point 
such that PeN (\ N’. If the co-ordinates 2’ of PeN’ are complex analytic 


functions of s*, the co-ordinates of PeN, with a non-vanishing Jacobian, 
then the manifold is called a complex analytic manifold of real dimension 2n 
and complex dimension n. In case of an analytic manifold equation (3) can 
therefore be written as 

(4) z = y"(z), 2 = (2) 

where Y% denotes the complex conjugate of the function ¥(z). Putting 2* = 2 
we can write (4) as 

~ ° ° k 

(5) gs = f(s’) 


whose Jacobian is 


(6) J = = real > 0. 





dz | 
a 
dz | 
Such a transformation is known as pseudoconformal. 


Can. J. Phys. Vol. 38 (1960) 


145 











146 CANADIAN JOURNAL OF PHYSICS. VOL. 38, 1960 


All the usual concepts of tensor-calculus, such as vectors, tensors, affine 
connexion, and metric, can be defined in a complex manifold in the same way 
as in the real case. 

A particular complex analytic manifold which has been studied in great 
detail is the Kaehler manifold (1933). A Kaehler manifold is a complex analytic 
manifold with a positive definite metric 


(7) ds? = gydz‘dz*, 
where the gx are self-adjoint and satisfy the condition 


(8) Lan’ = 255° = 0. 


The above condition is invariant under the pseudoconformal transformation 
(4). The Kaehler metric therefore reduces to the Hermetian form 


(9) { ds” = 2 gagd2°dz", 
Lab 


where = £3a = Bab = Ba: 


CONFORMAL MINKOWSKI SPACE-TIME 


We consider a two-dimensional space-time manifold with x! = x and x? = 
ct and we introduce 
1 + 2 
=x+ix 
ae 
=x —1x 


(10) 1 


uy 


for the co-ordinate of the one-dimensional complex manifold. A proper Lorentz 
transformation of x! and x? 


Jxt = y(x'— Bx’) 


, = = (1— R*)\-3 
(11) Le a v(x? — Bx’) B =v/c,y = (1—6') 


corresponds to a non-analytic transformation of z 


(12) 
ly = (z-#,) 
4 
78 
. : : Y ‘se. : ; 
The transformation matrix 8 is unitary with the determinant 
7 
ine y 
1 
7B 
Y 7 = : 
J= =y (1-6) =1. 
78 
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The improper transformations are given by 


(i) space reflection: 2’ = —2, 2’ = —z,J =1; 
(ii) time reflection: 2’ = 2,2’ = 2, J = —1; 
(iii) space-time inversion: 2’ = —2z, 2’ = —2,J = +1. 


Space-time inversion is rather distinct in the sense that it corresponds to an 


analytic transformation of z. 
A vector (é, €) transforms like (z, 2) and a tensor of rank two will transform 


as 
(13) Eve = AvArein 


where Alis the transformation matrix of (12). 
One can introduce a metric in the complex manifold by means of a sym- 
metric tensor gy which is in general complex: 


(14) ds? = gudz‘dz*. 

Since z‘ takes only two values z and 2, we denote the corresponding metric 
tensors by g,z, g:z, and giz, so that (14) can be written as 

(15) ds? = g,,d2*+gz77d2"+2g,3d2dz 

£22(dx? — c*dt? + i2cdxdt) + gzz(dx? —c*dt? —i2cdxdt) + 2g,;(dx?+c*dt*), 


The metric reduces to the conformal Minkowski form 

ds? = 2g,,(dx?—cdt*), £22 real 
if 
(16) 8ez = Sia 82 = 0. 


However, one must show that the above condition is an invariant relation 
under Lorentz transformations. This can be seen from the transformation 
formula of the metric tensor. We have from (13) 


2 
v'twt(24) en-+2y( 24) Beis 
2 
(18) a7 = (-28) tuty'en+2( -28) VE:i) 


( - x8) Zeet 4) veut2y'¢.2. 


t 


(17) Za's’ 


(19) Z2'2' 
From (17), (18), and (19) it is evident that, if 


B22 = B22 gas = 0, 


then also 


Bee = Br, Be = 0. 
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Thus the condition (16) is Lorentz-invariant. In fact, the only condition 
that g.; be zero in all Lorentz frames is sufficient. Because from (19) we then 
have 

2 


0 = vf zz zz 
; (g22— zz) 


or 


And it follows therefore from (17) and (18) that 


Sea's’ = 823": 


So a sufficient condition for the metric (15) to reduce to the conformal 
Minkowski form is the vanishing of g,z in all Lorentz frames. It is interesting 
to note the difference between (16) and (8), the condition for a Kaehler 
manifold. 

In general, however, the manifold would be characterized by a non-vanishing 
affine connexion [!(gm,) and curvature tensor Kf, 
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THE NUCLEAR MATRIX ELEMENT RATIO IN THE 0° — 0° 
BETA TRANSITION OF Pr'* 


J. M. PEARSON* 


On the currently accepted view of the beta-decay interaction having the 
VA form, the 0-0, YES, transition is of considerable interest in that the 
spectrum shape is determined completely by a single nuclear matrix element 
ratio, \ = 2 (y5)/(é.r). Such a transition is believed to take place in the 
ground state decay of Pr'#4, which recently has been well measured by Graham 
et al. (1958) and by Porter and Day (1959). Both have found A to be positive 
with the latter quoting \ = 542. Porter and Day point out that this is in 
clear-cut conflict with the estimate of Ahrens and Feenberg (1952), whose 
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method gives \ ~~ —12. The purpose of this note is to make the following 
two observations: 

(i) Rose and Osborn (1954a) have claimed that the non-relativistic limit 
(ys) ~ —(1/M)(6.p), used by Ahrens and Feenberg, has an incorrect sign, 
whence A ~~ 12. 

(ii) Ahrens and Feenberg base their evaluation of the nuclear matrix elements 
(in their non-relativistic limit) on the semiempirical mass formula. The author 
has made an essentially different calculation in which a single-particle j-7 
coupling model is assumed. The same non-relativistic approximation for 
(ys) is made, using the sign of Rose and Osborn and following the methods 
of Rose and Osborn (19546). It is found that ’ = 2.5, if the single particle 
potential is taken to be an infinitely deep square well. With harmonic oscillator 
potentials \ will still be positive, provided that the well strength of either 
parent or daughter does not exceed the other by a factor of 13/9. If the two 
well strengths are equal, then A = 8. The two different methods of calculating 
the matrix elements are thus in reasonable accord as to absolute value. 

The question of agreement or non-agreement between theory and experiment 
therefore depends essentially on the sign of the non-relativistic limit of (ys): 
Rose and Osborn’s form is clearly favored by experiment. 
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LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 


previous to the first day of the month of issue. No proof will be sent to the authors. 


Frequency Measurement of Standard Frequency Transmissions": ? 


Measurements are made at Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
as reference standard (with an assumed frequency of 9 192 631 770c.p.s.). Frequency deviations 
from nominal are quoted in parts per 10. A negative sign indicates that the frequency is 


below nominal. 





GBR, 16 kc/s 
Date, — 
October MSF, 7-hour 24-hour KK2XEI, 
1959 60 kc/s average* average 60 kc/s 
1 — 164 -171 — 164 N.M. 
2 —170 — 167 — 168 N.M. 
3 —176 -171 -—177 N.M. 
4 —176 -177 -177 N.M. 
5 —173 — 168 — 169 — 69 
6 —178 — 166 — 168 —63 
7 —170 — 167 -—172 -77 
8 — 169 -171 — 167 —73 
9 — 168 -—170 —-174 —70 
10 —168 —172 —176 N.M. 
11 —170 -171 —176 N.M. 
12 -—170 -173 —170 —72 
13 -171 —175 —174 —86 
14 —170 —170 — 168 —76 
15 — 167 -- 168 — 167 —69 
16 — 164 -171 — 168 —83 
17 — 169 —170 — 168 N.M. 
18 — 168 —172 — 168 N.M. 
19 —170 -171 — 167 —74 
20 — 163 — 166 — 164 —76 
21 — 169 — 163 — 164 —76 
22 — 165 — 167 — 165 —75 
23 — 168 — 165 — 165 -77 
24 — 158 — 166 — 164 N.M. 
25 — 167 — 162 — 164 N.M. 
26 — 162 —165 N.M. —75 
27 — 168 N.M. N.M. N.M. 
28 — 166 — 164 — 162 —74 
29 —161 — 159 — 158 N.M. 
30 N.M. — 162 — 163 N.M. 
31 — 169 — 164 — 163 N.M. 
Midmonthly 
mean — 168 — 168 — 168 —74 
Midmonthly 
mean of WWV — 104 


Note: N.M. No measurement. 


*Time of observations: 00.00 to 05.30 and 22.30 to 24.00 U.T. 
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“Spin-flip Narrowing’’ of Paramagnetic Resonance Lines 


Copper, cobalt, and zinc potassium sulphates are isomorphous, and mixed crystals can be 
made with any desired proportions of copper, cobalt, and zinc. The paramagnetic resonance 
spectrum of the copper ions consists of four lines, with an over-all separation of some 250 
gauss. These lines are broadened by spin-spin interaction so that the four lines are unresolved 
at all temperatures in copper potassium sulphate, and are further broadened at high tempera- 
tures by spin-lattice relaxation. In order to resolve the four lines, it is customary to reduce the 
spin-spin interaction by diluting the copper with diamagnetic zinc, dilutions of 1 part copper 
to 1000 parts zinc being usual. The paramagnetic resonance spectrum of cobalt is similar to 
that of copper, except that the spin-lattice relaxation times are much shorter. As a consequence 
of this, only a very broad line can be seen at liquid air temperature, and no line is visible at all 
at room temperature. 

We have observed the paramagnetic resonance spectrum of copper at 3-cm wavelength at 
room temperature, and at liquid nitrogen temperature (77° K) for the following two compounds: 
(5% Cu 95% Zn) K2(SO4)2.6H2O and (5% Cu 95% Co) K2(SO,)2.6H20. In the case of the 
copper-zinc specimen, a broad line, 300 gauss wide, is seen at room temperature, and at 77° K 
this line narrows somewhat and the four hyperfine components start to be resolved. For the 
copper-—cobalt specimen, we find a line 300 gauss wide at room temperature, just as in the case 
of the copper-zinc specimen. On cooling the copper—cobalt specimen to 77° K, the copper line 
almost disappears, becoming very feeble and at least 2000 gauss wide. 

The behavicr of the copper-zinc specimen is easily understcod; the line broadens at room 
temperature due to the shorter spin-lattice relaxation time of the copper ions. For the ccpper- 
cobalt specimen we suggest the following explanation. At 77° K, the copper line is broadened 
by the spin-spin interaction between the copper and the cobalt ions. At room temperature, 
the cobalt ions have such a strong interaction with the lattice that they remain in any one 
spin state only for a very short time. The interaction magnetic field, produced by the cobalt 
ions, is thus fluctuating so rapidly that the copper ions see only its average value, which is 
zero. A somewhat more quantitative explanation of this effect fcllows from the treatment of 
relaxation by Pines and Slichter (1955). According to these authors, the contribution of the 
cobalt-copper interaction to the width of the copper line should be 6H or g8(6H)*r/h according 
as t > or < g86H/h (equations 6 and 4 respectively of their paper). In these formulae 6H is 
the average interaction field produced by the cobalt, g is the g factor for copper, 8 is the Bohr 
magneton, and 7 is the correlation time for the cobalt, that is the average time for which a 
cobalt ion stays in any one spin state. There is a close relation between 7 and the spin-lattice 
relaxation time 7). (For instance in a two-level systém 7; ~ 1/2. In a calculation as crude as 
this, the difference between the 7's of the two states can be ignored. It should be noted, how- 
ever, that 7 depends on exchange as well as spin-lattice relaxation.) 

We have thus demonstrated that in the case of paramagnetic resonance of compounds con- 
taining two paramagnetic ions with very different values of the spin-lattice relaxation time, 
the width of the line of one component decreases when the spin-lattice relaxation time of the 
other becomes very short, a situation analogous to motional narrowing in liquids. The study 
of the line breadth as a function of temperature should enable the spin-lattice relaxation time 
of the other component to be measured. Further, this should permit the measurement of spin— 
lattice relaxation times of an order of magnitude smaller than has been possible with other 
techniques. It is possible that there would be a resonant effect on the signa] strength when the 
correlation time of the other component is equal to the Larmor precession time of the first. 
These investigations are being continued. 

The apparatus which was used to measure these lines has some novel features. It is basically 
a 3-cm transmission cavity spectrometer with 60-cycle modulation and crystal—video display. 
With such a system, 60-cycle pickup can be very serious, and the 60- and 120-cycle com- 
ponents of the output from the crystal are usually eliminated at an early stage using two twin 
T filters. Such a system does not produce a faithful oscilloscope display of a broad line, since 
such a line has important components at low frequencies. (We may define a broad line as one 
whose important Fourier components in the video circuit are all below 1 kc.) To overcome this 
difficulty, we adopted the following system. A 15-kc square wave was applied to the reflector 
of the klystron to switch it on and off. The output of the crystal detector was passed through 
the usual twin T’s and through a filter which passed only frequencies between 10 ke and 20 ke, 
and whose phase characteristic was linear between 14 kc and 16 ke. A simple Fourier analysis 
of this process shows that the output of this circuit is a 15-kc sine wave carrying the paramag- 
netic resonance absorption signal as an amplitude modulation. This output is then added to a 
15-kc reference voltage tc reduce its amplitude, amplified more, detected, and finally dis- 
pare’ on the oscilloscope. This system enables us to filter out the pickup without distorting the 
ine shape. There is also, in theory, a slight enhancement of the signal/noise ratio. (A system 
similar to this using 1-Mc frequency modulation of the klystron and 2-Mc detection was used 
by Bagguley and Griffiths (1952) to enhance signal/noise ratio.) 
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